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The Electronic Structure and Bond Lengths of 
Coronene and Pyrene 


By W. E. MOFFITT anp C. A. COULSON * 


Physical Chemistry Laboratory, Oxford 
* Now at King’s College, London 


MS. received 1 September 1947 


ABSTRACT. Molecular orbital and valence-bond pairing treatments are given for the 
bond lengths of coronene and pyrene. In both cases the agreement of the molecular orbital 
treatment with the experimental values is good, and in both cases it is rather better than the 
agreement of the valence-bond method, though certain characteristic features are common to 
both treatments and are found experimentally. The Fries rule, which governs the relative 
importance of the various Kekulé structures in valence-bond theory, is shown to be only 
partially applicable. A more valid rule is suggested for condensed hydrocarbons generally. 


Sa LINER O DUCTION 


on the dimensions of the polynuclear hydrocarbons coronene and pyrene 

is of great interest. It discloses a new testing ground for the various 
theories of the electronic structure of these aromatic systems. In this paper 
we give the results of two alternative theoretical treatments of these molecules. 
In the first, using the molecular orbital method of Hund, Mulliken, Hiickel 
and Lennard-Jones, we express the electronic wave function as a product of 
non-localized one-electron molecular orbitals (called MO), constructed as linear 
combinations of atomic orbitals (LCAQO). In the second (HLSP) method, 
due to Heitler, London, Slater and Pauling, the total wave function is obtained 
as a linear combination of many-electron functions or structures, in each of 
which the electrons are perfectly paired; since these latter calculations would 
be extremely tedious if followed to the full, we examine an approximation to 


them. 


TT recent very accurate x-ray work of Robertson and White (1945, 1947) 


§2. THE MOLECULAR ORBITAL METHOD 


The characteristic properties of aromatic systems are attributed to those 
electrons in orbitals which are antisymmetrical with respect to reflections in the 
planes of these molecules. We shall consider the behaviour of these (7) un- 
saturation, or mobile, electrons alone, and make the usual assumption that their 
interaction with the remaining o electrons is negligible for our purposes. The 
LCAO MO procedure has been fully described elsewhere (e.g. Lennard-Jones 
and Coulson 1939), so that we shall merely quote the results of our calculations 
for pyrene. The treatment for coronene 1s of greater interest, however, owing 


to the high symmetry of the molecule. 


PROC. PHYS. SOC. LX, 4 aT 


310 W. E. Moffitt and C. A. Coulson 


The symmetry elements of 
coronene, like those of benzene, 
belong to the point group Dg). 
The 24 2p, atomic orbitals, one for 
each carbon atom (figure 1), which 
we call *the 6(@—[02 aca), 
constitute a reducible representa- 
tion I’ of this point group. The 
characters of I’ are found to be 


x(E) =24, x(C2) = x(Cs) = x(Ce) 

= (Ca) 0 x0Ca)) eee 
x(n) = —24, x2) = x(S6) = x(S3) 

=(o,~)=0, x(a, ) = —+- 
The notation is fully explained by 
Wilson (1934). Accordingly, by 
the theory of group characters, [ 
has irreducible components 

T= Aja 3Aeq + Big -3 Boy 4h (1) 

Following a method analogous to that used by Bloch (1929) in his discussion 
of metallic structure, we define the functions 


K(w) =$, + wd, + wd; + wd, + wtds + ode, | 
\(w) = by a wg i wby Rs w*h19 a: why at: wo, 
p(w) = b43 + why5 + wd 17 + w*d49 + w4ds, +. w°do3, f 


V(w) = bt Odig + O*Prg + who twhpy, + wo. J 
The numbering system of the ¢; is as in figure 1, and w is a number, in general 
complex, of unit modulus. It may now readily be shown that the required 
LCAO forms for the molecular orbitals are = 
Ay: (1) — (1), 
Agu: hx(1)+IA(1) +m{u(1) + r(1)}, 
1 M—D+o(=0), 
set Ha =) 4EM +m, fT % 
> R''k(e) + 1X) + m'"wl(e) +n" v(€) ‘ 
yy 2 mr tr yr Se eu OF ee 
au: RM K(e?) 41 Ae?) +m!" (e2) +n!" v(e2), J 
In these formulae k, J, m, k’, ..., n'’ are constants, one set being found from |} 
each solution of the secular equation. These sets of coefficients are easily 
calculated, for, according to (1), the secular equation is a cubic for the levels of _ 
symmetry A,, and Bog, it is linear for A,, and Biz, and a biquadratic for the ] 
degenerate levels Ej,, E,,. The secular equations may be written down when | 
the matrix components of #—E in the representation (2) are known, where J] 
H is the Hamiltonian and £ is the energy, as usual. We have | 
(|? —E| x) =6xB + 6(e0 +08, 
(«| 9 —E|)) =A —B| py) = (lo? —E |») =68, 
(XH —B|\) = (u| 0? —E |p) = (10 —E| v) = 6x8, 
(10 —B| v) = (v1? —E! p)* = 608, 


Figure 1. Coronene C2,Hjp. 
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where f is the resonance integral between adjacent 2p, orbitals, E, is the atomic 
energy term (Coulomb term) defined by E)=(¢,|¥%| ¢;), and «=(E£,—£)/p. 
‘The secular equations thus reduce to 


Ain: 0=x—1, 
ASL eae 0 
AS = 1 4 be za\ed 
0 2 2x +2 
B,, 0O=x-+1, 
Xe 0 
Dic = 1 oy z. = —N(— x); 
0) ey, 2x —2 
x+1 1 0 0 
: 1 3 1 1 
ie 10 ay ‘ : =—}(x+ 1)A(—x- 1), 
0 1 Coens, 
ve teil 0 0 
< 1 x 1 1 
Fou: US 0 1 . 2 ee ING aly) 
0 1 2a 
‘Thus the roots are all either +1 or simply ee ee ae 


expressible in terms of the roots of the cubic 


equation for the levels A,,. The resultant a eee 
spectrum of energy levels may be represented ae A 
graphically as in the accompanying diagram Tors 
(figure 2). In the ground state we fill the Sins 
. . aur 
orbitals of lowest energy with electrons, two 
- 0539 e 


exclusion principle. ‘These electrons, whose 
spins are anti-parallel in pairs, are represented 
by crosses on our diagram. The final electronic 
configuration may thus be written Deer fomte 


at a time, in accordance with the Pauli = 
— 


TA ig? (Gou)?(C1g)*(€2n)*(@ou)*(b2)*(B1g)*(ig)*(Cou)*- = = 1-539 acy 

The explicit LCAO forms for each MO 
are readily calculated at this stage, since the XXX— 2214 eg 
solutions of the secular equations enable us =a 2-675 any 
to determine the values of the sets of coefficients 
k, l, m, nin (3). That is, each MO may be 
written in the form 


Figure 2. Energies of molecular orbitals 
in coronene 


PO ev bn (7, 4= 15 2, 2225.24), 
where the c*” are known, and satisfy the normalizing condition &| c\2—1, The 


mobile (or 7) bond order p,; of a given bond between atoms 7 and j may 
now be calculated, for it is the sum of separate contributions from each of the 
various 7 electrons. Following Coulson (1938), the calculation 


( 
Pye Upp = UE Ccy”* + cL) Fey”) 
ue 


21-2 
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is readily accomplished. ‘The final values have been tabulated under column I 
for coronene (table 1) and pyrene (table 2). The internuclear distances which 
these would lead us to predict are given in column 4. Part, but not all, of the 
information in table 1 has already been given in a preliminary note (Coulson 
1944). 

§3. THE PAIRING METHOD 


The HLSP, valence-bond, or pairing, method of treating aromatic molecules 
consists in selecting that linear combination of the structures ys, which gives 
the minimum energy to the particular molecule under consideration. The 
procedure is well known and has been applied to many molecules by Pauling 
(1933) and others. The complete treatment of our molecules along these lines 
would be an extremely arduous task. For coronene we should have to take more 
than 200,000 non-ionic structures into account. ‘Therefore most of these 
structures are neglected, and only a reasonable number are considered to con- 
tribute in any marked extent to the ground states of these molecules. It is 
customary, especially in qualitative discussion, to make the simplification of 
neglecting all ionic and excited structures containing formal bonds (non-adjacent 
pairing). Accordingly we confine ourselves to the remaining non-excited, 
or Kekulé, structures (figure 3), which number 20 and 6 for coronene and pyrene 
respectively. 


bo Ho Yai 


wy 
Ve 72 


Figure 3. Canonical structures for coronene (top line) and pyrene (bottom line). 


By symmetry we may write the two wave functions of the ground states as: 
coronene ‘P= a(yy +2) + b(p3 +a) + c(h + ... +449) +d (by + P12 +45) 
+ ey t ... +39) +frboo, . 

pyrene = ‘P= r(hy + pp) + 5(h3 + by) + Abs + Yo). 
The secular equations are of degree 6 and 3 respectively. They may be set up 
with the help of the symbolic technique devised by Pauling (1933), and they 
are easily solved. For the two ground states 
(coronene) = Q + 6:64J, and a:b:c:d:e:f=0-50:0-60:0:76: 0-75: 0-88: 1-00, 
E(pyrene) =Q+4-70J, and) "57 :s 7 f= 153191272 1-00: 
This corresponds to resonance energies 3-64J and 2:20./. Q and J are the 
Coulomb and exchange integrals as usual. 
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The bond orders may now be defined in one of two alternative ways. The 
first is due to Pauling, Brockway and Beach (1935), which we shall abbreviate 
as (PBB). The second is due to Penney (1937), using a technique devised by 
Dirac: we shall abbreviate this as (PD). The latter definition of Piz as 
P= — (3) §;. s;, where s; and s, are the spin matrices of the electrons on atoms 
z and j, is fundamentally more satisfactory than the earlier one which is merely 
the weighted mean of double-bond character as revealed by the weights of the 
various structures. ‘The (PD) bond orders are shown in column 2 of tables 1 
and 2; the (PBB) bond orders are shown in column 3. In columns 5 and 6 
we give the (PD) and the (PBB) predictions of the bond distances (in a.), using 
appropriate curves of order against length in each case. The rows labelled 
p and r in the tables give the average bond order and internuclear distances 
respectively. 


‘Table 1. Coronene 


1 2 3 4 5 6 7 
MO Pairing MO Pairing es 
Coul Pe) 
(Coulson) (PD) (PBB) (Coulson) (2D) eee (eB) SRWinte) 
1 O75 0-816 0-742 1es72. 1-362 1-362 1-385 
Q 0-538 ose —  (Quneys: 1-411 1-450 =: 1-438 1-415 
R 0-538 0-629 0-484 1-411 ISSRS7/ iL asxoy 1-43 
Ss. 0-522 0-328. 770-2538 1-415 1-446 =1-438 1-43 
Dat, 0-576 0-481 0-400 1-405 1-419 1-414 1-415 
Table 2. Pyrene 
1 2 3 + 5 6 a 
MO Pairing MO Pairing Pee 
(Coulson) (PD) (PBB) (Coulson) (PD) (PBB) and Wie) 
A 0-536 0:280 0-231 1-412 1-458 1-445 1-45 
B 0-524 0:516 0-384 1-415 1-408 1-410 1:39 
c 0-594 0:587 0-500 1-401 1371390 1-42 
D 0-669 0:587 0-500 1-388 Heyy essere 1-39 
E 0-503 0-141 0-116 1-420 1-495 1-484 1-45 
pres Ona 0-931 0-884 1-370 1:347 1-348 1&39 
p,7 0-592 0-498 0-421 1-403 1-409 1-412 1-41 


Columns 1-3 are bond orders, columns 4-7 are bond lengths. 
‘The last decimal in columns 4, 5 and 6 for the bond lengths (A.) is unreliable on an 
absolute scale, but is useful for comparisons. 


§4. DISCUSSION 


Robertson and White (1945, 1947) have determined the dimensions of 
coronene and pyrene by means of x-ray crystallography, and we have reproduced 
their values in tables 1 and 2. First, we note that for coronene the agreement 
between the LCAO MO predictions and the empirical data is very close indeed. 
The average discrepancy is only 0-01 a., and in no case is it more than 0-02 a.— 
quite within the limits of accuracy of the x-ray determinations. The pairing 
method gives somewhat poorer results, although the average internuclear distance 
agrees very well with the measured value. But there are relatively large variations 
in the lengths of individual bonds, which are apparently not found experimentally. 
Thus both the x-ray and MO methods gives the bonds R and S a practically 
equal value, whereas the pairing method expects these to be separated by some 
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0-05 or 0:06 a. As Robertson and White have remarked, using only the non- 
ionic non-excited structures, the empirical data can be reconciled with this 
representation if the relative weights of the individual structures follow a scheme 
which is the exact opposite of that which we have calculated; the failure of this 
set of structures to distinguish between the bonds Q and S is also to be remarked. 

It has been suggested (the Fries rule) that those Kekulé structures which 
show the greatest number of benzenoid rings should have the greatest weight 
in the superposition diagram. In table 3 we show that this is roughly valid for 
coronene. Certainly the greatest weight is for yy), which has the greatest number 
(6) of benzenoid rings; and the least weight is for ¥, and ys, with the smallest 
number (1). But apart from these extremes, the weights do not support the rule, 
which should therefore only be granted a relative validity. A better rule seems 
to be that those structures which show the greatest number of double bonds in 
the ‘‘exposed”’ positions (e.g. P of coronene) have the greatest weight. Table 3 
shows that there are effectively no exceptions to this in coronene. It seems 
as if there was a tendency for such outer bonds, in which both ends of the bond 
are only attached to one other aromatic carbon, to be double; both phenanthrene 
and pyrene show the same characteristic. For pyrene, shown in table 4, the two 
types of structure with most exposed bonds (here only 2) have a similar weight, 
which is greater than that of the remaining structures with only one exposed 
bond. But the situation regarding the numbers of benzenoid rings is much 
the same as in coronene. 


Table 3. Coronene Table 4. Pyrene 
Structures Coefficients Np Ne Structures Coefficients iN pn ele 
ty, be 0-501 1 3 ty, fo 1:31 D, 2 
bs, tba 0-603 4 3 ts, wy 1:27 3 2 
es os o Why 0:757 3 a ds, be 1-00 1 1 
tu, tre, vis 0-749 Z 4h 
tis. -+ Pig 0-879 4 5 
boo 1-000 6 6 


Np=number of benzenoid rings ; Ne=number of exposed bonds. 


A similar comparison between theoretical and experimental bond lengths: 
may be made for pyrene as for coronene. We notice, however, that the difference 
between the MO and x-ray values is larger here than in coronene. This is 
possibly due to a lack of certainty in the x-ray results, since only nine of the sixteen 
carbons are “‘seen”’ in the diffraction pattern of the crystal. Again the LCAO MO 
method gives the more likely of the theoretical values. The high double bond 
character of the exposed bond F is, however, expected by all these methods—in 
agreement with the chemistry of this molecule. 

The relatively poorer values given by the HLSP treatment are readily 
attributed to the lack of a truly representative set of structures. In conjunction 


with the findings of Daudel and Pullman (1946) for a rather different set of 


polynuclear hydrocarbons, we may therefore conclude that in these condensed- 
ring molecules the contribution of excited structures to the ground states is. 
considerable, since the LCAO MO and the full HLSP calculations using only 
non-ionic structures generally lead to very similar results for those simpler 
systems where they can be carried through completely. 
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ABSTRACT. The variation of current with voltage at an idealized contact between two 
crystals of a semiconductor is calculated on the assumption that electrons penetrate the 
surface barrier by tunnel effect. Comparison with experimental curves for silicon carbide 
powders leads to the conclusion that these crystals are covered by an insulating surface film. 
The existence of such a film affects the rectifying properties of the crystal when it is in 
contact with a metal. It is shown that this can account for the discrepancies observed 
between experimental curves and those deduced from the simple Schottky theory of recti- 


fication. 


§1. INTRODUCTION 


HE non-ohmic resistance between a metal and a semiconductor can be due 
either to the existence of a contact potential between them (Schottky’s 
theory) or to the existence of an insulating layer of different composition 
(Wilson’s theory). In this paper we examine the effect of an insulating layer 
present in addition to the Schottky layer, the current penetrating the insulating 
layer by tunnel effect. This model seems correct for silicon carbide with thin 
quartz layers on the surface. Rectification is in general in the same direction as in 
‘the Schottky theory, but the rectifying properties are less good. We consider 
also the effect of these layers on the contacts between particles of a silicon carbide 


powder. 


§2. THE CONTACT BETWEEN TWO SEMICONDUCTORS 


The resistance of the contact between two crystals of silicon carbide is of the 
order of 101° times the resistance of the main bulk of the crystals for potential 
differences of less than 10-2 volt. As the voltage is increased, the resistance falls 
rapidly, as shown by figure 1, reproduced from the paper by Braun and Busch 
(1942). We conclude from this and earlier work (Kurtschatow et al. 1935) that 
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this is due to the existence of an insulating film on the surface, and show that 
it is impossible to explain the magnitude of the effect otherwise. 

Asa model for calculation, the contact is represented by two spheres, of radius 7, 
in point contact. The area of contact of the surfaces is then zero, but electrons 
can pass through the air-space in a narrow zone surrounding the point of contact. 
For a given value of the angle 6 between the common axis OX and the radius vector 
r, the width of the interspace measured in the direction of the common axis is 
(figure 2): 


x =2r(1—cos 6). adie ch) 
5 |- 
‘125 
— ut 
Oo 
as 
= 
<x 
“0S 
025+ 
in bas 
00 200 «300» 400 »=— $00 v/em. 
| 2 3 4 5 V/contact 
Figure 1. The current-voltage char- Figure 2. Idealized contact between two spheres. 
acteristic for silicon carbide crystal 
contacts. 


Figure 3. The potential energy of an electron Figure 4. The potential barrier in an 
in the interspace. external field. 


The potential V of an electron in the interspace is represented in figure 3. Ey is 
the activation energy of the semiconductor; thus the conductivity depends on 
temperature according to the formula 


o = oe Bol2RD, eas 


e¢ is the work function of the semiconductor; the fraction of incident electrons 
having sufficient energy to surmount the potential barrier is then 


ian oe IRB) 


Since ed is estimated to be of the order of 2-4 ev., this is negligible at room temper- 
atures (RT~2-5 x 10-*ev.). The chief mechanism by which electrons cross the 
interspace is, therefore, that of ‘‘tunnel effect’’. 

On this mechanism the current may be calculated as follows. Under an 
applied voltage V the potential barrier takes the form shown in figure 4. 
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The finite probability, P, that an electron with insufficient energy to surmount 
the barrier should penetrate it is given by 


P me. eoecee (4) 
where 
x= $nl($— V6] y= Amt, 


and h=2zh is Planck’s constant. 
The number of electrons incident on unit area of the barrier per second from 
left to right is 


N,=43NV(2kT/nm). ne (5) 


From right to left the number with energy great enough to enter the second 
conductor is . 


Na =a NAAR Dire Be 1h hae (6) 


where WN is the number of free electrons in unit volume of the semiconductor. 
The element of area, normal to the field in the angular range 0, <0 <0,+d0, is 


dA 272s, ay COS 0s.) 1 9 eae. (7) 
The current crossing this area is 
dI=e(N,—N,)PdA 
= N{/(RT/2nm)}(1 — e eV RP )e—* cos IV | Dar72 sin Oy C08 Oy d0....... (8) 


Integrating over @ the total current crossing the contact is 


rO= In 
T= ie ai=N{ V(R T/2nm)}(1 = e- VIPD) Dry f(ar | V), Mee (9) 


where 


V} (ar/V)® (ar/V)? © (ar/V)° 
If +>), where A1=2«¢"?, then f(«r/V)~V/ar, and the current is given approxi- 


mately by 
Te NGA (Rigi) (ler!) Arr Vio, | 2 aie eee (10) 


The integration over the surface from @=0 to 6=@, approximates closely to 
that over the whole hemisphere for values of 6, such that 


OF cos), 10a ee (11) 


This means that the emission of electrons is practically confined to the region 
0 <0 <6), and there is an effective area of contact 


INA 0 eter (12) 
where 
a® = 2r?(1 —cos 69) = 107A. 


For such a small area the classical resistance due to the constriction of lines of 
current flow will be considerable. ‘The contact may be represented by a cylinder 
of the semiconductor of area 7a? and length L, where 


L=(10rA/N)/(2nm RT) . (a/2m)\(1—e PAA, a, (13) 
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The total resistance of each contact is then 
R=(L4+Ba)/naa,) a ee (14) 


where B~0-82 (J. H. Jeans, Electricity and Magnetism). 
The specific resistance p of a mass of crystals of mean diameter 275, each im 
contact at a point of average curvature, will be given by 


p=2r,R =(27,/7a*o)(1 + 8a) = ae ee (15) 


Curve (a) in figure 6 is drawn in the case 
r=1, b=4-5 volt and c=0°'5 ohm cm ™ 
For V<kTvie this gives a resistance per 
contact of the order 10-2ohm, so that it is 
clear that p is of the correct order of mag- 
nitude only for p:py~10-8, which is not 
the case. Also the form of equation (10) 
shows a much slower increase of current 
with voltage than is actually observed. It 
therefore appears necessary, in order to 
obtain a sufficiently high resistance, to 
assume that the surfaces are separated by 
an. insulating film. The contact is now Figure 5. Contact between two semi- 
represented by figure 5. conductors through a surface film. 

The thickness of the barrier 


x=2r(l—cos0)-+-05,. 7 ener (16) 


where 4, is the thickness of the film on each surface. On substituting this value 
for x in equation (8), the current is given by 


T= (e/N)s/Qaum/RT) (1—e--"™ )\2arfler|V eae (17) 
In the approximate form this is 
I =(o/N)4/(2am|RT) (1 = e-°"!"7) (2arV ajo, sss aes (18), 
and the specific resistance is 
p=(2%)e)(Li-- Ba\ianas ee (19) 
where 
; L' =O” N(R 2am) (er Ae eee (20) 


Curves (6), (c) and (d) in figure 6 correspond to values of 59=1, 2 and 10a. 
respectively. ‘The resistance for V<kT/e is of the right order for 8)~10a. 
Figure 7 shows how the shape of the theoretical curve compares with experiment. 
The agreement is not good. _ Busch’s results (Braun and Busch 1942) for rounded 
grains agree well, but Mitchell’s results (Mitchell and Sillars 1947) differ widely 
from the theoretical curve. Since, however, he finds some rectification even 
between crystals of the same type, his contacts cannot be of the simple symmetrical 
form considered and agreement is not to be expected. He does not observe the 
inflection of the curve at V~kT/e due to the factor (1—e-°”/*") and ascribes. 
this to “‘leakage”’ through the barrier possibly due to gaps in the oxide layer. It 
seems more likely to be due to the rectifying properties already noted. 
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Evidence on the temperature-dependence of the resistance of the barrier is 
inconclusive. According to tunnel effect mechanism it should be the same as 
that of the specific resistance of the crystals. 

Electron diffraction experiments indicate that there probably exists a surface 
film of the thickness required. 


10°? 10? 107! 2 10' 10? 10° Volt /cm. 


10? io"! 10° 10* Volt/cm 
T 
10" vl , 
10"? 10'}- 
103 107 
a 103. 
€ -5 
a = 104- 
& i= 
10° So} 
= Mi(l) 
107 10° 
| 
ooh 4 107 
/ 
| Ae, : 
lo"? YB / 10°°|- i 
iy Bi!) / 
BS (ANS: / 
1 Busch’s E ( ( 0° M(3) ve 
Values eee 
10° 6(d) -7 
Figure 6. Silicon carbide powder. Theoretical Figure 7. Silicon carbide powder. Experi- 
curves. mental curves. 
50 contacts/em. Work function d=4:5 v. Braun and Busch: : 
Film thickness B (1) Region of validity of Ohm’s law. 
(a) = 0. (c) d=2a. B (2) Rounded grains. 
(Oy) C=, (Gh) Beal. B (3) Angular grains. 


Mitchell and Sillars: > 
M (1) Bonded disc ; particle diam. 0-005 cm. 
M (2) 45 3 45 LO OLS crm: 
M (3) i es ss EO z025.cme 
6 (2d) 6=10 4. ; d6=4-5 v. (50 contacts/cm.).. 


We conclude that the theory that conduction at these contacts takes place by 
tunnel effect through a surface barrier accounts for the magnitude, and approxi- 
mately for the variation with voltage, of the resistance of the contact. The problem 
to be considered next is that of the nature of a contact between such a crystal anda 


metal. 


§3. THE CONTACT BETWEEN A SEMICONDUCTOR AND A METAL 


In this section the conduction is always taken to be due to the presence of free 
electrons. The results can be applied to defect semiconductors, in which 
positive holes carry the current, by a change in sign of the electronic charge. 
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The resistance of the contact between a metal and a semiconductor depends not 
only on the magnitude but also on the direction of the applied voltage. The early 
theory of rectification, proposed by A. H. Wilson (1932), in which he assumed the 
existence of a surface barrier penetrated by tunnel effect, gave the wrong sign to the 
rectification. The more recent theories of Schottky (1939), Mott (1939) and Bethe 
(1942), which give the right sign, assume that no tunnel effect occurs. These 
theories give rectification in the observed direction. It will be shown here that a 
combination of the two theories—a Schottky barrier with an insulating barrier in 
addition, penetrated by tunnel effect—still gives rectification in the right direction. 

Wilson (1937) took a barrier of the form shown in figure 8. He assumed 
exchange of electrons to take place entirely by tunnel effect, and obtained an 
expression for the current 


j= NeNa/(RT/2nm) (eo 1a (21) 
B=2ng'. 


V is positive when the metal is at the higher potential. According to equation (21) 


the current is greater then, than when an equal, negative voltage is applied. ‘This 
is contrary to experience. 


Figure 8. | Wilson’s potential Figure 9. Schottky’s space- Figure 10. Surface film and 
barrier between a metal and charge barrier. space-charge barrier. 
a semiconductor. 


Schottky (1939) assumed that the only potential barrier at the contact is that 
set up within the semiconductor by space-charge, when electrons are lost to the 
metal by the impurity centres nearest the surface. This type of barrier is shown in 
figure 9. e¢ is the height of the conduction band of the semiconductor above the 
Fermi level in the metal, and is given by the difference e(¢, —¢,) of the work 
functions of the metal and the semiconductor respectively. In general e6>E,, 
so that for simplicity it may be considered that all the impurity centres at 


distances up to x» from the surface are ionized. Then by Poisson’s equation, 
a?V |dx* =4rp/k, 


V = (2a Nehey(2nny x2 ee eee (22) 
where the total potential drop 
Vyo=¢4+V =27Nex,?2/k. PAD as) 


In order to calculate the current, either “diode theory” (Bethe 1942) -o8 
diffusion theory can be used. In the former it is assumed that the barrier is thin 
compared to the mean free path of the electrons. Then onthe average an electron 
will not suffer collision in passing over the barrier, and the current will be given by 


i= (o/N)4/(2rm|RT) (1 — e-eV RT) eed RE, 
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When the barrier thickness is very much greater than the mean free path the 
diffusion theory must be applied. The current is then obtained from the equation 


i=e{ eee - FI, Ai (25) 


RT idx -dx 


where the mobility v=eD/kT, and D is the diffusion coefficient. Using the 
. . y 7 ib . 8 
integrating factor e®’’*”, this gives for the current, 
i| Sei e610 rl) ae (26) 

¥ Q 
where A is a constant. 


Wihenw—0./(a,)—o and m(05)=.V,,so.that 
A=Ne-ev/RZ, 
When 740, V(x») =¢+ V, therefore 
if * peVikT dx =eNDe-*#!*"(1 — ¢- eV IRD) 
0 
= loecere (len ercd 2) cain leita alt (27) 
The integral on the left hand side depends chiefly on the field at the surface and 
varies much more slowly with V and 7'than the terms on the right. In fact, for 
V <¢, when the barrier is altered little by the applied voltage, equations (24) and 
(27) are indistinguishable. With the usual assumption that e6>E,, V has the 
form given by equation (23), 1.e. 
iP DevIBU) Tye . e2INeX(a?—2ate)/KRL yy 
i) ¥ 0 
The vaiue of the integrand is a maximum at x=0, and decreases rapidly as x 
increases. When x 92>1/47Ne?, the integral has the approximate form 
ie eeVIRT dx = he ee efx y/KRT dx. 
0 ~ 0 
Since practically the whole potential drop takes place in the blocking layer, 
values of xx) are equivalent to x=oo. The integral can then be evaluated 


and gives 
| oe AAW ebato/KRT dx = KR T/417.Ne?x>. CORE (28) 
0 


The current is then 
t=ar/{8rNex(P+V)}(L—e Rte eee, (29) 
Both equations (24) and (27) show rectification in the direction commonly 


observed. 

The contact between a metal and a crystal of the type considered in § 1 must take 
place through a surface film. Both types of barrier will then be present, and the 
contact may be represented by figure 10. 

The total potential drop across the barrier falls into two parts, V, from x= —6 
tox=0, and V, from x=0 to x=. Since there is no change on the surface, the 
field at x =0 must be the constant field in the insulator and therefore 


Vi=FySy=47Nexy/x, 


V,=| Fdx=2nNex,2/«, i 
19) 
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where 


V,+V,=¢+V=9¢,—¢2+V. 


¢, and ¢, may be assumed to be so large that transfer of electrons takes place 
entirely by tunnel effect. For small voltages, i.e. | V| <d, the current calculated 
on the diode theory is 


i= (o/N)a/(2 [RT )e-2*82"8 [ee VadlRE _ geVRT] (31) 


‘On diffusion theory it is 
i= {ss"[(s-+91)} [e-tP~ PORT _ pV BH 030) 


where 


s=av/{47Ne*x($ + V)}, 
s’=0/N4/(2z). 


Equations (31) or (32) reduce to 
(24) or (27) when V,=0, i.e. 6 =0, and 
to (21) when V,=0, i.e. 6, ~,, and, 
therefore, x»=0. In the intermediate 
case, rectification may be in either 
direction, but will not be so pro- 
nounced, and the current in the direc- 
tion of high resistance will not saturate 
at K=kT. ‘This is" illustrated spy 
figure 11. 

The equations (21), (31) and (32) 
are derived on the assumption that 
an appreciable current is obtained 
under a small potential difference. If 
either ¢ or 6 is large this is not so. 
Then a voltage must be applied to the Figure 11. Alternative types of rectifying curves. 
contact, which is comparable to the (@) z=conse eeviRT a), 
contact potential difference before an (6) i=const. {e— e(O— Vi/RT = e— eV [RT 
observable current is obtained. The (8/A=1; x9/A=5). 
current fromametal, throughasurface ©) ?=const. (1- e~eV/k7), 
barrier in strong fields, is given by the 
Nordheim-Fowler formula (1928) 


6 WP (VY? ses ena 
t= G+nee (5 ) € ae ae a (33) 


From the semiconductor it is 
t=(c/N)\/(2nm/kT) grate CATS ~Panoie (34) 


When ¢, <¢,, the latter current will be the greater. Wilson’s theory then gives 
the right direction of rectification. 

Since the space-charge is eliminated when V <—4, the resistance of the 
contact vanishes in Bethe’s theory and there remains only the resistance of the 
a For V>RT/e, equation (31) shows that the current reaches the 
value 


i=(c/N)\/(2nmkT), == st. (35) 
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which is independent of the applied voltage. In strong fields, however, the barrier 
is distorted by the image force, so that for V~d 


t= (a/N)\/(2nm/k ih eat 


A voltage such that V, = —¢, applied to the surface film contact, removes the 
space-charge barrier, so that the current is limited only by the barrier of the 
surface film. For —¢,<V<-—4, 


i=(o/N)/(2nm|RT e292" V — 9/462), 


and for V< —¢4, 
EIGEN (Canukeije 8 ro) 2 9 Lee (38) 


when V,>¢, electrons are drawn from below the Fermi level of the metal. 
The number of electrons crossing against the field is negligible, and the current is 
therefore given by 
a 
=e | N(w)P(z) dz, (39) 


u— (p + Vx) ne 
where 


P(w) we-'2k(c— wy" —9/4(e—w))}, 


-N(w) ~(4rmejh?)(— 2). 
Integrating, this gives 


: 47 me?2 e “1 — Vi/4Px) 5 y _V/ad NY) ry 
L_= a . (1+ V,/4d,) (2¢,)? [1 — Ca (1+ V i/4 DV a: P2044 +. x(1 + V,/4¢,)}], 


where 
a = 2Kh, 10. 


When V,>4,, equation (33) may be applied, with V replaced by V,, ie. 


ane Cas Pag Vi\? ,4x6,270/37, 
Sopa Spal 5) é tee roe (41) 


where, from equation (30), 
. V/V =2%6/2?. 


This set of expressions, together with (31), covers the whole range of voltage 
for the surface film contact. Figures 12 (a), (6) and (c) have been drawn for suitable 
values of the constants. The number of unknown parameters, N, ¢,, ¢, and 4, is 
too great to allow a quantitative comparison with experiment, but certain conclu- 
sions can be drawn from the shape of the curves obtained. 

Figures 13, 14 and 15 are examples of experimental curves obtained for three 
different kinds of crystals. They are compared with the theoretical curves to which 
they most closely correspond. If the values assumed for N, ¢, and 4, were of the 
right order of magnitude, it would appear that any film (if it exists) on germanium 
crystals is less than 5 a. thick and has no influence on the rectification. On the 
silicon crystals it is between 5 and 10. and on silicon carbide between 10 and 20a. 
thick. The effect of increasing film thickness, for| V| <¢, is to decrease the degree 
of rectification and to prevent the saturation of current in the direction of high 
resistance. 


(0'Volt,/contact 
i 


10° 


Figure 12. The current-voltage 
characteristic of the contact be- 
tween a metal and a semicon- 
ductor, for three thicknesses of 
the surface film. 
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Figure 14. Curve showing rectification in the region 
V at the contact between silicon and a metal. 
Expeximental data supplied by British Thomson- 
Houston Co. Ltd. 


(b) 6/A =1-0 
(c) 6/A =2-0 


Curve showing rectification in the region V-~¢ at the | 
contact between silicon carbide and a metal. 
data reproduced from a report by E. W. J. Mitchell and 
Sillars, by 
Electrical Co. Ltd. 


Experimental 


permission of Metropolitan - Vickers } 


(a) 6/A=0-5 


Figure 15. Curve showing rectification in the regic# 
V at the contact between germanium and |f 
metal. Experimental data supplied by Briti 
Thomson-Houston Co. Ltd. 
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If the temperature 7 is changed, the thickness x, of the space-charge barrier, 
and therefore the critical ratio 6: x, changes proportionally. | Where the therm- 
ionic factor e~**/*? also enters, the total temperature effect is complex and is 
not considered here. For the germanium crystals, where § may be taken as zero, 
one has more simply 


i ocoe*P/FF — const . e~Pel?k? g-ed/kT. 


as in Bethe’s theory. 
§4. CONCLUSION 
The theory, applied in §1 to the contact between identical semiconductors, 
accounts qualitatively for the behaviour of the contacts between a semiconductor 
and a metal, if the effect of space-charge is taken into account. It follows that both 
the Wilson and Schottky mechanisms are necessary to a complete treatment o/ 
crystal rectifiers. 
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ABSTRACT. The specific heat at moderately high temperatures has been calculated 
for a one-dimensional lattice in which both electrostatic and repulsive forces are taken 
into account. The Thirring expansion method, which does not necessitate a knowledge 
of the density of normal vibrations of the lattice, has been applied and the resulting values 
of the specific heat are compared with those obtained by graphical integration of the 
frequency spectrum. outs 

Thirring’s method involves expanding the specific heat in a series in which the terms 
include summations over the 2th powers of the normal frequencies. It is found that 
values of n greater than 4 lead to computations which are excessively laborious. Nevertheless, 
with four terms, it is shown that reliable results may be obtained over a considerable range 
of temperature. At lower temperatures the Thirring method becomes progressively 
less accurate, showing that the number of terms used is insufficient, and a correction is 
derived from considerations of the remainder of the series after four terms. 
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von Karman (1912) who worked out the distribution of normal vibrations 

on the assumption of quasi-elastic binding forces between neighbouring 
particles. he specific heat of such a lattice was evaluated by Blackman (1935), 
who compared his results with those derived from the continuum theory of 
Debye (1912). The case of the ionic lattice was examined by Broch (1936) who, 
using the generalized zeta-function, worked out the vibrational spectrum of a 
one-dimensional lattice taking into account the electrical interactions of the 
charged particles and the repulsive forces. 

In this paper the specific heat of a linear ionic lattice is calculated by two 
methods. If the vibrational spectrum is represented by a function p(v) giving the 
number of frequencies between v and v+dy, the specific heat at constant volume 
can be obtained in the form 


x properties of a one-dimensional lattice were first investigated by Born and 


vmax hv 
cy=h | POE (gr) dy, Neo (1) 
xe” 
where E(x) = (e@— 1p > 


Vmax, 18 the maximum frequency, A is Planck’s constant, R Boltzmann’s constant 
and 7 the absolute temperature. Once the function p(v) is known, the solution 
of (1) using graphical integration is relatively straightforward, and greater interest 
attaches to the application of a method due to Thirring (1913) by which the specific 
heat of moderately high temperatures may be obtained without a knowledge of 
the vibrational spectrum. 

In this method the mean energy of a linear oscillator 


hy/RT 
exp (hv/RT)—1? 


is expanded in powers of hv/RT, giving 


i= thv+kT 


a=RT{1— & (—1)"(B,|(2n)!)(hv/ RT}, 


where 8,, are the Bernoulli numbers; the specific heat, obtained by differentiation, 


is thus expressible as a series of terms in v”, v4 etc., which have to be summed over _ |} 


all the frequencies of the lattice. ‘These terms become progressively more 
laborious to compute, and it is desirable to determine whether a comparatively 
small number of terms will yield useful results. In particular, it is important to 
know the temperature range over which the method is reliable, and the degree of | 


accuracy to be expected. Some information may thus be obtained for the three- | 


dimensional case. 

We consider a chain of particles which may be taken to consist of a number of |} 
cells each containing two particles. The coordinate of the kth particle in the /th_ |} 
cell in the equilibrium configuration is given by | 


x, =, + la, 


where x, is the coordinate of the kth particle in the base cell and ais the lattice] 
constant. The distance between the k’th particle in the base cell and the Ath | 
particle in the /th cell is denoted by x4,,, and the total potential energy of two |] 
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particles with coordinates -x;, and x, is denoted by ¢4,. The equations of 
‘motion then assume the form (Born 1923) 


fome[! ee] ee. | 


2 ; ae - ae (2) 
4p o0 | he lee | sake a aa aa 
kk’ 

‘The symbols vx |e defined by the equations 

freee 7, eee 

x (Giy)oe CXP{—2rthxy lay, sane (3) 

x x l= —a 

where ees= (Beda gy eee (4) 


Further, h/a is the wave-number, w is the angular frequency, U, and U, are the 
amplitudes and m, and m, the masses of the two particles. 

The potential 4/,, may be conveniently resolved into two parts (Broch 1937), 
subsequently referred to as B.) : one due to the Coulomb forces and the other 
due to the repulsive forces. ‘The van der Waals attractive forces are neglected. 


/ 


The functions fe | are therefore similarly separable, and for the electrostatic 


parts one obtains (B., equations (10), me and (27)) : 


ee Ss See om aoe (5) 
ae, . a : ue 24e 13) 28 ys, 
x x aan Ve 
oe (6) 


where the accent on the summation sign signifies the omission of the term /=0. 
Here {(z) is the Riemann zeta-function and the charge on the jth particle is +e 
or —e according as 7 is odd or even. 

Assuming that the repulsive potential varies as the inverse sth power of the 
distance, and taking into account nearest neighbours only, the repulsive parts 


IR 
i a may be expressed in the form (B., equations (34), (35), (37) and (48)) 
nM 


E al =e al = = cos-mhi(s-+1)log2 merce (7) 
= he ala -|. Al = — *F (e+ I) log?. eens) 


Returning to the equations of motion (2), we see that the normal frequencies 
are obtained as roots of the determinantal equation 


wm, +A(h) —_B(h) 
=| Oy SE 6 (9) 
Bh) wm, + A(h) 
Be? 22 +e 2nihl 
where A(h) = = [3£(3) —(s-+ 1) log 2] + roe a Ce (10) 
Be? 22 ., t% exp (2nihl) 
and B(h) = 4 cos mh(s + 1) log2— Fe ora, SP iad e (11) 
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In order to consider the distribution of the frequencies we specialize to a 
chain consisting of N cells, ie. of 2N particles. Application of periodic 
boundary conditions then gives 


exp (27zNh) = 1. 


Now 0<h<1 (since the wavelength cannot be less than the lattice constant) ;. 
consequently we have 


h=piN. (le a eee (12) 
where p is any positive integer such that 0<p<N. The allowed values of h 
have a constant spacing (=1/N) and the normal frequencies are distributed 
uniformly in h-space, each value of h corresponding to two frequencies. 
In the method of Thirring (1913) the heat capacity of a linear oscillator 
(of frequency v) is expressed in the form 


ere oes | Bike ee 
cmki lt SIS (er) 


and the series is convergent if hv/kRT<27. In this expression B, are the Bernoulli 
numbers. 

In order to derive the heat capacity of the one-dimensional lattice, it is necessary 
to sum this expression over all the 2N frequencies of the lattice, 1. e. 


ate Bit yea 
=k{2N+ © (-—1)"———— (=) 2 
cs R{2N+ E ( 1) (any! (en) ee } ie sees (13) 


We thus have to perform the summations over all the normal frequencies using 
equations (9), (10) and(11). For the case in which the two types of particle have 
equal masses, we write 7, =m =m and the solutions of (9) become 


An?mv? = — A(h) — B(h), ) 
2= — A(h)+B(h). J 


From these solutions the sums = v?” may be evaluated, and the cases corresponding 


47 Mov 


to n=1, 2, 3, 4 are treated individually in the next four sections. For n>4 the 
labour involved is prohibitive and a discussion of the error incurred in neglecting 
higher order terms is given later. ‘The equations (14) show how two frequencies. 
are associated with every value of 4; the frequency curve splits into two branches. 
but these two solutions are identical apart from a difference in phase (B., p. 493). 


S2- EVALUATION ORS sr? 
From (14) we obtain 
vit v}= —(1/272m,)A(h), 
and, summing over both branches to include all frequencies, we can express the 
required sum as 


ty? = —(1 [2n*mg)® A(h), | ees (15) 


where the right-hand side is to summed over all the permissible values of h. Now 
from (12) it is clear that the summation over h is equivalent to a summation over fe 
and, since NV is very large, this may be replaced by the integral 


N 1 
| AWN) dy, ive. by | NA(h) dh. 
0 
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Reference to (10) shows that A(h) consists of two parts, the first independent of h 
and the second an infinite series involving h. The latter is uniformly convergent 
in the interval 0<A<1 and may therefore be integrated term by term. This, 
however, gives zero for all integral values of / and we are thus left with 


Dv? = —(4Ne2/72mga?) [30(3)—(s-+1)log2]. sea (16) 


To obtain the numerical value of the expression in square brackets (subsequently 
denoted by C) the exponent s of the repulsive forces is given the value 10. Further, 


(3)= = 1/n8=1-20205...., 
1 


“= 


so that Cra 4018452 = Ee ee. (17) 
Substituting in (16) gives 
Pee OD IEN eV igds oN we Te <nssee (18) 
$35 BVAL UAT TON “OF =r 
From (14) we obtain 


me vail rl 
Det = N/8a4m? | | A2(h) dh + | Bh) ah ee PY eden (19) 
v 0 0 
Now A?(h) = 64e4a~ 6C? + 32e4a~°Co, + 4e4a- a2, see (20) 
where Cj a a : 
i} 


Of these three terms the first is independent of h, the second provides a zero 
contribution and the third may be expressed as a double sum by means of Cauchy’s 
theorem concerning the multiplication of absolutely convergent series. Upon 
integration the non-zero terms reduce to 


+0 Se 
Y’ 1/n8=2 > 1/n’, 
il 


and, since Pie OL I34 x 
n=1 
1 
we have | A*(h) dh=1:042 x 10% <*/a8. ww sw ee (21) 
0 


Squaring expression (11) gives 
Bh) = 64e4a-®[(s + 1) log 2]? cos? wh — 32e4a~*(s + 1) log 2 cos wh.o + 44a" *.0”, 
NOROEE (22) 
Secs oe OP [ath(21+ 1)] ; 
<a, [f+3\? 


Here the first term is straightforward, and the second splits into two sums both of 
which give, when integrated, a numerical value of 8. Application of Cauchy’s 
theorem to the third term in (22) leads to the expression 
+n 1 
x ———=2 2 — jy 
na—ol@+3|° ; n=o(n+3)° 


= 126 £(6), 


using (71) with x=6. 
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1 
Hence { BXh) dh =421-2 4/25, eee (23) 
0 


and, inserting (21) and (23) into (19), we have 
Lv=1871 Nee) ee (24) 


§4. EVALUATION OF 2% 
From (14) we obtain 
1 1 
Lv8 = —(N/32a®ms) | | A3(h) dh +3 [ A(h)B*(h) ah | noes (25) 
v “0 “0 
Now 
A8(h) = 5128a-C3 + 3848a C20, + 96%a- °C oo? + 88a 04°... (26) 


Here the first term is constant, the second vanishes when integrated, and the value 
of the third may be obtained from §3. The fourth term reduces to the double sum 


Yd’ [| m+n)3| m3 | n/3]-, (m+n<0), 


which is shown in Appendix 1 to have the numerical value 0-816. Integration 
of (26) thus gives 


“1 
| A%(h) dh= —3-400x10%8/a% sae (27) 
0 


For the second term in (25) we find 
A(h)B*(h) = 512%a-®C[(s + 1) log 2)? cos? rh + 32€8a°C a.” 
— 256¢8a°C(s + 1) log 2 cos zh.o, + 128€8a~*[(s + 1) log 2]? cos? zrh.o1 
+ 8¢8a~®o,057 — 64e8a-%(s+1)log2cosmh.oyog, == snare (28) 
Of these six terms the first is straightforward and the second and third may be 
evaluated using the results of §3. The fourth term splits into three sums: one 


gives a zero contribution and the other two have the same value, namely unity. | 


The fifth term may be expanded as a triple sum which reduces to 


22 [lanl A lan ee (m+n#A-—1), 


m nt 
and by adding the fifty largest terms a numerical value of 139 is obtained. 
The sixth term in (28) may be written as two double sums which lead to the 
expression 


4 | Sim 19] 2m 1-245 [l om 2m ete | : 
It is easily seen that each of these terms is equivalent to 


x [m3(2m+1)8}-2+ LD [m3(2m—1)3}, 
m=1 


m=1 
and the value of this combination (obtained from the general case in Appendix 2) is 
80 — 48¢(2), 
1.€: 1:0432 - ((2)=1-6449....,.). 


Integration of (28) thus gives 
1 
| A(h)B*(h) dh= —1:278x10%e8/a®, sid. (29) 
0 


and substituting (27) and (29) into (25) yields 
Lvt=2351N (e/a. eee (30) 
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§5. EVALUATION OF 218 
From (14) we obtain 


Ey = (N/128n%m§) lI. | An) ath+6 | A®(h)B%(h) dh + ls BA(h hy dh |. 


Now 
At(h) = 4096<8a-2C4 + 4096c8a-¥2C 8g, + 1536c8a-2C20,? 
200cas “Co, + leas to; 7 wees (32) 
Here the first term is constant, the second vanishes when integrated, the third is 


evaluated in § 3 and the fourth in § 4 (using Appendix 1). The fifth term reduces 
to the triple sum 


L’ DD’ [13] m|3| n|3| 2+ mtnP},  “(l+m+n<0), 
~mn 


and addition of the forty largest terms gives a numerical value of 6:79. Integration 
of (32) thus gives 


; | 
i A(h)dh=11-178x1088/a® aa (33) 
0 


The second term in (31) is obtained from 
A*(h)B*(h) = 4096e8a ?C[(s + 1) log 2]? cos? zh 
+ 256¢8a-l2C2o5 — 2048e8a-12C(s + 1) log 2 cos zh.oy 
+ 2568a—12[(s + 1) log 2]? cos? rh. + 16c8a-"a 73 
— 128¢8a-l*(s + 1) log 2 cos zh.o,0,? 
+ 2048«8a-12C[(s + 1) log 2]? cos? zh.0, + 128¢8a-!2Co,03 
—l024e"a-C(s4+ I )log2cosahiayg, ~~ = ne nye (34) 
Of these nine terms the first may be integrated straightforwardly and the second 


and third evaluated using the results of §3. The fourth term splits into three 
double sums which lead to the expression 


| 2 [| m2 |3| 2 + 1|3]-2 +X’ [| m9} m— 12423" [mis | : 


This reduces to z [m3(m + 1)3}-*+ x m* (cf. equation (63)) 


m=1 m=1 
and, using the results of Appendix 1, the numerical values are found to be 
0-1304 + 1:0173 = 1-1477. 
The fifth term in (34) may be expanded as a quadruple sum which reduces to 
ee al Ras [E+m+n+ PP", 


and by adding the ae cae terms a numerical value of 309 is obtained. ‘The 
sixth term in (34) splits into two triple sums and gives 


4 | = D’ [| m|3| 2 [3] 2m +204 1/3}-14+ 2’ d’ [|m|9| 2|3| 2m + 2n—-1 P| ; 
m n mn 

in which both terms have the same value, since changing the signs of m and m in 

either sum gives the other. By direct addition of the forty largest terms this value 


is found to be 2-345. 
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The seventh, eighth and ninth terms in (34) are similar to the fourth, fifth and 
sixth terms respectively in (28) and hence integration of (34) gives 


{ * 4%(h)B(h) dh=3-900 x 1058/a, eee (35) 
0 


Lastly we consider the third term in (31), which is obtained from 
B4(h) = 4096e8a-!2[(s + 1) log 2]4 cost wh — 4096e8a-[(s + 1) log 2]? cos* zh . oy 
+ 15368a-?[(s + 1) log 2]? cos? th . a? — 256e8a-(s + 1) log 2 cos wh . 09° 
+ 16e8a=12 G4, (36) 
The first term here is straightforward and the second, written as four separate 
sums, may be shown to have the value 
1/8 [8/27 + 24 +24 + 8/27] =6-074. 
The third term may be expressed as three double sums which reduce in the usual 
way to 


16 | Sri2e+10 21+3)3}2 +5 [| 2/+1/3| 21 1)8}4425/21+1/ | 
1 1 


In this expression the first two terms have the same value (the second emerges 
from the first on replacing / by —(/+1)). 


Now S [2+ 1)]§| 21+ 3|§}#=142 © (21-1)3(21+ 1), 
U t=1 
and SE (21—1)-3(21+41-3 = 5 (42-1) = (32 —3n2)/64 = 0-03736 
t= = 


(Knopp 1928, p. 269). 
Further, from (71) 


| 21+ 1|-8 = 63 £(6)/32, 
1 


so that the numerical value of (37) is 98-483. 
The fourth term in (36) splits into two triple sums which lead to 


256 | SE{2m+1)% 2n-+1|8| 2m + 2n+ 3|3}4 


mn 


+22 [| 2m+1|8| 2n+1)3| 2m+2n +1) | ; 
and by adding the forty largest terms the value of each sum is found to be 3-12. 
To evaluate the fifth term in (36) we have to consider the expression 


2 SUX [| W+1|9| 2+ 1/9] 2n +13] W+ 2m + 2n+ 3p, 
1 


mn 


and by direct addition, ignoring all terms in which J, m, are greater than unity, 
an approximate value of 2-86 x 104 is obtained. 
Integration of (36) thus gives 


1 
i Bh) dh=2-957x10%8/a®, =... (38) 
0 


and substituting (33), (35) and (38) into (31) finally gives 
2? = 3-091 N(e*/7,a°) (39) 
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§6. CALCULATION OF SPECIFIC HEAT BY THIRRING’S METHOD 


The specific heat, which may now be computed from (13) by inserting the 
results given by (18), (24), (30) and (39), is given by 


1-629 1-877 2-351 3-091 
CNRS 2 (= SE yen eaten We At 
{ ( 12 Ja+( nao) (eae) + (F300) \, 


ae (40) 


where A=(h/RT)? (€2/mya°). 

Since values of the specific heat obtained from this formula are to be compared 
with the results of graphical integration it is essential to ensure that the maximum 
frequency is the same in both cases. ‘The angular frequency may be expressed in 
terms of the wave number by means of the relation (B, p. 493) 


Myw? = 8e2a~*(s + 1) log 2(1 — cos zh) — 24e?a 33) 
2 oo) -e} a 2 
~4209) — = Tap) SS jee ara ay | cos 27hl 


2 
=1 
+ 00 


ta [exp(—afh+lP) + alk +P Bi (—a[h+ me) 


: 2 <e — {4+ 1]? G(,/nrf44+1 
+4ea3 {2 = | spore ( a ale Dy Se B 73 2] cos 27h(} +1) 


tn (—1)!' [exp (—a[h +1}2) +a[h + I? Ei (—afh +IP)]}. 


=—o 


In this expression G(x) =1—-— F(x), 
where F(x) is the Gauss error integral, and Ei( —~) is defined by 


0D pt 
—Ei(—a)=| Sas (0 <x <0). 


The w-h curve (B., p. 498) shows that there is a maximum value for the fre- 
‘quency when A is approximately 0-88*; the density of normal vibrations (given 
by dh/dw) therefore becomes infinite at this point. ‘The behaviour of the density 
function in the neighbourhood of the maximum frequency, though not trouble- 

some from the physical point of view, 

_ provides difficulties in graphical integra- 

tion, since it is not possible in this region 

to estimate accurately the area under the 
density curve. 

To avoid this difficulty we approxi- 
mate to the exact frequency spectrum by 
replacing the high-frequency end by a 
single “line” (i.e. by an Einstein term) 
whose frequency is equal to that corres- 

- ponding toh =1. Figure 1 shows this 
spectrum together with the exact relation 


according to equation (41). The range ok rT 
of h over which the Einstein term extends Figure 1. The angular frequency as a 
is decided by reference to the density function of h. The scales are identical 


curve, and the upper limit for the with those of B, figure 1. 


graphical integration has been taken to correspond to h=0-675. 


* The exact position of the maximum depends on the exponent s of the repulsive forces. 
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The frequency v,, of the Einstein term and the maximum frequency Vo of the 
continuous spectrum may be determined from (41). The series in this expression 
are rapidly convergent and numerical calculation is a fairly simple matter. 
Further, the functions F(«) and Ei(—.) are both available in numerical tables. * 
The calculations give 

Voy = 1:1769(e?/(tga®))?, wen wes (42) 
Vp = 1:1208(e?/(atga®)) we wes (43) 


In the Thirring calculation the maximum frequency is adjusted to have the 
value given by (43). At this stage the parameter , is introduced, defined by 


Oy)=hvglR)- 9. 5 eee (44) 


and the variation of the specific heat with temperature is obtained by plotting c, 
as a function of 6)/T. From (43) and (44) we have | 


h\?/ 2 -( 6 \2 
RT) \ma?] \1-121T) ’ 


and this expression, substi- 
tuted in (40), enables the 
specific heat to be calculated 
for any value of 0)/T. 

Table 1 shows the values 
obtained by this method over 
a range of temperature extend- 
ing down to 6)/T=4-0. The - 
table also shows the results 
of the graphical method dis- 
cussedin §8. ‘The results are 
shown graphically in figures 2 
and 3, in which the specific 
heat is plotted against 0,/T 
and 7/0) respectively. At 
high temperatures the specific 
heat tends to the classical ai 
value 2Nk. For values of 
9)/T up to 3-0 excellent agree- 
ment is obtained between the 
two methods, the maximum 
discrepancy (approximately °° 
0-5°%) occurring in the neigh- 
bourhood of 6)/T=0-7. For 
§,/T>3-0 the Thirring curve °* 
diverges from the graphical 
curve and gives higher values 


for the specific heat. This Figure 3. The specific[heat as a function of T/@. The | 
indicates that the formula (40) corrections given in § 7 are shown for 8)/T=3-0, 3°5 and 4-0. 


* 


° ‘ 2 3 = s 


e. 
& 
Figure 2. The specific heat as a function for 0)/T. The 
corrections given in § 7 are shown for 6)/T=3-0, 3:5 and 4-0. 


i 


— GRAPHICAL MeTHoD 


——— THIRRING MeTHoD | 


° o2 o4 os os ro 


a 
6. 


e.g. Publications of the Federal Works Agency, Works Projects Administration (New York 
1940-41). 
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is valid up to 6)/T'=3-0; above this value the results become rapidly less accurate, 
as might be expected, since the infinite series in (13) has been cut off at the fourth 
term. 


Table 1 
6o/T 0) 0:3 0-7 1:0 LoPs) Aled 2:0 2°5 3:0 3°5 4-0 5:0 


Cy/2Nk 1-000 0-995 0-974 0-948 0-921 0-890 0-818 0:741 0-667 0-611 0-604 
(Thirring) 

Cy/2Nk 1:000 0-993 0-969 0-945 0-918 0-892 0-817 0:740 0-667 0:587 0-519 0-407 
(graphical) 


§7. CONSIDERATION OF REMAINDER IN THIRRING’S METHOD 


Having found the range of 6)/T over which equation (40) can yield useful 
results, it is of interest to investigate the remainder after four terms in the original 
expression (13). In this section a rough estimate is made of the upper limit to 
the remainder in terms of 6,/T. 

An extension of Euler’s summation formula is used, starting with the identity 
(Knopp 1928, p. 528) 


Me es (a yes tel (ip ae df 46 

e*—1 TD aes ) (2n)1° 4 (Be i | 2 41%)? eee ( ) 

where Poy. 3(*)=(—19-? SE (2sin Qmx) (2m), oo. (47) 
m=1 


Now the mean energy of a linear oscillator is given by 
a =thv + hy/[exp(hv/RT)—1], 


and hence the indentity (46), with «=hy/RT, shows that the remainder after A. 
terms is 


1 
kT {arte hr | Pass(ae* ax 


Differentiating this expression we obtain the remainder R,= in the expansion for 
the specific heat : 


1 
R,= ka®*2(ex— 1) [ Poy y(x)e* dx 
J0 
1 
—(hhv]T) [2A+2)(— 1o**— ea2**| (A —1)-* | Py, galaer a 


—(hv/T) o2*+? (ex 1) { Poi (wje@de, =n anaes (48) 


To determine the maximum value of | R,|, with A=4, we consider the first two: 
terms together and write them as 


1 
Re {o10 (e% — 1)-1 — (1029 (e* — 1) — atte) (e*—1)-? } iF PA ee Ait ee are (49) 
An upper limit to the value of the integral may be derived from the inequality 


| Px) <2/2n)) & m™, 
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which follows from (47). Since 
Sm < Sms, and  m-%=n8/9450, 
m=1 m=1 m=1 


D 


A 


| P,(x)| is certainly less than the quantity P, where 


P=(056 4 94507)-2 (50) 
rl 


Thus for the limiting value of | P,(x)e*" dx we may write 
0 


lel ‘ 
P | ee dey = Px-Ye%—1), 
0 
which, substituted in (49), gives 


keP{—909 +029 + 09/(eX—1)}. eee (51) 
Since the last term contains e“—1 in the denominator, we may write e*~—1=« to 
obtain the form 
kPx(x—8), 
which indicates that (51) is negative, since the convergence of formula (13) requires 
a<6. 
We now examine the third term in (48), viz. 
fi rl 
fej a(e —1)+* | HE ale as 52) hs senate (52) 
0 
and seek to establish the sign of the in- 
tegral since this term has to be subtracted 
from the first two. ‘The function P,(x) 
has the form sketched in figure 4 (Knopp 
1928 p07): 
In view of the fact that the integral 
contains the factor e*”, where « is positive, Figure 4. “Form of the function P,(x). 
it is clear that the integral itself, and 
hence the expression (52), must be positive. 
To obtain the maximum value of| R),_4, therefore, the maximum value of (52) 
must be subtracted from (51). The limiting value of the integral in (52) may be 
expressed as 


P [se dx= Pa-2e%(a —1)+1}, 
which, substituted in (52), gives 
kP{oMee (e*— 1)-4 — 9}, 
Subtracting this expression from (51) we obtain —8k Px? so that we may write 
|R| <8kP(hy/kT)°, 


which, summed over all possible frequencies, becomes 


| R| <8kP(h/RT) Dv. EARN CE) 
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Now v° <vpXv8 and hence, using (39) and (43), we can define a quantity 
|R| max, by the equation 


| | mas, = BN(1-12 x 3-09) P(/RT)%(€2/mga®)9?. 
Finally we substitute from (45) and (50) and obtain 
(| R | max,/2NR) = (1-82 x LOs? ii) d= tz ee lal (54) 


Since R has been shown to be negative, this equation gives an upper limit to the 
correction which must be subtracted from the values given in table 1. 

Values of the remainder calculated from (54) are given in table 2 and the correc- 
tions are indicated by vertical lines in figures 2 and 3. 

The correction increases rapidly with decreasing temperature, since it is 
proportional to the ninth power of 6)/7, and although negligible for 6,/T<2:0 it 
is nearly 2°, at 6)/T=3-0. This is not borne out by the results obtained (table 1) 
and it is clear from the somewhat drastic approximations that the upper limit to| R| 
given by (54) is considerably in excess of the true correction. For the two cases 
6,/T =3-5 and 4-0 the values given in table 2 are approximately twice the discre- 
pancies between the ‘‘ Thirring”’ and “‘ Graphical” curves in figures 2 and 3. 


Table 2 
CAI 1-0 2:0 3-0 3°5 4-0 
[Rhmax./2NR 5x10-7 3x10 ~— 0-01, 0-05, 0-17, 


§8. CALCULATION OF SPECIFIC HEAT BY GRAPHICAL INTEGRATION 


The fundamental expression (1) may be evaluated by graphical means if the 
density of normal vibrations is known. In the present case, with the frequency 
spectrum illustrated in figure 1, the specific heat may be expressed in the form 


He hy hy, 
cy=k I, RON (Er) dv +2NRtgE (F) oe (55) 


where the function p(v) is normalized by means of the relation 


| Wea ones) ee (56) 


Here t, and ¢, are weight factors which represent the fractions of the total number 
of vibrations included in the Einstein and continuous regions respectively. Since 
equal ranges of A contain equal numbers of vibrations, we have 


EEO Gite r0:32550 1 eee (57) 


The frequency spectrum for the region 0<v<vy has been calculated from 
equation (41) and the values of w and h so obtained have been used to calculate the 
density function p(v) in terms of w by means of the relation (B, p. 499) 


2 ce Bs 2 a 
ees { (¢+1)log2sin mh+2 = [eter 2: sr | sin 2rhl 
l=1 


pian a Ni, 


<n E [1=(~1yN(e+) Bi(—afh+ 1) 


> | exp(—7ls+4)) | Giv7l+4))].. 1 
—2 2 joe a | sin 2rh($ a p}. Pareles (58) 
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For any value of 6)/T we find can E (hv/RT) as a function of » ; computation is 

facilitated by means of the relation 
E(«)=x7e7/(e? 1)? = 327) (cosh'x 1) ar (59) 
Hence the product p(v)E(hv/RT) can be obtained as a function of w. In this way 
curves have been constructed for values of 6/7 varying from 0-3 to 5-0, and the 
areas under the curves determined by counting squares. The second term in (55) 
may be calculated immediately from (59), since 
Gee =a lO Sion iin 

and the sum of the two contributions gives the specific heat, as listed in table 1. 

Referring to figures 2 and 3, the slight discrepancy in the neighbourhood of 
6,/1 =0-7 is rather puzzling, but could probably be decreased by raising the 
frequency v,, of the Einstein term or increasing its weight fy. At the low temper- 
ature end of the curves one would expect the graphical method to give reasonably 
accurate results, since in this region the specific heat is controlled by the low- 
frequency end of the spectrum and hence the effect of the approximation (figure 1) 
is least noticeable. 
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APD END Pixel 


In this appendix a method is given by which an approximate numerical value 
may be obtained for the expression 
8, =2) 2) [lense ale ee eee (60) 


mn 
(mA0, nA~0, mA~—n). 
Considering first positive values of m only, and putting n=«, we obtain series 
of the form 
8” [| 272 |3| me +o |8}-3, 
™m 


and in the two regions m>0 and m<—«, we have 
ao 
a? 2 (at 0; (61) 
m=1 
while between the two prohibited values of m (0 and —«) there are («—1) terms, 
_ whose sum may be represented by 


2 EPG Sot a, he eee Btn Geer ene (62) 
taken over all non-zero values of € and 7 satisfying the condition E+n=a. 
The same reasoning will apply to negative values of m._ For if we write n= —« 


and consider the series 
SBI bes 
8% '[} m9 mo), 
m 


it is clear that the two regions m>« and m <0 will each yield (61) and the remaining 
(«—1) terms will be given, as before, by (62). 


The specific heat of a linear ionic lattice 339 
We thus have the result 


ic @) 
2X" [| m|3| m+ «|2}44 2D’ [| m|3| m—a|8}1=4 ¥ m3(m + a)-3 42 DE 37-3, 
mM 


uk m=1 


At ecee (63) 
which we can apply to (60) and obtain 
Sy =4 X [mr 3(m+1)-7] a= | 2+4 x mn-4 2) | 
m=1 m=1 
+ ee 4 y [m-3(m + 3)-3] 
27.2 m=1 
+ gat gp tte rtm 64 
Gallons sp be” (ae | Ace ts: aaitieee (64) 
By expressing 1/[m°(m+«)*]-in partial fractions, we can derive the result 
ee ee in eel 
mai MP(m+a)j® oF Lm mt+a af Lm? (m+«a)? 


i 2 pal 1 
+35|a- @ae 6 Mo DOK (65) 


in which the three terms on the right may all be calculated easily for small 
values of «. The original expression (60) converges sufficiently rapidly for the 
sequence (64) to be cut off at, say, «=4. 

The sum in the first term in (65) is simply 


1+4+4+....41/c, 
while that in the second term is equal to the sum of the first « terms subtracted 
oO 
from 2%1/m?. In the third sum the terms cancel out with the exception of the 


1 
first « terms which donot overlap. Values of the expression in (65) for «=1, 2, 3, 4 
are given in table 3. 


Table 3 
a 1 y 3 + 


ao 
2 m(m+a)- > 0:13040 0:03940 0:01687 0:00864 


m=1 
Substituting these values in (64) we obtain 
SUS Oe ato. (66) 


Referring to (26) it will be seen that 8S, is small compared to 512C® and 
therefore, in this instance, a more refined calculation of S, is unnecessary. 


APPENDIX 2 
In this appendix a closed expression is developed for the sum of the series 
SK ( 3) ei 02 9) (67) 
m=1 e 


(GieciaNGREEN ern) 
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Using the method of partial fractions we find 


. le ee | 52 | a+ oa | 
mai ne(2m+ B)> — BiLm m+p/2 B47 Lm (ant B/2)° 


| 

telo-osaop | 

B24 Lm (m+ 8/2) | 

The three terms on the right will be considered in turn and it will be convenient 


to deal with positive and negative values of 8 separately. For the first term, 
with the notation 


<i) Pall 1 
== © Ss ae 
1S pel, E m+ | 5 | ; 


we have 
Mr] 1 Mr 1 1 
TET y pS ena 
Oy) cae 1 E m+ x | 1 E m+ |} 
Teles M+KG-1) 4 
= Lin) i} 
mM>o (1 mre? B41) Nts 

where n=m+4(P—1). 


Thussfor bei; 
1 WS” 1 2 1 1 1 
s(8) —s(1) = : ae (Grete +3): 


| Coed Care| 

Now 1)=2(5-3 re ey) =2(1-log2), 

and hence s(8)=2 (; + : - += —log2) wo 3° @ Oeeaee (69) 
For B<-—1, we have ; 

: 1 Leal eet 1 
—s1)=- & =—2(--=-—=—- oo +> 
s(P) < ) s(p+1)N+4 € 1 3 % on 
and hence s(B) =2(1 —log 2) —2 tees men (70) 
[9 3°. 6a oe 


The second term in (68) may be split up into two parts, of which the first is 
simply 


y 1/m? = (2). 
m=1 


For the second part we use the relation 


(1 =2-*)C(e) 1 3 oi (71) 
which, for x=2, gives : 
e 3(B+1) 
2. (a+ 8/2) 2 = 302) — | ao (72) 
1 1 


This holds for 8>1, while for B<—1 the corresponding expression is 
4 
3¢(2)+ & (m+4)?. 2 sssa5 (70) 
a(B+1) 
The third term in (68) is similar and we have, for the first part, 


ZI /m? = (3), 
m=1 
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while for the second part, using (71) with x =3, we find 


3(B+1 
0 7e)= bs ‘(m—4y°8, (RES ley Where (74) 
D (m +B/2)-8 4 sf 
1 [ =763)+ = (mb ays, CEA SY Nae (75) 


From (68) the original sum for 8>1 is obtained from (69), (72) and (74): 
Soro aol 1 6 MB+1) | 
95 a ery (Me ~, a Pa eCNO Sar lac? aa = 7 eames Fees 
= BT + 3454+ + g—lo5?2)—e(8O)—"2 
1 (= 1) 1 603 
ae SS PR) ae | Mn ae Ore tence. 
50 2 ee HO) (76) 


(m—} 


ila ike 


~ 5,(#42)+ 9 cas) 


(B41) (m+ 3)? 


ies eo! 
+ (= ao +6E(3)) 2 te RR ee (77) 


For the special case 8 = 1, (76) reduces to 
S,( + 1) = 80 — 6[8 log 2 + 42(2) + (3)], 
and for B= —1, (77) gives 
S,( — 1) = 6[8 log 2 — 4¢(2) + (3)]. 
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The Interpretation and Application of Electron-Diffraction 
‘Kikuchi-Line ’ Patterns—Part I. The Determination of the 
Crystal Unit Cell, its Orientation and the Crystal Symmetry 
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ABSTRACT. ‘The interpretation of electron-diffraction Kikuchi-line patterns from single 
crystals is developed as a powerful independent and general means of determining the 
crystal lattice type, dimensions and orientation and the Laue symmetry of the atomic arrange- 
ment. This aim is realized by a construction of the reciprocal lattice by methods which are 
described and illustrated. The satisfactory accuracy of the method is due to the relatively 
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precise definitions for the location of the Bragg reflection positions, derived and stated here 
for the first time, which enable the method to be used even when the crystal lattice constants 
are large so that the bands are narrow and diffuse-edged, as with rhombic sulphur. The 
patterns shown also illustrate features of the line intensities, especially the reduction in 
intensity when the bands pass near the undeflected-beam spot. 


Si IN LRO DUG TLON 


INGLE-CRYSTAL diffraction patterns provide much more complete thougn 

usually also more complex information than is yielded by the ring patterns 

obtained from polycrystalline material, with respect to the structure of the 
material, i.e. lattice type, symmetry and atomic arrangement. It is therefore 
especially important to be able to interpret and utilize single-crystal patterns in 
such a way as to reveal all the information available in a form easy of application 
to all cases. The interpretation of the single-crystal Kikuchi-line pattern is 
developed below as a powerful independent technique by which the structure 
and orientation of single crystals can be determined. 

Electron-diffraction patterns of black and white lines from a single crystal of 
mica were first observed in 1928 by Kikuchi, who showed the relationship between 
the line positions and the crystal net planes, and the rough agreement between the 
relative line intensities and the relative x-ray reflecting powers of the orders of 
reflection from a series of net planes. Since then, although such patterns have 
been used as a means of defining the orientation of crystals of known structure 
(e.g. Thirsk and Whitmore 1940) and of checking the structures of graphite and 
cadmium iodide (Finch and Wilman 1936, 1937 a) they have otherwise only been 
made use of in a general way, by inferring from their presence that there are 
relatively large thicknesses (several hundred a. or more) of coherently diffracting 
crystals in the specimen, and from their diffuseness or displacement the degree of 
fine-scale surface roughness, mosaic structure and distortion of the crystal. 

In order to use these patterns to determine the lattice structure, orientation and 
Laue symmetry of a crystal, all that is required is to construct the reciprocal lattice, 
and this can be achieved in the way described below. Moreover, it can be done 
independently of whether or not a clear pattern of Laue spots is also present in the 
pattern, and independently of any other source of knowledge. The crystal 
lattice and its orientation can be determined from a single photograph if this 
records a wide enough range of diffracted rays, i.e. the usual 10 to 20° from the 
primary beam direction. ; 

The method is also of general applicability, for it has been demonstrated 
(Finch, Quarrell and Wilman 1935, Finch and Wilman 1936, 1937a, 1937 b) 
that clear Kikuchi-line patterns can be obtained by transmission through thin 
cleavage flakes or natural lamellar crystals in favourable cases, or by reflection from 
(a) smooth natural faces, (b) cleavage faces, (c) ‘ parting’ fractures, (d) conchoidal 
fractures, (e) polished crystal surfaces (plane or curved), or (f) polished and 
lightly etched surfaces. It is the only stationary-crystal electron-diffraction 
method applicable if the patterns contain few or no Bragg or Laue spots owing to 
the narrowness of the Laue zones. In such cases the crystal surface can usually 
be slightly roughened by careful abrasion or etching (see examples of FeS, and 
ZnS respectively: Finch and Wilman 1937b) so as to yield more extensive spot 
patterns, or the crystal can be rotated during exposure of the plate, but the accuracy 
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of location of the lattice rows from the spacings and curvatures of the spot rows, 
i.e. the Laue zones, is in any case often very much less than in the Kikuchi-line 
method, in which angles between net-planes are directly measurable. 


$2. THE BASIS OF THE INTERPRETATION OF 
KIKUCHI-LINE PATTERNS 


(1) The origin and general nature of the lines and bands 


For the present purpose it suffices to apply, as Kikuchi did, Rutherford and 
Andrade’s (1914) basic conception of the line pairs as originating by selective 
Bragg reflection by the crystal net-planes from the diffuse background scattering, 
which has a radial fall of intensity with increasing angle of deviation from the 
primary beam. A practically parallel pair of black and white lines in the pattern is 
thus the intersection made with the photographic plate by the two halves of a 
cone round the net-plane normal, with semi-apical angle 90°-0; @ is the Bragg 
angle, related to the net-plane spacing d and the electron wavelength A by Bragg’s 
law, 2d sin 6=nd, where 7 is an integer, the order number of the diffraction. Itis 


found (Kikuchi 1928, Shinohara 1932, Kirchner 1932, Emslie 1934, Finch, 


Quarrell and Wilman 1935, Finch and Wilman 1936, 1937 a, 1939) that the 
measured line-pair separations correspond well (to within about 1 to 2°) with 


this simplified interpretation, the wavelength of the electrons which contribute 


to the lines being practically identical with that of the primary beam. 

In cases so far examined (Kikuchi 1928, Kirchner 1932) the relative 
intensities of the different line-pair diffractions also show a general correspondence 
with the relative structure factors and electron or x-ray reflecting powers of the 
net-planes when allowance is made for the different inclinations of the planes to the 
beam. ‘The relative intensities of orders of reflection from a net-plane are, 
however, appreciably modified by secondary scattering, or more correctly by 
dynamic interaction of the scattered wavelets, so as to even them out partially. 


_ In certain definable cases orders of reflection may appear which are forbidden 


according to the theory of single-scattering, but this cannot affect extinctions 


_ characteristic of the lattice type (Brandenberger 1937). 


Kirchner (1932) has pointed out that the elementary diffraction theory for 


, single scattering also leads to a Kikuchi-line breadth proportional to the corre- 
sponding plane spacing d/n, and inversely proportional to the thickness of the 


diffracting part of the crystal in a direction normal to the net-plane. ‘The lines 
corresponding to d/n about 2 A. or more are in fact generally very broad and diffuse, 
while those of higher orders from the same plane or from planes with lower d/ are 


progressively much sharper. The dynamical theory of scattering has been found 
by Lamla (1938) to give an expression for the line breadth proportional to 1/sin 28, 


which is .\early the same as Kirchner’s result, since 29 is small for fast electrons. 


Though the theory for selective single scattering from the diffuse background 
goes so far to account for the position and form of the lines it fails to explain the 
existence of the prominent bands with symmetrical intensity distribution, which 
_ still appear, except near the primary beam, even when the net-plane is parallel to the 


primary beam. Hayasi (1934), however, showed that such bands are to be 
expected if from each scattering centre spherical waves spread out with amplitude 
proportional to the cosine of the angle of deviation from the primary beam, 


and these are scattered again once by all the other atoms in the crystal. ‘The more 
complete dynamical theory, applied by Shinohara (1932a), also yielded an 
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intensity distribution across the symmetrical bands which was similar to that 
observed (figure 1()). For present purposes, however, these details of the 
intensity distribution need not be considered further. 


(ii) Location of the Bragg reflection positions and measurement of their separation D) 


For present purposes it is first desirable to define more accurately than has. 
hitherto been done the positions of the Bragg reflection directions relative to the 
Kikuchi lines and band edges, which are very diffuse when they arise from planes. 
whose spacings are large. 

In the case of the symmetrical bands associated with planes that are parallel to 
the primary beam it appears to have been generally assumed that the Bragg 
positions are at the band edges; these, as the author now defines them, are the: 
points where the intensity is zero relative to the general background in the region 
considered, as it passes from a positive value (excess over the background) near the 
maximum on the side within the band, to a negative value which reaches a minimum. 
of intensity slightly farther out. For example, Thomson and Cochrane (1939, 
pp. 110-116) seem to indicate this to be their opinion, and in indexed patterns of 
Emslie (1934), Tillman (1935), Finch, Quarrell and Wilman (1935) and Finch 
and Wilman (1936, 1937 a), the lines denoting the consecutive orders of Bragg 
reflection from such planes are drawn in accordance with this view. No close 
quantitative test yet appears to have been made, and this is now supplied by 
measurements on patterns taken at not too small a camera-length (~22 cm.) and 
obtained with a finely-focused narrow electron beam (~55 kv.) reflected from a_ | 
perfect magnesium-oxide cleavage face (cube face (001)), so that the lines are as. 
sharp as possible. | 

Figure 1(a@) shows part of~such a pattern where the median of the strong |] 
vertical band system due to the (100) plane lies in the plane of incidence, and the 
‘undeflected-beam spot lies on this median 3°8cm. below the lower edge of the 
figure. ‘The intensity distribution across this band system in a region where 
obliquely crossing bands are not prominent, is of the form shown in figure 1 (6). | 
In this figure are also shown measurements made on the plate at the 400, 600 and 
800 band edges, corresponding to points where the intensity was estimated visually 
to be the same as that of the general background in the same neighbourhood. 

It will be seen that there is very good agreement in the spacings for the strong 
and relatively sharply-defined 400 and 600 band edges, and only slightly 
less so for the faint but sharp 800. The positions for the very diffuse 
200 band edges calculated from these measurements were found also to | 
agree excellently with the points where the intensity was the same as that of 
the background. ‘These regularly spaced positionsare thus evidently the succes- 
sive orders of the Bragg reflection positions. It is also clearly shown in figure 1 (a). 
that the vertically elongated but laterally very sharp Laue spots lie in every case | 
within 0°001 cm. of these positions when the net-plane passes exactly through the 
primary beam. (Actually, in figure 1 (a) the Kikuchi line pattern is not exactly | 
symmetrical about the plane of incidence, but is displaced sideways by 0°007 cm. 
relative to the Laue-spot rows, due to an azimuth deviation of 0° 1’ from the 
symmetrical one, but the quantitative correspondence is clear.) This fact | 
confirms that these are the Bragg reflection positions for electrons which have | 
practically the same speed and energy as those of the incident beam; in the case. |} 
of the very sharp higher-order band edges such as 800 and 10, 0,0 the electrons | 
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contributing to the Kikuchi line pattern are calculated to have lost not more than 
50 electron-volts energy. 

In the case of net-planes not parallel to the primary beam, black and white 
line-pairs arise in general, but here also there is a band of increased intensity 
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Figure 1 (6). 


between the lines, which is strong when the line-pair separation is small, corre- 
sponding to planes of high spacing. This can be seen prominently in figure 2 (a), 
which shows also the sector-like effect caused by the crossing and superposition 
of parts of two or more such bands. In these unsymmetrical bands and line-pairs 
it seems hitherto to have been usual to measure the distance between the centres of 
the black and white lines, assuming these to occur at the Bragg reflection positions. 
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To test this, measurements have now been made on several successive orders of 
line-pairs or band edges such as those in figure 2 (a), which is also part of a pattern 
from an MgO cleavage face at an azimuth near that of figure 1 (a) ; thus the pattern 
adjoins on the right that recorded in figure 1(a). An example of the measure- 
ments, for the orders of diffraction from the (101) plane, is given in figure 2 (d), 
at the points indicated relative to the intensity distribution shown qualitatively in 
the figure. It will be seen from the upper series of measurements that the points 
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where the intensity is the same as the background intensity again forma regularly” 
spaced series as in the symmetrical bands. On the other hand, the lower series 0 
measurements shows that while the distance apart of corresponding main maxima 
and minima are virtually identical with corresponding spacings of the upper series, 
nevertheless as a whole these points do not form a regularly spaced set. The 202 
maximum and minimum are’so broad that it is unequivocally shown that they are 
considerably displaced (0-028 cm.) from the positions (at 0°826 and 1°718) which 
would form a regular sequence with the 404 and 606 maxima and minima. 

It is thus clearly seen that (i) the Bragg positions must occur in this case also, 
as in symmetrical bands, practically where the intensity is the same as that of the 
background; (ii) the same spacing is obtained whether the measurements are 
made between the minima and maxima (centres of the absorption and intensified 
lines) or between the above-defined Bragg reflection positions (inner edge of the 
absorption line and outer edge of the intensified line), though the measurement 
is more definite and easily made in the latter case, especially when the band-width 
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Figure 3. 

D is small and the band edges diffuse ; (iii) the line of intersection of the net-plane 
with the plate is relatively accurately obtained from the above-defined Bragg 
reflection positions. For the present purpose these results are mainly important 
in showing how greater accuracy of measurement of D can be achieved, thus 
permitting satisfactory analysis even of patterns in which many diffuse-edged 
bands occur, from crystals whose lattice axial dimensions are large. Measurement 
between points where the intensity is the same as the background also often 
allows relatively accurate measurements to be made, even when some of the lines 


are much broadened due to a curvature, mosaic structure or distortion of the crystal. 
§3. CONSTRUCTION OF THE PROJECTION OF THE RECIPROCAL 
IDA AIMED. IONE MENS) IANS. IPOS IAP IOIN'S 
(i) Calculation of reciprocal-lattice point positions from the line positions 


In figure 3 let the crystal origin be at O, near the centre of the diffracting part | 
of the crystal, and let OO,, of length L, be normal to the plate and meet it at O,. 
Let PO,P’ and OP represent in section the plate or the ‘positive’ print, and a 


‘Kikuchi-line’ patterns—Part I BAF, 


plane through O parallel to a net-plane system (k/) in the crystal, inclined at ¢ to 
O0,. Then the radius vector r* from O to the reciprocal-lattice point [[AA/]* 
is ha* +kb*+Ic* and is along the normal OP’ to this plane, with length 
1/d,,,,. To define the position of this point [[Ak/]]* and the symmetrical point 
[[—h, —k,—I]]* in space, it is only necessary to specify the direction and the 
magnitude r* cos ¢ of the projection of r* on the plane through O normal to OO, 
(parallel to the plate), together with the distance r* sin ¢ of the point from this 
plane, which is considered positive if measured from the plane along OO, and 
negative if along O,O. 

These magnitudes can be obtained easily from the camera-length L and 
measurements of the perpendicular p from O, to the line-pair median, and of the 
separation D of the line-pair K,K,, measured along the perpendicular O,P. Thus 


r* =1/(d/n) =sin0/A~6/A= (D’ cos ¢)/AL = (D cos? 4)/AL. 
cosp=Li(L2+p%)!; sing=pi(L?+p%)!;— tand=p/L. 
Whence 
Ty* =1* cos¢= (Dos? $)/AL 


WHY SY LONG = NY GEN Sey Gio 

(iz) -{-a(¢) + ez) ~Te(2) + 

=(1/AL).D. f(p/L), 

z=r* sin @=(p/L) . r*cos¢. 
Figure 4(a) can be used to find f(p/L) for any values of p and L likely to be en- 
countered normally, since under the horizontal p/Z scale are marked scales of p 
for a range of Z values. The projection of the reciprocal lattice is drawn by 
plotting as radius vectors AL . r* cos¢, i.e. D . f(p/L), to a convenient scale, for 
example a five (or sometimes ten) times scale, the corresponding distances z . AL 
of the reciprocal-lattice points from the projection plane (p/L).D. f(p/L), being 
written by each point as the actual distance in centimetres relative to this same scale. 
An exact expression for D, i.e. PK, +PKag, is obtained from triangles OPK, 
and OPK, in figure 5: 
PK, = L’' sin @/cos(¢+ 8); PK,=L’ sin @/cos(¢—8), 
thus D=PK,+ PK,=L' sin 20 cos ¢/(cos? ¢ — sin? 6) 
™ L’ sin 26/cos¢ ™2L’ sin 6/cos¢ 
= 21 sin 6/cos? ¢. 

The last approximate expression for D is accurate to within 0-5°% even when ¢ is as 
much as 20° (p~7:7cm. when L=22cm.) and @ as much as 4° (D~3cm. with 
L=22cm.), thus the same accuracy, sufficient for present purposes, is obtained in 
r* cos¢, because from this D value we have without further approximation : 


r* cos =(n/d) cos d = (2 sin 6 cos ¢)/A=(D cos? ¢)/AL’ = (D cos? 6)/AL, 
which is the expression already given above. 

Similarly, to find how far the trace of the net-plane on the plate is displaced 
from the actual line-pair ‘median’, midway between K, and K, and normal to 
O,P, we have 

(PK, —PK,)/D=L’ sin@ . [{1/cos ($+ 4)} — {1/cos (¢ — 4) }]/D 
=2L’'sin?@ . sin¢/(cos? d — sin? @)D 
=tan@tand 
=(p L)tano=D),/2L( L(+ p*); 
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by using the above accurate expression for PK; + PK,=D. Since 6 is never likely 
toexceed 6° in normal practical cases, it follows that even if p/L is as high as 0-3 and 
§=6°, the plane trace is displaced from the ‘median’ by only (0-3 x 0-1)D/2 
i.e. 0:03(D/2), which can be neglected for present purposes. 

A further consideration sometimes required is the variation of line-pair 
separation D,, with distance y along the ‘median’ or net-plane trace, since the 
distance apart D, i.e. K,Kg, at the hyperbolae vertices (y=0) is sometimes not 
directly measurable because this region lies outside the recorded photograph. It 
can be shown that 

D;, = D*® + 4y?/(cos? ¢ cot? 6 — sin? 4), 
whence 
(D/D,)? = 1 —(4y?/D,?) sin? 8/(cos’ ¢ — sin? 6). 
As shown above, 
D=~2Lsin 6/cos? ¢, 
thus, neglecting sin? @ in comparison with cos? ¢ we have, when D,=D, 


(D/D,)? = 1 —(9?/L?) cos? 6 = 1 —(y?/L*){L*/(L? + p*)5.- 
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Figure 4 (a). Figure 4 (0). 


Figure 4 (4) shows the values of D/D,, corresponding to this expression, over the 
range of y/L and p/L concerned, from which can be interpolated in any usual case 
the factor D/D, by which D,, must be multiplied to estimate D. 


(ii) The relation between a ‘positive’ print and the initial ‘ negative’ recording 


It is convenient to measure D and p of the line-pairs on either a ‘ positive’ or a 
negative’ reproduction of the initial photographic ‘negative’ recording of the 
electron diffraction pattern. 
it a negative reproduction is made by first printing a positive on to a plate and 
taking a print from this on paper, it can be considered in its geometry as identical 
with the initial negative recording. 
A positive reproduction, however, is usually preferred owing to its greater 
clearness of detail, and it bears in its geometry the relation of a mirror image to the 


\ 
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initial negative. More exactly, if the plate and print are viewed face side up,. 
with top edges roughly in line, and then by rotation of one a line BA in the print is 
placed in line with the corresponding line AB in the plate, then the print is the 
mirror image of the plate (apart from interchange of black and white) in a plane 
which is normal to the line chosen, and cuts it midway between the pair of ‘A’ 
or‘ B’ points. 

Alternatively, if the positive print and the plate are put in register with the 
Sensitive plane faces together and are then opened apart by rotating the print about 
any line in their plane until the sensitive faces are again coplanar, then the print 
and plate are mirror images of each other relative to that line. It follows that, in 
respect of its geometry the ‘ positive’ print can be regarded as identical with the 
initial ‘negative’ which would have been obtained if a crystal had been used 
which was the mirror image of the crystal actually used, in respect of its lattice, 
orientation and atomic structure. If the crystal lattice and the atomic structure 
both possess a plane of symmetry (or have a Laue symmetry containing one) then 
the details of the positive print are congruent with those of an initial ‘negative’ 
recording obtainable from an actual crystal if it were in an orientation which is the 
mirror image of the actual orientation in a plane perpendicular to that of the 
negative and its positive print (both with sensitive sides up), and passing through 
thelineabout which the print was rotated to being it facing the same way up as the 
negative. 

It is this equivalence which has been invoked hitherto in most cases where 
‘positive prints of photographs from a crystal known to have a plane of symmetry 
have been indexed, both Kikuchi-line and spot patterns having been indexed as if 
they were the initial ‘negative’ recording (apart from interchange of black and 
white) viewed along the beam direction from behind the crystal (e.g. in Finch, 
Quarrell and Wilman 1935, and Finch and Wilman 1935, 1937 a, 1937b). 

| However, it is preferable in such cases, and it is necessary in the case of crystals 
-of unknown structure or those known to have no plane of symmetry, to adopt the 
following unambiguous system which will be applied in the examples given below. 
‘The pattern, whether a positive or a negative, can be thought of as situated so that 
‘the features of its intensity distribution coincide with the corresponding distri- 
bution of the electrons in the recording plane; a ‘negative’ print is then to be 
-considered as if it were the actual negative recording viewed along the primary 
beam direction, from behind the crystal! (figure 3), and a ‘ positive’ print is to be 
‘regarded (apart from interchange of black and white) as if it were the initial negative 
‘recording viewed from behind, in the direction opposite to that of the primary 
beam, i.e. looking towards the crystal. 


(iii) Practical details of pattern measurement and projection construction 


Location of O,. Though it is not absolutely essential to locate O,, the foot of 
the perpendicular from the crystal to the plate, to within more than about 2cm. 
when L is about 20 to 25 cm., it is obvious that the more accurately this is done the 
more consistent will be the estimated lattice-point positions and the accuracy of 
determination of the lattice symmetry and dimensions. One indirect way of 
locating O, sufficiently closely is from observations of the positions of the focused 
undeflected-beam spot F (figure 5), recorded in the photograph, and of the 
unfocused-beam spot U, the camera-length L, and the distances BA and BC of the 
-focusing-coil from the anode diaphragm A and the plate. ‘Thus, approximately, 
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by assuming that any change in beam direction caused by the focusing-coil occurs | 
near the central plane of the coil, we have 


CO,~BE+E,0, 
~» BE+(CF —BE). {1 —(O0,/BC)} 
= BE.OO,/BC+CF{1—(O0O,/BC)} 
= (OO,/BC)(AB/AC) .CU+ {1-(00;/BC)}. CF. 
AC, AB, BC, CU and CF can all be measured, and the vectorial form of the 
expression applies to both the cases (a) and (d) of figure 5, though that of (a) gives. 


the best beam focusing, making F as near U as possible, and it should therefore 
always be used. If F and U coincide, 


CO,=CF (AC— O00) AC. 


(a) 


Figure 5. 


Measurement of D and p. When a print of the negative has been prepared and |} 
stuck down on a suitable-sized piece of paper, the line pairs are sorted out and | 
extensions drawn round the print so as to label each of the line-pair medians |} 
with a serial number, as in figures 7(a) and 8(a). Corresponding black and |} 
white lines are practically parallel and of similar intensity and sharpness. | 
The fainter higher orders are not always visible continuously across the | 
whole photograph, but usually appear more clearly on regions of raised intensity | 
due to other band systems. ‘The line-pair separation D at or near the foot of the | 
perpendicular p from O, to the lines are then measured on the print, using a) 
thin-edged finely-graduated millimetre scale to avoid parallax, and measuring D 
usually to an accuracy of 0-1mm. with the aid of a hand lens magnifying about | 
five times. The perpendicular distances p of the medians from O, are also|f 
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measured, usually to 0-5 or 0-3 mm., or if possible, to 0-1 mm. in cases where p is 
only a few millimetres; and a table of D and p values is constructed from which 
D . f(p/L) and (p/L)D . f(p/L) values can be calculated. 

Where possible D should be derived from measurements between the sharper 
higher orders, divided by the order number. If it is necessary to measure between 
faint lines, the average of several estimates should be taken. Note especially that 
when the main line pairs have been sorted out the sharpness of any intensified line 
can be used to give an indirect estimate of D and so usually leads to identification 
of the practically parallel absorption line (with the help of a lens), even when this is 
too faint to be noticed immediately. 

If D cannot be measured directly because the region near the foot of the per- 
pendicular from O, to the median lies outside the recorded photograph, the 
spacing D,, must be measured where possible at a distance y from this region, and 
D estimated by use of figure 4(b), as described in §3 (i). 


Construction of the reciprocal-lattice projection. It is desirable to mark points 
on both sides of the origin, which is chosen some 20 cm. or more to one side of the 
print, and the geometrical positions of several higher orders should be marked 
even if they are too faint to be clearly distinguished in the pattern. The serial 
number of the net-plane is marked below each order, with the order number as 
suffix if desired. Accuracy of location of the radius-vector normals to the line- 
pair medians is essential and is assisted by careful use of transparent 60° set-squares 
of 10-inch hypotenuse. Allowance has to be made for the appreciable opposite 
curvature of the components of a line-pair in estimating the direction of the median, 
by finding the direction midway between chords or tangents in the same region of 
the line-pair. 

The final stage is to write above each point in the projection the calculated 
distance of the point from the projection plane, with the appropriate sign. The 
projection is to be considered always as viewed from below the recording plane, 1.e. 
in a direction opposite to the primary beam, so that the positive z distances (along 
OO,, as defined in §3 (i)) are from the projection plane towards the observer, and 
negative distances are on the opposite side. With the sign convention stated 
above and if a positive print is used, the projection points whose radius vectors are 
in the same direction as that in which p is measured from O, will have negative 
distances from the projection plane (cf. figure 3). 

For the purpose of determining the lattice shape, dimensions and orientation, 
it is desirable to mark the reciprocal-lattice points as small points surrounded by a 
ring, but it will often be useful to make copies showing the points as black discs 
with diameters a measure of the intensities of the diffractions. 


(iv) Use of scales to read off AL .r,,* and XLz directly from line positions 


When the line-pair separations D are mostly between | and 2cm., and the 
patterns not too intricately criss-crossed, time and trouble can be saved with but 
little loss in accuracy by the use of measuring scales from which D -S(p/L), i.e. 
AL .r,* and D(p/L) . f(p/L), i.e. ALz can be read off directly as the difference in the 
scale readings at the black-and-white-line Bragg positions when the scale is placed 
normal to the lines with its zero at O,. 

The possibility of constructing such scales is due to the fact that the factor 

f(p/L), by which D is multiplied to obtain AL . 7,,*, is not far from unity in all cases. 
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likely to be needed, and that its variation is small over a range of p/L of approxi- 
mately 0-1 corresponding to band widths D up to 2cm. when L~20 cm. Under 
these conditions if a function, y, of p is derived such that 


dy|dp =f(p/L), 
we can take as a close approximation 
Ay=Ap.f(p/L), or (Y2- 1) = (P21) - (Pi + P2)/2L5 
even if py —p, is made so large as to take the value of the line-pair separation D, so 
that 
2-1 =D . f(D + Py)/2L}. 
Correspondingly, to construct the scale to read off AL. x directly, we require a 
function y’ given by 
dy’ |dp =(p/L) . f(p/L). 
We have thus: 
f(p|E)= cos" = [1 + (p?/L2)P% 

y=\f(o/L) dp =p . fol p/L) =p — pt —fal0/L)}, 

y =\(/L) . f(p/L) dp = Lt —fe(p/L)}.- 
The values of f2(p/L) calculated are given in table 1, whence values of (p —y) and 
y’ were calculated corresponding to a range of p values from L=15, 20, 25, 30, 35, 
40, 45 and 50cm. From these data scales were drawn which were then photo- 
graphically reduced to the required size. A typical scale is shown in figure 6. 
Similar scales for any required intermediate value of LZ can be constructed if 
necessary in the same way. 


Table 1. Variation of 1—/.(p/L) with p/L 


p/L 05 il 15 2 25 3 35 “4 45 3) 
1—f,(p/L) :00125 -0050 -01106 -0194 -0299 -0422 -0561, -0715 -0880 -1054 


A trial of such scales over a wide range of p and D showed excellent agreement 
with the calculated values of D. f(p/L) and (p/L)D.f(p/L). The best scale for 
practical use is a reproduction on a photographic plate, but it is also possible to use 
a paper print of the appropriate scale, stuck on the underside of a piece of trans- 
parent plastic sheet, if the shrinkage of the paper after processing is found to be 
negligible. The zero is placed at O, then the scale moved slightly if necessary so 
that a scale mark comes opposite one Bragg position, whence only a single inter- 
polation to a tenth of a division is necessary at the other Bragg position and the 
difference of the two readings is easily found and is then written down in the list 
of pattern data. When the list is completed the projection can then be constructed 
as described above. 


§4. CHOICE OF AXES, DETERMINATION OF SYMMETRY AND ORIEN: | 


LATION, INDEXING OF, PATTERN AN DECAL GUD Awl ONE OIE 
LATTICE CONSTANTS 


When an adequate number of points are located in the projection and the 


distances of the lattice points from the projection plane are marked, it is easy to | 
recognize the main reciprocal-lattice rows, because any straight line in space is also _ 


represented by a straight line in the projection and, further, the distance ALz of the 
lattice points from the projection plane must be proportional to their distance along 
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the line in the projection. The rows with densest lattice-point population are of 
course the most prominent. Similarly, it is easy to recognize which are the main 
lattice planes because they contain two or more systems of parallel densely-popu- 
Jated lattice rows. As a convenient ‘scaffolding’ or reference of coordinates to 
assist in this process, it is helpful to draw lightly the main lattice lines through the 
coplanar set of points (forming a two-dimensional lattice) due to the planes of a 
main zone, usually one whose zone axis is near the primary beam direction, so that 
the points in the reciprocal lattice plane correspond to the Laue spot pattern in the 
circular spot zones. 

Having tentatively drawn in pencil some coplanar sets of main lattice rows in a 
favourable part of the projection, so as to show part of the cross-grating set of 
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lattice points in such planes, a different coplanar set can be explored and delineated 
in the same or another part of the diagram. From such initial trials it is usually 
clear whether angles of 90° or 120° occur, and which is the appropriate crystallo- 
graphic system of axes and lattice symmetry, and the whole diagram is then mapped 
out in lattice rows and planes accordingly, so as to outline the neighbouring unit 
cells of the reciprocal lattice; the indexing of the lattice points then follows easily 
as the simple process of writing down the coordinates of the points relative to the 
chosen axes. 

From a consideration of the indices of the observed diffractions, supplemented 
by similar data from photographs obtained at other crystal settings, any systematic 
extinction of definite types of diffractions can be recognized and used in the usual 
way to determine the possible space groups characteristic of the atomic arrange- 

nent in the crystal. Allowance must be made, as in x-ray diffraction, for the fact 
that the ‘Laue symmetry”? shown by diffraction patterns includes a centre of 


354 H. Wilman 


symmetry even if there is none in the symmetry of the crystal structure; and other 
data such as observation of crystal forms and etch figures is required to decide this 
point. In electron diffraction, as already noted, there is also the possibility in 
certain cases that orders of reflection from a net-plane may occur with appreciable 
intensity, due to secondary or dynamic scattering, which would otherwise be 
forbidden according to the more approximate single-scattering theory. The 
further consideration of Kikuchi-line intensities will be left to a later contribution. 

The angles between the reciprocal-lattice rows in planes inclined to the main 
projection plane, and the spacings between the points in the rows, can be estimated 
by the usual methods of projective geometry. Side projections in suitable direc- 
tions, on planes perpendicular to the initial projection plane are especially useful 


yocm. 


0 5 
Scale of Actual Drawing and its Vertical Co-ordinates 
Figure 7 (0). 


when these can be drawn so as to give desired angles and spacings directly (see 
figures 7(b) and 8(b)). In awkward cases it may sometimes be most direct to 
construct the lattice-point positions in space by mounting the projection on 
cardboard and sticking pins through the projected points so that their heads are at 
the appropriate heights. For convenience all the points can be brought on to one 
side of the cardboard base by adding a constant height to the calculated distances 
of the points from the projection plane. 

The reciprocal-lattice angles and the lengths proportional to its axes are 
estimated from the projection or the model by taking averaged directions and 
lengths most consistent with the constructed lattice point positions over a region 
comprising as many unit cells as possible; the actual reciprocal axes are then 
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Figure 7(a). NaCl, polished and etched face; L~22 cm., 50 kv. 
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Sulphur conchoidal fracture; L~22 om., 54 ky, 
of 4 and 37 band medians. 


Vigure 8 (a), ; O, at intersection 
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calculated by dividing the above lengths by the scale factor of the drawing or model 
and by the estimated AL for the photograph. From the a*, b*, c*, «*, B*, y* s0 
obtained, the crystal lattice constants a, b, c, a, B, y can be calculated front the 
usual relationships. ‘The symmetry of the crystal lattice is of course of the same 
kind as that of the reciprocal lattice. 

Although the determination of the lengths a, b, c involves the value of AL, 
which is not usually known with great accuracy, the axial ratios and the angles «, 
B, y can be obtained fairly accurately as above from the reciprocal lattice measure- 
ments, or also from the angles between net-planes by the usual methods used with 


| 
focm. 
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Figure 8 (2). 


goniometric measurements obtained from crystal faces. In any case, measure- 
ment of the angles (180°—4,,) enclosing O,, between suitable line-pair medians 
in the pattern, can be used to calculate the angles 6 between the corresponding 
net-planes, as check comparisons with the 5 values calculated from the lattice 
constants obtained from the reciprocal-lattice projection. The angles (180°—%,) 
can usually be estimated to within about 0° 10’ by use of a normal set-square and a 
transparent rotating-arm set-square having a 4-inch radius scale of degrees. To 
obtain 5 from the measured angle (180°—3,,) and the p, and p, values for the two 
line-pair medians, we have simply 


‘cos 5, = {cos 6 — cos (90° — 1) cos (90° — z)}/sin (90° — d;) sin (90° —4,), 
go that cosd={cos5, + (PipelL*)} - {1+ (pilL)}> «(1+ (Pal L)2}* 
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§5. EXAMPLES OF RECIPROCAL LATTICE CONSTRUCTION 


Two examples are given to illustrate the results and the accuracy obtained in. 


typical cases and to show some features of the line intensities. In these examples, 


since the camera-length and position of O, were not known very precisely, the 


accuracy obtained in the projections is not the maximum obtainable, but even so, 


there is no difficulty in recognizing the periodic distribution of the lattice points. 


and the lattice symmetry. 


(i) Rock-salt, NaCl, figure 7 (a) 


This relatively simple pattern was obtained from a rock-salt crystal face which 
had been cut in an arbitrary direction then polished and lightly etched with water. 


Lines from about fifty planes were visible, some with strong second and | 


higher orders, and as the pattern has a plane of symmetry denoted by the line 


m—m, corresponding planes were numbered similarly, one set with distinguishing | 
dashes. (Since the rock-salt structure is known to have (100) and (110) types of | 
mirror planes the plane ‘“‘1”’, i.e. m, must be of this type, and the analysis shows | 
that it was a (110) plane.) The most prominent zone axis where many of the | 
narrower bands cross, slightly below the centre of figure 7 (a), was used for con- | 
venience as O,, though not a very close approximation to O,, and the reciprocal- | 
lattice projection shown in figure 7 (b) was constructed from the D and p measure- | 
ments. ‘The projection, figure 7(b), shows the neighbouring body-centred cubic _ 


cells of the reciprocal lattice, and the side elevation on the right demonstrates the 


rectangularity of the lattice. Probably a better consistency of the coordinates of | 


the points in the outer parts of the projection would have been obtained with the 
closer estimate of O, at about 2cm. nearer the undeflected-beam spot, and in 


general it is advisable to use the actual estimated position of O,.. The Laue indices 


of the line diffractions relative to the axes chosen are shown in table 2. 


The line intensities depend on various factors, including the structure factor, | 
S’, which is 16(Na + Cl)? if h, k, Jare all even, 16(Na— Cl)? when A, &, /are all odd, | 
and zero when A, k, / have other values; h, k, ] are the Laue indices and Na and Cl | 
denote the atom-form factors. As a first approximation the atom-form factors _ 
can be considered to be represented by the atomic weights multiplied by a function | 
which decreases with increasing deviation of the diffracted beam from the primary | 
beam direction; thus it is natural that the line intensities should decrease with | 


increasing D for both groups, i.e. h, k, J all odd and all even. As an example, in 
the orders of diffraction from 111, 111 has medium, 222 strong, 333 nearly zero, 
and 444+ medium intensity. 

The asymmetrical position of the mirror plane m relative to the primary beam 


enables the relationship between the line intensities and the inclination of the | 


reflecting planes to the beam to be demonstrated clearly in a single reflection 
pattern by many pairs of band systems. It is possible in this way to see that the 
expression for the line intensity must contain a factor which reduces the intensity 
rather rapidly as the angle between the plane and the primary beam (roughly 
proportional to the perpendicular from the undeflected-beam spot to the band 
median) approaches zero. ‘Thus, for example, the difference in intensity of the 
lines of 8 and 8’, 7 and 7’, 10 and 10’ etc., is not large, while in 3 and 3’, and 4 and 4’ 
it is larger, and in 5 and 5’, 9 and 9’, and 11 and 11’ it is quite high because the 
planes 5, 9 and 11’ are so much nearer to the primary beam than 5’, 9’ and 11 


————— 
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respectively. Although the more strongly developed narrow bands almost dis- 
appear near the undeflected-beam spot, as in the case of 5, the intensified lines 
(white in figure 7 (a)) of wider bands are often strengthened near the foot of the 
perpendicular from the undeflected-beam spot as in 24’, 25’ and 26, apparently 
because of the narrowness of the selective-reflection range (corresponding to the 


line sharpness) in such cases, and the steep gradient of the background intensity 
near the undeflected-beam spot. 


Table 2. NaCl line indices and intensities in figure 7 (a) 


ee ee ee 


Laue Laue 
Bond indices Intensity pang indices Intensity 
sa hil a hil 
1 230 VS —_ = = 
440 M 
: 660 F 
2 111 M == = ae 
223 Ss 
333 ~mnil 
444 M 
3 234 ME 3” 524 MS 
4 042 MS 4/ 402 VS 
084 VE 804 F 
5 131 MF 5/ 311 M 
262 M 622 M 
6 443 Ss = a 
fi 062 MS of. 602 S 
8 423 S 8’ 243 S 
844 VVE 484 VVE 
9 282 VE 9’ 822 MF 
10 420 S 10’ 240 MS 
840 F 480 VVF 
11 640 M ile 460 VE 
12 462 M 127 642 M 
13 242 Ss 13% 422 M 
14 824 F 14’ 284 F 
15 464 MF 15’ 644 MF 
16 311 MS 16’ 131 M 
622 MS =— 263 M 
17 620 M = gts ee 
18 331 ME — — — 
662 EF 
19 513 F 19’ 153 F 
20 531 a == a 7a 
21 802 MF Bile 082 
22 51l F — = = 
23 842 VE Day 482 
24 064 F 24’ 604 M 
25 244 F 257 424 M 
26 646 F 26/ 466 F 
27 664 F == =< = 


S=strong; M=medium ; F=faint ; V=very. 


A further prominent feature of most Kikuchi-line patterns, conspicuous in 
figure 7 (a), is the peculiar weakening or step-down in intensity of certain sections 


ea 
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of some lines; for example both the white and black lines of 3’ almost disappear 
between the second-order band of 1 and the adjoining first-order band of 4 (or 
second-order band of 5), while the line-pair of the symmetrical plane 3 behaves 
similarly across 1 and 4’ (not 4). Such conditions help to make the parabolic 
line envelopes more prominent, for example that (concave to the left) formed by 
the white lines of 3’, 4’ (first appearing order), 5’ (first-order), 1 (first-order), and 
the corresponding black lines of 5, 4 and 3, and the mirror-symmetrical parabola 
(concave to the right) formed by the black lines of the first group and the white lines 
of the second. 
(ii) Sulphur, figure 8 (a) 

This pattern was from an irregular conchoidal fracture of a sulphur crystal of 
about 1cm. diameter. The bands are nearly all narrow and correspondingly | 
diffuse-edged, but by applying the criteria described above for the Bragg reflection - | 
positions it was possible to measure the band widths to an accuracy of 0-1 mm. and 
hence to construct quite satisfactorily the projection of the reciprocal lattice, 
figure 8(b). Here again it is easily seen that the pattern has a plane of symmetry 
m represented by the median of band 1, with respect to the band positions and 
intensities; but band 2 represents a plane which is at right angles to 1 but is not a 
mirror plane. 

The lattice-point projections in figure 8 (0) lie close to a series of equally-spaced 
lines parallel to AB. They will not all lie exactly on these lines because the crystal 
is such that the lattice planes of (0R0) type are not quite normal to the plane of pro- 
jection, i.e. the photographic plate ; in this case, it is sufficiently close for indexing 
purposes to use these lines without adding all the closely neighbouring lines which 
would strictly be required in the projection. ‘The points are also grouped on or 
near lines parallel to CD, which represent projections of lattice rows inclined 
downward from left to right at about 84°, as the z coordinates of the points show. 
The side elevation represents the points of intersection of all these lattice rows with 
a plane through O at right angles to the projection plane (represented in section by 
EF) and to CD, and shows the centred rectangle distribution decisively. Strictly 
the projection of the lattice rows should be done upon a plane which is normal to 
the rows, not to their projections, but the small inclination of the plane of the 
illustrated side elevation to this plane introduces little error. The data in figure 
8 (6) thus show that the Kikuchi-line median positions and the line-pair separations 
correspond to a face-centred rhombic unit cell, and the ratios a:b:c estimated 
by calculations based on the measurements of the projection combined with the z 
coordinates were 0-82:1:1-87, Using the approximately known AL value of 
1-13 a.cm. a, b and c were found to be 10-6, 12-9 and 24-2 a. respectively. This 
is in good agreement with the previous x-ray results, since Mark and Wigner (1924) 
found rhombic sulphur to have a face-centred rhombic cell with a=10-61 a., 
b=12-87a., c=24-56a. thus a:b:c=0-820:1:1-900, with space group V2! 
(Fddd) and 128 atoms per cell; Warren and Burwell (1935) who confirmed this 
space group and calculated the atom positions, found a=10-484., b=12-92a., 
c=24-55 A., i.e. a:b:c=0-811:1:1-900. Trillat and Oketani (1938) concluded 
from their electron-diffraction transmission patterns that a=5-2a., b=6-4A., 
c=12-2a., but these values do not take account of spots which occur in some 
regions of the patterns and show the true Laue-zone spacings, which corre- 
spond to cell dimensions close to those of Mark and Wigner (1924) and 
Warren and Burwell (1935). 
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‘Table 3 shows the Laue indices of the observed lines in figure 8 (a) and gives, 
where available, the x-ray structure factor and spot intensity in the rotation 


pattern, as listed by Warren and Burwell. Unfortunately only eight diffractions 


Table 3. Sulphur line indices and intensities in figure 8 (a) and x-ray data 


iba eave Si Gk 
Gl eke x ray data er pee x ray data 
a3 indices eat Struct. int indices Int Struct. Int 
hkl factor hkl = factor a 


1 040 VS 472 S (29 44,12 M 
080 MF 3 nil J (113) MF 
0,12,0 VE 30 4,4,12 ME 
0,16,0 VVE i! (103) ee} 

2 206 VS 440 VS (31 44,10 M 
OIG 1a) < 8,8,20 F 

3 408 M 427 MS (32 4,410 VF 

4 131 M 33 .4,2,10 M 
262 MS ae 4,2,.10 MF 
393 M 85 357 VS 
412540 O ibe 357 S 

5 175 S 37 6,2,14 MS 

6 175 8 403 MS tee 6,214 8S 

7 Ose? vag) [39 488 MS 

Sa On1022) 2k 40 488 MS 

9 062 Ss 610 Ss 41 395 M 
Ono eh} 42 755,15 ME 

3 56,3511 ™ 


BWA yy 
| & 
Z 


11 5 6,0,14 M 
12 ; 
13 Nil F dip ST) AP 
14 val F sf) UAL 48 6,6,10 M 
15 284 MF ay Seillge AR 
16 284 MF BOM TON lO mmmneE! 
17 268 ME 51 135 MS 
AGUS ame 2085 0° -F f52 -224 M 
1) ALR Te 53 224 F 122 F 
Oe Ae 14 VE 54 337 MF 
f21 359 MS [tbe = 105) ee 
“| 22 359 ME 56 195 F 
Sf 28 133 F 57 355 M 
e224 133 FE Re Sei 
ay AAI 59 466 MF 
26 4,212 F f60 438 $M 
27 319 MS {61 428 F 
Ae 319 MS J 62 193 EF 
(63 193 F 


are common to both lists, but these are all in good agreement except in the inten- 
sities of successive orders from the (010) plane, where (080) is evidently much 
stronger in the electron-diffraction pattern than in the x-ray pattern, presumably 
due to appreciable secondary or dynamic diffraction. Figure 8(a) also shows 
further examples of the great reduction in intensity of bands which are parallel to 


the primary beam or nearly so. 
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§6. SUMMARY 

The interpretation of electron-diffraction Kikuchi-line patterns from single 
crystals is developed as a powerful, independent and general means of determining 
the crystal lattice type, dimensions and orientation and the Laue symmetry of the 
atomic arrangement. This aim is realized by a construction of the reciprocal 
lattice by methods which are described and illustrated. ‘The satisfactory accuracy 
of the method is due to the relatively precise definitions for the location of the 
Bragg reflection positions, derived and stated here for the first time, which 
enable the method to be used even when the crystal lattice constants are so large 
that the bands are narrow and diffuse-edged, as with rhombic sulphur. The 
patterns shown also illustrate features of the line intensities, especially the reduction 
in intensity when bands pass near to the undeflected-beam spot. 
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Thermionic Emission Constants and Band Overlap 
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ABSTRACT. A general treatment is given of the effect of band overlap on thermionic | 
emission. ‘The results are applied to nickel, for which the experimental value of the current _ 
constant A of 1380 amp. cm~* deg~* differs markedly from the standard value of 120 amp. 
cm * deg~*. The case of nickel has recently been discussed by Sun Nien T’ai and Band | 
(1946), who derive an expression for the variation of A with temperature. It is shown 
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that their relation strictly applies only to ‘“‘ holes ” in an isolated d-band and is reliable only 
over a temperature range such that kT/¢9q is considerably less than unity, ¢9q being the 
unoccupied width of the d-band. With an overlapping s-band, as in nickel, there is an 
increasing transfer of electrons from the d-band as the temperature rises. The effect is 
strongly dependent on the ratio r9= €os/€9q, where €os is the occupied width of the s-band. 
Expressions are derived and precise calculations carried out for the transfer effect and the 
variation of A over a range 0-0 <kT/ €9q < 1-0 (including temperatures up to about 2000° k.), 
and for ry values of 10, 20 and oo. The last corresponds in effect to an isolated d-band. 
The results of the calculations are shown in figures 2 and 3. A comparison with the 
experimental value of A is made using a value of ¢ q estimated from thermal and magnetic 
data. The satisfactory agreement previously obtained is due to the adoption of a value of 
the bandwidth which is incompatible with the experimental evidence, and it is shown that 
the calculated and observed values of A actually differ by a factor of about 6. Various 
reasons for this discrepancy are discussed, especially the effect of surface impurities. 


§1. INTRODUCTION 
N a recent paper Sun Nien ‘T’ai and Band (1946) attempted to interpret 
anomalous values of thermionic emission current constants on the basis 
of the electron energy band theory of metals. They consider the effect 

of band overlap in the vicinity of the surface of the Fermi distribution, and show 
that under certain circumstances this may lead to abnormally high values of the 
current constant. In the actual treatment, however, no account seems to be 
taken of the possibility of a transfer of electrons between the two overlapping 
bands. The relations obtained are simply those for separate single bands, 
in each of which there is a fixed number of electrons (or “‘holes’’) determined 
by the degree of overlap. ‘The calculation differs essentially from that for a 
standard single band only in that the energy density of states is much greater 
and the width of the unoccupied part of the band much smaller than for a “ free- 
electron”? band containing a corresponding number of electrons. ‘To the 
approximation used, the expressions given appear to be correct, but in substituting 
a numerical value for the d-band energy width for nickel, the value used differs 
by a factor of about 10 from that which the most direct interpretation of the 
available experimental evidence seems to indicate. The agreement between 
the observed value of the current constant A for nickel (1380 amp. cm? deg? 
as compared with the standard theoretical value of 120 amp. cm~? deg?) and 
that calculated would thus appear to be largely fortuitous. For this reason, 
and also because the ideas presented have a bearing on many other physical 
properties of metals, the problem of thermionic emission constants and band 
overlap is reconsidered in this paper, with special reference to nickel. 

Variations in the value of A, the current constant occurring in the Richardson 
equation, may be due to a temperature dependence of the thermionic work 
function y (Fowler and Guggenheim 1939, p. 484). ‘The Richardson equation, 
written in the form 

I=A,(1—7)T? exp (—x/RT), = are (1.1) 
gives the variation of the emission current per unit area with temperature, 
7 being the reflection coefficient, which will be very small for most clean metal 
surfaces, and A, being a universal constant having a value of 120 amp. cm™ deg™?. 
Hence if x» is the value of x at absolute zero, and A the value of the current 
constant at a temperature 7, then 

PA 7) de exp (—xyRL),- wees (1.2) 
where A AgeseXp (Vax) 9 0) Baeees (i 
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Usually the experimental results are expressed in the form (1.2), with y, constant, 
so that a decrease of the work function with temperature will be manifested by 
an increase of A with 7. Such a variation will occur for the band form proposed 
by Mott (1935) and Slater (1936) for nickel. In this assembly the electrons 
are distributed in two overlapping energy bands, one of which, the “d-band”’, 
corresponding roughly to the 3d electrons, is very much higher than the other 
one, the “‘s-band”’, which corresponds to the 4s electrons, and which will approxi- 
mate to a “free-electron” distribution (see figure 1). Fermi-Dirac statistics 
are readily applicable to this assembly if for the relevant portions of each band 
v(e), the number of distinguishable states per unit energy range, is proportional 
to the square root of the energy «. The diagram given by Slater on the basis 
of numerical calculations indicates that this assumption is justifiable for nickel 
as giving a reasonable approximation (see§4). For other metals there may 
be considerably greater deviation from the standard parabolic band form. 

For nickel the equilibrium state at y@ 
absolute zero is such that the number of 
electrons in the s-band is equal to the 
number of absent electrons, or ‘‘holes’’, 
in the d-band. Evidence drawn mainly 
from a study of the magnetic properties 
of nickel and its alloys indicates that at . <3 € 
absolute zero approximately 0-6 electrons cS 
per Ni atom are in the s-band, and 9-4 Figure 1. Collective electron assembly for 
electrons (i.e. 10—9-4=0-6 holes) in the ee 


: . i i 7. , density of 
Aiband  aeehou €, electronic energy. v( ©), energy density 
; ow Dn figure 1. It is states. € 9s, occupied width of the s-band;, 


through the presence of these holes in the éyd, unoccupied width of the d-band, both 
d-band that most of the thermal and mag- at absolute zero. The diagram is illustra- 
netic properties of nickel have received tive only (cf. Slater 1936). 
satisfactory interpretation. 

€o, and «9, are the respective widths of the relevant portions of s- and d-bands 
at absolute zero. For nickel 


| . v,(e)de = 


0 ~ gs 
at all temperatures, and it is seen from the diagram that 


Re sT, >i ee (1.5) 


As will be shown below the equilibrium appropriate to absolute zero will be 
displaced as the temperature increases, when there will be a transfer of electrons 
from the d-band to the s-band. As a result the number of holes in the d-band 
and the number of electrons in the s-band both increase with a rise in temperature. 
A method of estimating the extent to which electrons are transferred in this way ~ 
is given in §2, and it is shown that the current constant as calculated by con- 


sidering an isolated narrow d-band has to be modified if the transfer effect is taken 
nto account. 


d-band 


s-band 


rEgs t £od 


vi(ede = Ap An eee (1.4) 


§2. THEORY OF THE TRANSFER EFFECT IN NICKEL 


In this section is given the quantitative treatment of the electron transfer 
which takes place in nickel. In the analysis the notation employed is in the 
main that proposed by Stoner (1939 a), whose paper should also be consulted 
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for fuller details regarding Fermi-Dirac statistics. The distribution expression 
applicable to free electrons, or, as here, to the collective electrons in parabolic 
bands, may be written in the form 

dN 3 Nj (e?/€92) 

de 2exp(e/RT—n)+1° 


Here N, is the number of electrons per atom present in the band and « its occupied 
width, both at absolute zero. Alternatively N j may be the number of holes 
per atom, in which case ey refers to the unoccupied portion of the band. The 
statistical parameter, 7, is related to the chemical potential, or mean Gibbs free 
energy per particle, u, by 


iene Se eg) ee eee: ere (2.2) 
On transforming (2.1), putting x=«/RT, integration over all positive values of 
gives 
kT\? (* xPdx 
N=3N (=) | ees eaeel(Q3 
aN ot) 6 CRP) oe 
Cr (Rige, ha. 9 |e Saas (2.4) 


F(m) is the Fermi-Dirac integral given in (2.3). The relation (2.4) states 
essentially that, as the temperature T increases from absolute zero, the total 
number of particles per atom also changes, unless 


ATMEL IAG) Sik oa alee We) See (2.5) 


This relation holds for a single band, not overlapped by any others. In the 
present case it is satisfied in the limit 7)=€9,/egg>00. The total number of 
electrons in the band is then necessarily constant. For the assembly to be 
considered, however, with 7 finite, the relation (2.5) does not hold for either 
band. 

For the s-band let N, and No, be the number of electrons per atom in the 
band at temperature 7 and at absolute zero respectively. Also let «9, be the 
occupied width of the band (cf. figure 1) and RT, the chemical potential. Then 
from (2.4) 

Nie eNew ie 1ien.) EF ae)ual) ao plat aeons (2.6) 


A similar expression is applicable for the d-band, but here Nz and No, are the 
respective numbers of holes per atom. Let «9, be the unoccupied width and 
kT», the chemical potential as measured downwards from the top of the d-band. 
Denoting the chemical potential of the d-electrons, referred to the same level 
as that of the s-electrons, namely the bottom of the s-band, by RT'7’z, one obtains 


RI 5 =(oe oq) Riga ee 27) 
The relation corresponding to (2.6) is then 
Na=*NodRTleoa)E (qa)? te (2.8) 
For nickel the condition (1.4), which holds at all temperatures, gives 
Ngee Nog Vo ee 2-9) 


The condition of equality of the chemical potential (i.e. of the Gibbs free energy) 
for thermodynamic equilibrium gives 


n'a => Ue eccece 
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Insertion of the condition (2.7), (2.9) and (2.10) into (2.6) and (2.8) finally yields 


the relation 


F {(€oa 4 €9;)/R [ie natl(€os)® = F(1a)/(€oa)?s 


or, writing 7p =€o./€9q and Tp =RT/eoq, 


F{(1+79)/to—na} =To'F a). eee (2.11) 
The number of holes in the d-band, Nj, is thus given by (2.8) as 
NalNoa= 270 F(na); Bae 0-2) 


nq being determined from (2.11). Owing to the conditions (2.9), the ratio 
N,/Nos is equal to Nz/Noq and is thus also given by (2.12). It will be shown 
below that as T increases the ratio N,/Noqg=N,/No, increases from unity at 
T=0. With a rise in temperature the d-band thus loses electrons which are 
transferred to the s-band. During this process the thermionic work function 
also changes, and it easy to show from the definition of this quantity (Fowler 
and Guggenheim 1939, p. 477) that 


(Xo- x)= —(RT7q—€oa) =(RT75— 0s): BOIS (2.13) 


A/Ag= exp (79* —7a)- (2.14) 
Precise numerical calculations of Nj/Noqg and of A/Ao, using (2.12) and (2.14), 
have been carried out for a series of temperatures (§3). ‘These may be usefully 
supplemented by series expansions appropriate for t)<1. For 7)<1, ng>1, 
{(1+79)/7) —nq}>1, the asymptotic series expansion of the Fermi-Dirac integral 
may be written in the form 


F(a) 7 ang{l =. (7?/8)nq? ar ood iS 


~ 


Similarly, and to the same degree of approximation, writing 
2={((1+71)/7o— na} 
F(z) = 324 {1 + (?/8)2-}. 
Hence from (2.11), after some rearranging, 


Pe he mal rolTo— Na 
” “12 Lyall +79)/79 — na} _| 
To a first approximation the right hand side of (2.15) is zero, so that nyg=7)7 


=€oq/RT. Inserting this value into the right hand side the second approximation 
is obtained in the form 


Hence from (1.3) 


aye? ih = (a7? 79/12)(75 = 1 ines eee ene (2.16) 
Hence from (2.8) and (2.9) 
Nal Noa=NolNos=1 £07? /8%o\ro*. ee (2:E7) ae 
The relation (2.17) shows that for 7>>0, Ng>Nog and N,>N,,, so that there is _| 
in fact a transfer of electrons from one band to the other. The temperature ] 
variation of the current constant is given by (2.14) and (2.16) as 
. AjAy=expaar* Ty 12) 6 — 1) Foe (2.18) || 
Since T> will in most cases of physical interest rarely exceed unity (see § 4) 
it 1s unnecessary to derive any high temperature approximations, applicable : 
if 74>1. A point of interest does, however, arise in connection with (2.11)-41 
For 7) 00 this relation defines a certain lower limit of Na Which is here denoted 
by y, where | 


F(—y)=79 FO). ae eee ee (2.19) 
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For ry>1, y will be negative and will be a function of 7) only. The expression 
(2.14) now shows that the limiting value of A/A) (for ty) is very simply 
related to y by 


AmaxiAo=€Xp|¥l, save (2.20) 


Ajax being the maximum value of the current constant for any particular value 
of 79. Values of y and A,,,,/A, have been calculated for 7)=10 and 20 and for 
ro 0 and are given in § 3 as supplementing the curves shown in figure 3. 

As regards the relations derived by Sun Nien T’ai and Band (1946) it may 
be noted that, since t)=RT/e,, the low temperature expression (2.18) is exactly 
equivalent to their relation (13), which, assuming that 7,>1, was further 
simplified to give 


AlA,=expiGr/IZ)\RTlega)pe = wes (2:21) 
‘This expression was taken to be the basic relation and was employed for 
“complete calculations of current-temperature curves’’. Such a procedure is, 


however, open to objection. In the first place the low temperature approxima- 
tions used are reliable only if RT/egg<1. For nickel, however, ¢94/R is of the 
-order of 10? °K. (see § 4), so that, at the higher temperatures used in experimental 
determinations of the thermionic constants, RT/e),; may be of the order unity. 
For this reason precise calculations of the type outlined in § 3 have to be carried 
out. Furthermore, by stating their relation (14) in the form given, Sun Nien T’ai 
and Band seem to neglect the effect of the s-band altogether. This relation 
‘would in fact be obtained by considering the temperature variation of (¢)/RT —7) 
for the d-band by itself, the transfer effect being thus neglected. A second 
approximation to the Richardson equation had, in fact, already been derived 
in this way by Fowler and others (see Reimann 1934), although, since particles 
_-rather than holes were considered, a decrease of A with temperature had been 
predicted. 


§3. COMPUTATIONAL DETAILS AND RESULTS 


The equations to be solved are given in §2 and are here rewritten for con- 
venience :— 


F{((1+79)/to—nat=To F(a), wt eee (al) 
NalNoa = 37°F (nq), ccapererere (322) 
AjAg=exp {t> —Na}- wee (3.3) 


_As has been stated the approximations (2.17) and (2.18) are applicable only if 
-t)<1. Since, however, this condition is not satisfied if the temperature is high 
and if €9, is relatively small, it may be necessary to obtain exact numerical solutions 
of the equations for a suitable range of 7), say 0-0<7)<1-0, and for one or more 
‘values of 7). For nickel the magnitude of this parameter is rather uncertain. 
This limitation of the present treatment is discussed further in § 4, where estimates 
-of €9, and €9q are obtained. These indicate that 7, lies between 10 and 40. ‘The 
-calculations in this section were carried out for two values of ry, namely 10 and 20. 
Additional calculations for 7700, corresponding to an infinitely wide s-band, 
-can be made very simply, and the diagrams showing the results of the three 
sets of calculations almost certainly cover the physically relevant range. 

It would be most convenient to obtain values of N,/No q and of A/A, for 
-equally spaced values of 79 lying in the desired range. For the present purpose 
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it was found to be simpler, however, to proceed as follows. For a given value 
of rq a trial value of 7, is first chosen, satisfying the condition that na>y, where | 
y is given by (2.19). The magnitude of the corresponding Fermi-Dirac integral | 
F(7,) is then obtained from the tables of McDougall and Stoner (1938), so that 
F{(1+79)/7)—n} is calculable from (3.1). To obtain values of the argument 
inverse interpolation in the same tables would be necessary if 7) were small. 
For 7, equal to 10 or 20, however, the left hand side of (3.1) is large, and the value 
of {(1+79)/79 —7q} can be obtained from the inverse asymptotic series, which is. 
discussed fully by McDougall and Stoner. With 


a=[§F{(1+79)/79— naj} *= ro [3F (nal? 
it is shown that 
2 


{(1+179)/7p> —ngt =a Hl — aya? —ayat—agae—....}, sae ae. (3:4) 
where a, =0-822 467, a,=1-217 614, ag=9-161 386. 


For the purpose of the present calculations the three terms of the series. | 


given in (3.4) were found to give a sufficiently high accuracy. Having obtained 
{(1+79)/7) —na} from (3.4), nq is calculated using the same value of 79, 10 or 20, 
as before. If this value of 7, lies outside the range to be covered (0-0<7)<10), — 
a different value of y, has to be chosen and the calculation repeated. Witha jf 
suitable 7,, t) having been obtained by the method outlined, Nz/Noq and A/Ay, |} 
are determined from (3.2) and (3.3). About 10 values of 7, lying in the desired 


range were calculated for r=10 and 20 respectively, the computations being |} 


carried out with the aid of a Brunsviga calculating machine, using Barlow’s. 
Tables (ed. Comrie 1941) for the roots and the Federal Works Agency Tables. 
(1939) for the exponentials. ‘The procedure may be illustrated by giving the 
details for 79 =10, n7=1-0, the numbers being given to five significant figures.. 
Then 
F(q) = 1:3964, 7 =0-63510, 
F{(1+79)/t) —q}=44-157, Na/Noa=1-0601, 
{(1+79)/79 —ng}= 16-320, A/A, =1-8766. 
The final results are most conveniently given in graphical form and are shown in: — 
figures 2 and 3. Besides the results for 7) =10 and 20, those for ry+00 are also: | 
indicated. ‘These were calculated using the relation (2.5) which gives 


n= GF) 
NojNoa== lL, e Sa” g th Soe (35) 

A/Ay= exp [{3F (ya)}8 —Hal 
Figure 2 shows the variation of N,/Noqz with temperature. As is to be | 
expected from the relations (2.5) and (3.2), the transfer effect will decrease in | 
magnitude as 7, increases, until, when 7,0, the relation (2.5) is obeyed exactly | 
so that N,/Noq is unity at all temperatures. The physical explanation is that, 


since the number of electrons per atom in the s-band at absolute zero remains. |} 
constant at 0-6, an increase of €p, relative to <yq implies a corresponding decrease 


in the height v,(e9,), so that it becomes increasingly more difficult for a d-electron_ ]} 


to enter the s-band. For the smaller values of ry shown in the diagram the 
electron transfer may become quite considerable, For 79=10 and 7)=1-0, 
to consider an extreme example, N,/Noq is as much as 1-14, so that the number- 
of holes per atom at this temperature will have increased from 0:6 to nearly 0-7. 
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The possible applications of this effect to other physical properties have been 
discussed elsewhere. 

In figure 3 the three curves give the variation of A/Ay, with temperature for 
the values of 7p indicated. It will be noticed that A/Ay, is much less dependent 
on the value of r, than is N,/Noq, and, furthermore, that for any given value of 
T, A4/Ap increases with 7, reaching a maximum when 79> 00, Le. if the d-band 
is overlapped by an infinitely wide and low s-band. In that case the increase 
of the current constant is due entirely to a normal decrease, with temperature, 
of the work function as given by the relations (2.5) and (2.13). For finite values. 
of 7, however, the transfer of electrons from the d-band to the s-band causes 
a relative increase of the work function owing to the tendency of the Fermi limit 
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j Figure 3. ‘Temperature variation of A/Ap, 
Figure 2. Transfer effect for nickel, for nickel. 
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to shift to regions of lower energy (see figure 1). Hence the current constant 
will decrease relative to its particular value for rg>>0o. In the extreme case of 
7)=1 the increase in the work function due to the transfer effect just balances 
the normal temperature decrease, and A/A,=1 at all temperatures. This is 
at once apparent from the relation (3.1), which gives, with r9=1, na=7) ', so 
that, from (3.3), 4/A,=1. As has been stated in §2 an increase of the current 
constant will occur only if holes are relevant, while for particles A/Ay will decrease 
with temperature. For any finite value of 7, A/A, will tend asymptotically 
to a limiting value A,,,,./A9 as Ty 00, as is indicated by the discussion given in 
§2, where it was shown (equation (2.20)) that Amax/Ao = exP| | - . Calculations 
of y, using (2.19) can be readily carried out by employing a modified Taylor’s 
series method, the required derivatives of F(7) being given by McDougall and 
Stoner (1938). It is found that for 
Tay 2252), WAL | A, =9°5134 


and for 7220, y=—2:9844, Ansax/Ap=19-775. 
These figures indicate that A,,,./A) increases with 7, until, for 70, 
A max/Ag> ©. 


It is of interest in this connection to compare the results obtained by the 
methods of this section with those obtainable from the low temperature expres- 
sion (2.18). Selecting a particular value of 79, say 7) = 1-0, the three values of A/ Ap, 
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determined roughly from the curves of figure 3 are 2-35 (ry = 10), 2:56 (79 = 20) and 
2:79 (r5~ 0), while the corresponding values as calculated from (2.18) are ZO; 
2-18, and 2:28 respectively. The magnitude of A/Ag given by (2.21) (equivalent 
to equation (14) of Sun Nien 'T’ai and Band’s paper) is equal to that just given for 
ry >, ie. A/Ay=2:28. The discrepancy between the two values of A/Ay as 
calculated by the exact methods of this section and from the approximate relation 
is thus not very large, even forz7,=1-0. The considerations given in the preceding 
section show, however, that use of the approximation is nevertheless not justified, 
except when 7) <1. 


§4. COMPARISON WITH EXPERIMENT: CONCLUSION 


In order to compare the results deduced from the theoretical considerations 
of the preceding section with those that have been obtained by experiment, it is 
necessary to make an estimate of the bandwidth ¢9q for nickel. ‘Two methods of 
estimating this quantity from observational data may conveniently be employed. 
In the first place, if the low temperature specific heat of a metal can be expressed 
in the form 

CL. =yY T+ 6T?, 

the linear term gives the contribution of the electrons, and it has been shown 
(Stoner 1939 b), that for parabolic bands in ferromagnetic metals y may be 
estimated from the relative magnetization at absolute zero, ¢), and the band- 
width €9;. For nickel a consideration of the magnetic properties indicates that 
fy is very nearly equal to unity, so that from the experimental value of y, namely 
17-44 . 10-4 cals deg-? mol (Keesom and Clark 1935), €9;/R is found to be about 
2200° K., provided account is taken of the small contribution of the s-electrons. 
Alternatively <9;/R may be estimated from the observed Curie temperature, 
using the relation (Stoner 1938) 


€oq/R = O./f(o); 
applicable to parabolic bands, for which f(1:0)=0-41. Hence with 6,=631° k., 
€oq/R 1s found to be about 1600° k. The value of ¢9,/R may be estimated on the 
assumption that the s-electrons are “‘free’’, as in copper, giving «9,/R=6 . 104 °K., 
which has to be regarded as an upper limit. It is difficult to choose between 
the two estimates of €9q/k given above and also to see how far the premises under- 


lying the calculation of «),/R are correct. ‘The main difficulty appears to be the |} 
uncertainty as to deviations of actual band shapes from that postulated, and as | 
to the different effects such deviations may have on the various theoretically |] 


determined magnetic and energetic properties of metals. For the electron 


assembly considered both the d-band and the s-band may conceivably deviate | 
from a parabolic shape, and the extent of the deviations could only be gauged | 
if the calculations of Slater (1936) were to be carried further with an aim to | 


obtaining a higher accuracy in the final results. 


It is fairly certain, however, that the estimate of e)4/k given by Sun Nien T’ai_| 
and Band is far too low. In their paper it is stated that “in fact electronic specific _ 


heat data give €9q as small as 0-025 ev.’’, corresponding to €9q/k=300° K. The 
calculation outlined above, on the other hand, shows that the bandwidth deducible 


from the low temperature specific heat is 2200°k. Using their value of €oa/R, 7 


giving, for T= 1400° K., r7=4-7, Sun Nien T’ai and Band were able to explain 


the anomaly in the observed current constant for nickel quite well, However, | 


| 


| 
i 
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not only is the estimate of the bandwidth erroneous, but the low temperature 
expression for A/Ap will certainly not be applicable for ry~5. It is easy to see 
in fact that the agreement between theory and experiment is much more un- 
satisfactory. The most recent observed value of A is 1380 amp. cm? deg 
(Fox and Bowie 1933), the experimental results used for its calculation being 
obtained at a mean temperature of about 1400° kK. Assuming ¢,,/k =2000° x., 
i.e. tT, =0-7, figure 3 shows that A/A, will be of the order of 2:0, i.e. A240 
amp. cm” deg-®. ‘This discrepancy may be due to a number of causes, 
a summary of which has been given by Reimann (1934). Among these are: 

(1) Surface charging effects. Theoretical discussions of these effects have 
been given in a number of papers, references to which may be found in Reimann’s 
book. Since, however, the electrons emitted by hot bodies originate in the 
interior of the emitter, surface effects of this nature may be assumed to be of 
secondary importance only (see Reimann, p. 271). 

(11) Variation of emitter area. Since any metal surface will in general not 
be perfectly smooth the actual emitter area will be larger than that calculated 
geometrically. ‘This would imply an increase of the current constant above 
the value appropriate to the observed area. 

(111) Thermal expansion effects. These are too complex to discuss fully 
and reference should be made to the original papers. 

(iv) The main reason for the discrepancy appears almost certainly to be 
the imperfect purity of the emitter surface. As is well known even monatomic 
impurity layers may influence the thermionic emission very appreciably. This 
is borne out by the thermionic properties of platinum, which were first investigated 
by Du Bridge (1928) who found a value of A equal to 17000. Later work by 
Van Velzer (1933), who subjected his specimen to some very prolonged heat 
treatment, showed that for platinum A was very strongly dependent on the 
presence of gaseous impurities on the surface and in the surrounding atmosphere. 
The final value of A obtained by Van Velzer was equal to 170, and, it is stated,. 
“slightly lower values of A might have been obtained had the filament withstood 
further treatment”’. Fox and Bowie took great care to treat their nickel specimen 
thoroughly before the investigation. As Van Velzer had shown, however, 
a very stable state which resisted aging at 1650° k. was attained during the treat- 
ment of his platinum specimen. For this state A was abnormally high. Since 
all the heat treating processes for nickel were carried out at temperatures below 
1650° k., the anomalous result for this metal might also be partly indicative of 
such a state, corresponding to one of imperfect surface purity. 
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ABSTRACT. Anew method of investigating the discontinuous changes in magnetization 
(the Barkhausen effect) is described. The number of discontinuities occurring is counted, 
and a knowledge of the change in magnetic moment associated with each one, leads to 
information concerning the volume of material associated with a discontinuity. ‘The 
change in magnetization due to discontinuities and the distribution of number and size over 
the hysteresis cycle is obtained. Conclusions are drawn that the discontinuities associated 
with volumes of greater than 0:9 x 10-®°cm? only account for about 30% of the total change 
in magnetization. Possible explanations of these results are discussed, and suggestions made 
as to the further development of the method. 


SINR OD WC rOink 


T is well known that on the steep part of the hysteresis curve for ferromagnetics, 
the change in magnetization proceeds discontinuously (the Barkhausen effect) 
and that the discontinuous changes in flux can be amplified as voltage pulses, 

and displayed on an oscilloscope or made to operate a loud-speaker, when the 
familiar ‘‘rustling”’ sound can be heard. 

The discontinuities have been identified with changes in the directions of 
magnetization of the domains of ferromagnetic theory. The volume of material 
constituting a domain and the extent to which the discontinuous processes account 
for the change in magnetization, are matters of considerable interest, and have been 
investigated by various observers. 

Barkhausen (1919) in his original paper stated that the larger discontinuities 
corresponded to the reversal of the magnetic moment in volumes which can be 
estimated to be between 10~° and 10-§ cm’. ‘Tyndall (1924) also found that the 
larger discontinuities corresponded to a volume of 10-® cm®, and that the period of 
the discontinuity was about 10-* seconds. 

Bozorth and Dillinger (1929, 1930, 1932) in a series of comprehensive experi- 
ments, examined the size of discontinuities in wire specimens, together with their 
distribution over the hysteresis cycle. They came to the following conclusions : 
(a) The average volume of discontinuities varies over different parts of the field 
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range, the maximum occurring in the steep part of the hysteresis curve. (6) For 
annealed iron, hard-drawn iron, hard-drawn permalloy and annealed nickel, the 
discontinuities corresponding to volumes greater than 10! atoms (1.0; * em); 
accounted for most of the change in magnetization over the steep part of the 
hysteresis curve. (c) The maximum average volume ranged from 1:2 x 16-® cm? 
for annealed iron to 4-5 x 10-®cm3 for malloy. (d) It was suggested that these 
volumes could be identified with the ferromagnetic domains of magnetization. 

Bozorth and Dillinger’s method was to amplify the voltage pulses produced in 
a search coil by the discontinuities. From the output current of the amplifier, 
dB/dt, the rate of change in induction in the specimen caused by discontinuities 
could be calculated. 

Knowing dH/dt, the rate of change of magnetizing field, dB/dH could be 
found. dB/dH, given by this method of integrating the separate pulses, was 
compared with the value of dB/dH found from the hysteresis curve. 

The average volume of material 3, the magnetization of which changes from 
saturation in one direction to saturation in the other, is obtained from an expression 
giving v°/d in terms involving, among other factors, the mean square current of 
the amplifier, the reversible permeability and the conductivity of the specimen. 

With the development of modern methods of counting electrical pulses, it was 
thought that it would be possible to count the number of discontinuities of various 
sizes which occur when a ferromagnetic is slowly magnetized, so that distribution 
curves of number against size could be obtained. Further curves could be derived 
showing the contribution to the magnetization of discontinuities of various sizes, 
and also the total change in magnetization caused by the total discontinuities 
detected. This appeared to be a very direct method of approach to the problem. 

In this preliminary paper, a description is given of the apparatus developed for 
this purpose, together with some typical results. The work shows that the 
method is entirely practicable, and interesting results have already been obtained. 
The character of the results indicates the directions in which more detailed 
measurements are required and the developments of the method which are desir- 
able. The results of more comprehensive investigations will, it is hoped, be 
presented in due course. 


92. AP PARA WIS 


A block diagram of the apparatus is shown in figure 1. ‘The two search coils 
are connected in series opposition and the specimen is placed inside one, this 
assembly being enclosed in the magnetizing coil. By altering the current in the 
magnetizing coil by means of a current changing device, the magnetizing 


field can be slowly changed at an approximately uniform rate. 


The changes in flux caused by Barkhausen discontinuities are picked up by the 
appropriate search coil, and amplified as voltage pulses by the amplifier. The 
discriminator selects pulses of amplitude greater than any desired value, and these 


| are then counted by a scale-of-sixteen unit and a mechanical recording counter. 


The characteristics of the apparatus are as follows. ‘The magnetic field 
changing unit enables the current in the magnetizing coil to be slowly changed. 
This is done by charging and discharging a condenser, arranged to alter the grid 
bias of a pentode valve, thus causing the anode current to change at the required 
rate. The unit is based ona circuit described by Bozorth (1929). ‘The field can 
be changed at a rate of 2 to 20 milli-oersteds per second for a maximum field of 
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20 oersted. For smaller maximum fields the rate of change is proportionately 
less; in the case of hard-drawn mumetal it was 0-15 to 1-5 milli-oersted per second 
for a maximum field of 2:6 oersted. 


ARTIFICIAL 
PULSE 
UNIT 


RECORDING 
COUNTER 


DISCRIMINATOR. 

SCALE—OF— 
SIXTEEN 

UNIT 


SEARCH COIL 
MAGNETISING 
COILS, 

SPECIMEN. 


AMPLIFIER 


MONITOR 
C.R.O. 


Figure 1. Block diagram of the apparatus. 


The amplifier is of conventional design, consisting of five stages of resistance 
capacity coupling. The noise level with the search coils connected was about 
2-5 pv., the first valve being carefully chosen to generate as little noise as possible. 
The frequency response is shown in figure 2. 

The discriminator is based on a design due to Lewis (1942) and selects pulses 
above any desired size, enabling the various magnitudes of change in magnetic 
moment, corresponding to Barkhausen discontinuities, to be measured. The 
discriminator would differentiate between pulses whose amplitudes differed by 
only 0-5 volts. ‘The scale-of-sixteen unit is built up from the hard valve scale of 
two circuits originally described by Stevenson and Getting (1937). This scaling 
unit enabled the actual counting to be done by a recording counter of the watch 
type (Neher 1939), which had a minimum resolving time of about 1/300 second. 
A special driving unit was used to operate the recording counter. 

An artificial pulse unit was constructed which gave square pulses of about the | 
same time duration as those produced by 
Barkhausen discontinuities; this was used to 
calibrate the amplifier and to test the counting 
circuits. It is of conventional design; the pulse 
recurrence frequency is determined by a resis- 
tance capacity phase shift oscillator, the sine- 
wave output being fed into a squaring valve. 
The square wave so produced is passed into a 
circuit of short time constant for differentiation. ~o i0 io" a 
The negative and positive pulses produced are PO ara 
passed into a third valve which removes the 
positive pulse by grid current action, and ampli- 
fies and squares the negative one to give a positive pulse output. Varia- 
tion of the short time constant for differentiation provides a control of pulse 
duration. ‘The low impedance output is provided by a cathode follower. The 
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pulse recurrence frequency is variable from 100 to 200 cycles per second and the 
pulse duration has a range of 300 to 1000 microseconds. 

The apparatus can detect changes in magnetic moment of greater than 3 x 10-8 
E.M.U. Thus if these changes correspond to complete reversals of the direction of 
magnetization of domains, the apparatus will detect domain volumes greater than 
0-9 x 10° cm? in the case of iron where the saturation intensity J, is about 1-7 x 10°. 
The domain volume is given by the equation: 


215=AMIV, 
where AM is the change in magnetic moment in a volume V. Pulses can be 
counted at a maximum speed of 5000 per second. 

The essential limitation to the sensitivity of the apparatus is set by the noise 
generated in the amplifier, a large part of which is caused by thermal agitation noise 
set up by the resistance of the search coils. Hence in the design of the search coils 
it is important to keep the resistance as low as possible, consistent with the necessity 
of obtaining the maximum sensitivity.. Each search coil consists of 50000 turns 
of 44.s.w.G. copper wire, wound in the form of a multi-layer solenoid, 10 cm. long, 
0-68 cm. internal diameter and 1-63 cm. external diameter. With these dimensions 


MAGNETISING COILS 


SEARCH COILS 


SPECIMEN 


ALL DIMENSIONS IN CM$ 


Figure 3. Layout of specimen, search coils, and magnetizing coil. 


the response to the changes in flux over most of the length is sensibly constant. 
Two search coils in series opposition are used to reduce the effect of external 
interference. Figure 3 shows a section through the assembly of search coils, 
magnetizing coil, and specimen. The assembly is enclosed in a mumetal box 
which also reduces the effects of external interference and the effect of stray magnetic 
fields. With the specimen in the form of a wire 7-5 cm. long placed in the position 
indicated in figure 3 only one search coil is affected by Barkhausen discontinuities 
taking place in the specimen. 

A fuller description of the apparatus is given in a thesis by one of us (H.D.B. 
294-7). 


§3. CALIBRATION OF THE APPARATUS 


The time duration of a Barkhausen reversal, as seen in the output of the 
amplifier, is of the order of 300 microseconds, and depends on the time constant 
of the search coil and amplifier, and on the rate of decay of eddy currents in the 
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specimen. However, the actual change in magnetic moment, and resultant 
initial change in flux, take place in a time which is short compared with the time 
constant of the apparatus. Under these last conditions, the amplitude of the 
resultant output pulse, in the case of a linear amplifier, should be proportional to 
the change in flux in the search coil. 

A calibrating coil, 5cm. diameter, 10cm. long and two turns per cm., was 
placed over the search coil which normally contained the specimen. ‘The flux 
change produced in the search coil by a known current passing though the cali- 
brating coil could be calculated, making allowance for the variation in field along 
the length of the coil and over the cross section of the search coil. The relation 
between flux change in the search coil, and the change in magnetic moment 
resolved in a direction parallel to the applied field, for a discontinuity taking place 
in the specimen, is given by: 

Flux change = 27 AMN/(r? + 1?)}, 
where AM =change in magnetic moment; 
N=number of turns in search coil; 
r= effective radius of search coil; 
21=length of search coil. 


This formula applies strictly only for a discontinuity taking place at the centre 
of the search coil. However, as the ratio of length to diameter of the search coil is 
9:1, the mean variation in sensitivity over the length occupied by the specimen 
(3/4 that of the search coil) is only 1%. 

The change in flux was produced by passing a current pulse of square wave 
form (from the artificial pulse unit), through the calibrating coil, with various high 
resistances in series. In this way it was possible to detect the variation in gain of 
the amplifier caused by non-linearity, and any effect of leakage across the series 
resistance. 

The characteristics of the search coil were such that the rapid change of flux 
tended to set up an oscillation in the search coil, ata frequency determined by its 
self-inductance and the associated circuit capacities. There was sufficient 
damping across the coil to prevent more than one cycle of oscillation. Experiment 
showed that the amplitude of the voltage pulse so produced was proportional to the 
change in flux, but independent of its duration over a time range of three to one | 
(300 to 1000 usec.). Pulses from the amplifier were fed into the discriminator, 
and the setting of the discriminator required to cut off pulses of a given size | 
determined. In this way it was possible to obtain a calibration curve showing 
variation of change in magnetic moment with the corresponding discriminator 
setting. 

The counting system was tested for accuracy and speed of counting by means | 


of the pulse unit; with sinusoidal pulses, counting speeds up to 5000 cycles per | 
second are possible. 


§4. EXPERIMENTAL METHOD 


The discriminator was set at a given level; the specimen was then magnetized 
to saturation in one direction and the field gradually reduced, reversed and then | 
increased to saturation in the opposite direction. Readings of the recording 
counter were taken, in the case of iron at intervals of 2 oersteds, and in the case of | 
other specimens at other suitable intervals. This was repeated a number of times 
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at each level of the discriminator and an average taken in each case. A curve was 


| thus obtained for each specimen giving the number of discontinuities N above a 


given minimum size AM plotted against size, as in figure 8 for annealed iron. By 
differentiating this curve it was possible to obtain a size/number (AM vs dN/dAM) 
distribution curve and hence a curve giving size, AM, plotted against contribution, 
AM.dN/dAM;; see figures 4, 6, 8 and 5, 7,9 respectively. 

The area under the contribution curve 


dN 
Jam. Fam: 4M), 


with appropriate limits, is the total change in magnetic moment contributed by 


| discontinuities which can be detected above the noise level of the amplifier. This 
) total divided by the volume of the specimen gives the change in intensity, contri- 
| buted by the Barkhausen discontinuities measured, over the part of the hysteresis 


cycle considered. 
When the discriminator is set near the noise level there is a background count, 


} constant for each setting, which is appreciable and increases as the noise level in 


approached. Ifthis becomes comparable with the true count; the error in counting 
becomes large. It was found that accurate results could not be obtained with a 


| background count of greater than 10 per second at the recording counter. A lower 


limit to the discriminator setting was therefore put at AM=3-0x 10-*E.M.U.; 


| above the setting AM =3-75 x 10-*r.m.u. the background count was negligible. 


SS RE SOAS 


The following materials, in the form of thin wires 7-5 cm. long, were 


} examined :— 
Diam. cm. Composition 
Annealed iron 0:0376 99-92%, Ke, 0:03% C, 0-05°, Mn. 
Hard-drawn iron 0-0376 
Hard-drawn nickel 0-0458 D952 Ni 
Hard-drawn mumetal 0-0254 7 Op Nt e629 Cul, lean ai ve, 
06% Mn. 


The hysteresis curves for these materials are shown in figures 4 to 7. ‘These 


| curves were obtained using the standard ballistic method, with the specimen in 
| the form of a long wire, and a magnetizing coil 40cm. long and 1-5 cm. diameter. 


The change in flux in a search coil of 1000 turns wound round the centre of the 


|| specimen, was measured by a Grassot fluxmeter (1:43 x 10° maxwell turns per 
| radian). 


The distribution curves are shown in figures 9, 11,13. The number of counts 
per unit range of AM plotted as ordinates are reduced to number per cm* in order 


} to compare specimens of different sizes. In all cases the curves rise steeply near 
| the noise level. - 


Size contribution curves are shown in figures 10, 12, and 14, giving in the case 


‘.of annealed iron a reduction in gradient near the noise level. ‘This reduction 


may be incorrect because of experimental errors which occur near the noise level, 


» and the errors which are inevitable in differentiation. 


25-2 


Figure 4. Hysteresis curve for annealed Figure 5. Hysteresis curve for hard-drawn 
iron. iron. 
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Figure 6. Hysteresis curve for hard-drawn ‘Figure 7. Hysteresis curve for hard-. 
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Figure 9. Distribution curves for annealed iron, ©, and hard-drawn iron, A. 
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Figure 10. Contribution curves for annealed iron, O, and hard-drawn iron, A, 
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The contributions to the total change in magnetization which are from the areas 
under the curves 10, 12, and 14, are given in table 1. The contribution of dis- 
continuities lying above the noise level and below 20 x 10-%E.M.u. are shown 
together with the extrapolated contributions. The extrapolated contributions 
were found in the following manner. The number-level curve over the range 
which could be examined is clearly of the exponential type. A reasonably satis- 
_ factory fit is obtained with a curve of the form Ae~*", Differentiation of this gives 


n ap ee ms “ Lose * 
the expression 4Aax~e att and the contribution takes the form 1Axxe~%* 
p 2 2 ; 


378 H. D. Bush and R. S. Tebble 


Table 1. Change in magnetization contributed by Barkhausen discontinuities. 


Contribution of Barkhausen discontinuities 


Extra- 
Total From 3 to polated dH 
Material ify fy change 20X10-° % of total contri- % of total dt 
in I E.M.U. bution 10-3 oersted/sec- 
0 to oo 
Annealed 
iron 1700 140 1900 605 5 1300 68 3-15 
Hard-drawn 
iron 1700 102 1540 475 Sil 840 55 3-15 
Hard-drawn ; 
nickel 500 12-2 520 56 11 120 24 6-30 
Hard-drawn 
mumetal 750 430 1040 320 31 520 50 0-15-1°5 


The area under the contribution curve from AM=0 to AM= o is given by 


pe 3 
[gansta = 2A lor, 
0 


which was taken as the extrapolated value of the total change in magnetization 
contributed by Barkhausen discontinuities. 

As can be seen in the case of all the materials examined, the total change in 
magnetization is not accounted for by the Barkhausen discontinuities detected. If 
this comparison is restricted to the steep portion of the hysteresis cycle the agree- 
ment is no better. Further, the rate of change of the magnetizing field was varied 
within the limits shown, and this had no appreciable effect on the number of 


counts obtained. Some of the results for nickel, indicating this, are shown in 
table 2. 


- 


Table 2. Effect of variation of rate of change of field on number of counts 
obtained for hard-drawn nickel wire 


Discriminator setting (4M) : 
{Oxe RMA 3°74 4-14 6-9 
Time taken for one 
half-cycle I PA PA SP. BY) 23, 23, 25, 34, 40 27, 27, 30, 60 
(minutes) 
No. of counts in one 
Hee 106, 119, 109, 103, 110 89, 90, 88, 88, 95 39, 40, 44 40: 
x 10- 
Mean (x 107%) 110+3 90+2 3912 


With the results obtained it was possible to show the variation in the number of 
discontinuities per unit range of magnetizing field dN’ /dH, where N’ =dN /dAM 
with field H, for different sizes of discontinuities. In order to display these 
results on a convenient scale, the contribution AM.dN’/dH has been plotted 
against the field. ‘These curves show maxima near the point H=H, where H 
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is the coercive field, as might be expected, and the 
contribution falls in the region where reversible changes 
predominate (figures 14, 15 and 16). 

It should be possible to plot the number of dis- 
continuities per unit change in intensity, dN’ /d/, against 
intensity J for various sizes of discontinuity and thus 
obtain some information on the distribution of the 
discontinuities (number, size and contribution) as the Figure 18. Contribution per 
intensity is altered over the hysteresiscycle. Unfortun- unit change in H—annealed 
ately it is difficult to obtain a sufficient number of 9 °°" 
readings where the intensity is changing rapidly © AM= 32x10-6eM.v. 
(near HEIEOY ae Ot eee 

The “average size” of domains detected can be “ “<0 PY 

ge size’? of domains detected can be 

obtained by dividing the total number of discontinuities counted into their 
contributions for different field ranges. The variation of this “‘ average size” 
over the hysteresis cycle is shown for annealed iron in figure 18; this is a typical 
curve, showing a maximum near H=H,. ‘This curve is similar to those obtained 
by Bozorth and Dillinger. The significance of this curve is limited, as it only 
applies to those domains which are detected. If all the domains were accounted 
for the curve might be considerably different. 


§6. DISCUSSION 


For the materials examined it appears that discontinuous changes in magneti- 
zation in domains larger than 0-9 to 3-5 x 10-® cm according to the material, do 
not account for the total change in intensity over the steep portion of the hysteresis 
curve. It would seem that at least some of the discrepancy is caused by domains 
of a size smaller than those detected by the present apparatus. ‘The extrapolation 
of the experimental curves gives some estimate of the effect of these smaller 
discontinuities. Another factor which may account for some of the discrepancy 
is the extent to which reversible changes in magnetization take place over the 
hysteresis cycle. A quantitative estimation of this effect can be made by measure- 
ment of the reversible permeability at various points on the hysteresis cycle. 

There are two ways of increasing the 
sensitivity of the apparatus so that smaller 
discontinuities can be detected: (a) the [ 
noise level of the amplifier could be [ 
reduced; (b) the sensitivity of the search 2 a 
coil could be increased. Is 10+ 

In the case of (a) it is doubtful i 
whether the valve noise can be further 
reduced to any appreciable extent. It L 
does seem possible, however, that the [ 
search coils might be re-designed to give 
a greater sensitivity, without at the same Up pS eS, 
time increasing the thermal agitation 2 
noise. ‘This presents difficulties, as in- Figure 16. Contribution per unit change in 
creasing the physical dimensions of the H—hard-drawn iron. 
search coil will increase the effect of in- © AM= 32x10-6B.mu. 
terference from external sources. An [] AM= 60x 10-S Emu. 
additional method might be to reduce A, AM=12°3x 10-6 amy, 


xoe [7 T T | aa iy] le eles a1 at 2) ay tae T 


; a 
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the temperature of the search coils, since the thermal agitation noise decreases 
with both resistance and temperature. 

The effect of the rate of change of field could be investigated over a wider range, 
although an increase in this rate would appear to increase the possibility of coinci- 
dences in the discontinuities. "The reduction in size of the specimen would reduce 
the effective rate of change of field by reducing the number of domains counted. 

It is desirable that more results should be obtained near the region H=/,, 
possibly by increasing the demagnetizing factor of the specimen, so as to reduce the 
effective gradient of the hysteresis curve in that region. Annealed mumetal would 
appear to be suitable for this purpose. It would be further interesting to investi- 
gate materials of the mumetal type, as it is possible that domain size bears some 
relation to initial permeability. 

It is hoped to carry out experiments on materials covering a wider range of 
magnetic properties than has been possible so far. 
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Figure 17. Contribution per unit change Figure 18. Average domain size for 
in H—hard-drawn nickel. annealed iron. 


© AM= 3:9x10-6 EMU. 
[-] AM= 62x10-6EmMu. 
2, AM=13:2x10-6 pm.v. 
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DISCUSSION 

Prof. Batrs. I should like to emphasize the importance of this work, which I very 
much hope will be extended. The trouble is that the Barkhausen effect is comparatively 
easy to show on the lecture table, and people have therefore been inclined to think that 
everything is known about it. Therefore it is only when a serious and careful study is 
made that we find, as the authors have clearly shown, things are not at all satisfactory. 
Moreover, there has been a tendency, particularly in the American and German literature, 
to give magnetic domains a kind of rigid cubic-block structure that they cannot possibly 
possess. I very much hope that the extension of this work will enable us to get a clearer 
view of domain processes, and, in particular, help us to decide on the correctness of the 
views concerning fillet-shaped domains which have recently been put forward. 


Prof. W. SuckxsmirH. Do the authors believe that the discrepancy between the 
integrated sum of the intensity changes as determined by the Barkhausen discontinuities 
and the total change of intensity determined in the usual way is due to irreversible groups 
below the noise level, or to continuous and reversible changes ? 


Mr. McCaic. I should like to point out that fast counting rates can now be measured 
without any special form of mechanical counter. Counting speeds of more than 100,000 
per second have now been achieved in commerical scaling circuits. 


AuTHors’ reply. Measurements on reversible permeability are being carried out at 
Leeds, the results of which suggest that reversible changes account for about 10°% of the 
total change in magnetization over the range considered. It is possible that some of the 
discrepancy may be due to the assumption that the time constant of the discontinuity in 
the specimen is short compared with that of the search coil. ‘This is based on a formula 
of Wwedensky (1921). It is intended to carry out measurements on the time constant 
of the discontinuities, making use of a wide-band amplifier (1 kc/s. to 1 Mc/s.) with which 
it should also be possible to detect the smaller discontinuities. 


WwEDENSKY, 1921, Ann. d. Phys., Lpz., 64, 609. 


A Radiosonde Method for Atmospheric Potential 
Gradient Measurements 


By R. E. BELIN 
Dominion Physical Laboratory, Wellington, New Zealand 


MS. received 9 Fuly 1947 


ABSTRACT. Amethod is outlined whereby the National Bureau of Standards radiosonde, 
in a modified form, is applied to the measurement of the earth’s potential gradient. This 
paper describes the modification of the sonde circuit to measure the point-discharge current 
‘produced by a gradient and the result of a typical sounding of a cumulo-nimbus cloud. 


§1. INTRODUCTION 


XPERIMENTS on the variation in natural potential gradient have so far 
been almost wholly carried out on or near the surface of the earth. ‘The 
main exceptions were in the work of Simpson and Scrase (1927, 1935), 
who used a balloon to carry aloft, to heights of the order of 10 km., recording 
equipment capable of measuring pressure, humidity and potential gradient. 
‘The pressure and humidity were recorded, in the usual manner, by employing 
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an aneroid barometer and a human hair hygrometer respectively; the measure— 
ment and recording of the gradient was, however, more complex. Simpson 
and Scrase relied upon the fact that, if two sharp points with opposed senses. 
are connected by a straight conductor, the axis of which is parallel to a field 
of suitable magnitude, point discharge will take place, thereby tending to equalize 
the potentials between the respective parallel planes through the points. In 
their apparatus this discharge current was caused to flow between two electrodes. 
on polarity paper, and the magnitude of the gradient could be roughly estimated 
by measuring the width of the line traced by the positive electrode. Using 
a 20-metre length of fine copper wire to separate the points, it was possible to- 
measure fields not smaller than 10v/cm. and the main disadvantage was that 
results depended upon the recovery of the equipment intact. 

The method employed by the author (1946) differs essentially from that of | 
Simpson and Scrase in that a permanent record is obtained at the receiving 
station so that results do not rely on the recovery of the sounding equipment. 
The point-discharge method was used in conjunction with a radiosonde instead 
of polarity paper. This was so modified that it was capable of transmitting 
a carrier, the modulation of which was dependent upon the magnitude of the 
electric field. These signals were detected and recorded by the standard radio- 
sonde receiver developed by the Bureau of Standards (Diamond et al. 1940). 


§2. METHOD 
If two points be placed in an electric field as explained above, and a large- 
valued resistor be connected in series with the collector wire, at its centre, 


then the point-discharge current z flowing out of the points will be given (Simpson. 
and Scrase 1927, 1935) by the relation . 
t=a(F?— M7) Eee (¢) 

F is the existing potential gradient (greater than VM). Mis the minimum gradient 
required for point discharge to commence. a is a constant depending mainly 
on the exposure of the points (i.e. pointedness, number of points, distance between 
them etc.) and to a lesser degree upon the value of the series resistor. This. 
current z produces across the resistor a potential drop which is measured and. 
transmitted by the modified radiosonde. 

Thus by carrying aloft (with a neoprene balloon) the collector points separated’ 
by a length L of wire with the modified transmitter at its centre the measurement 
of gradient can be readily achieved. 


§3. MODIFIED RADIOSONDE 

The radiosonde developed by Diamond, Hinman, Dunmore and Lapham: 
(1940) consisted of two oscillators operating at frequencies of 1 Mc/s. and 70 Mc/s.. 
respectively. The 1-Mc/s. oscillator M (figure 1) was of the squegger type.. 
Its squegging frequency was determined by the R,C, combination introduced’ 
into the grid circuit. ‘The oscillator C was also of the squegger type being on: 
when M is off and vice versa, thus effecting a modulation. 

When M is oscillating it draws grid current which charges up the condenser 
C, and thus developing across R, a potential which builds up and eventually 
blocks the tube. The period of no oscillation is determined by the time constant 
R,C,. If a bias battery P be introduced, as is done in figure 1, the squegging: 
frequency will be altered. It is found that the squegging frequency has. an: 
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approximately linear relationship with the potential developed across R, (figure 2). 
This linear relation holds for positive bias, but only approximately for negative 
bias up to cut-off. Figure 2 gives steady grid voltage as a function of squegging 
frequency, when R,=2 MQ and C, =0-003 pF. 


Low reference conts. 


(ATT 


Commutator 


Aneroid 


Figure 1. The modified radiosonde circuit. 
R=1 MQ, R,=2 MQ, R3=1:25 MQ, C.=:005 ur., Ch=:003 pr. 


In figure 1 the small switch S, actuated by 
the standard relay S’ on the 5, 10, 20, 25 etc. ue re 
commutator contacts, was used to introduce 


periodically into the grid circuit of M a fixed a8: 
time constant R,’ in parallel with C,’. This > 
periodicity was determined by the rate of 2% 
change of pressure, or height, and hence was Sl 
only controlled by the baro-switch. (The S 
baro-switch consisted of a commutator com- & 


posed of 80 separately insulated silver contacts Wa 
of which only 5, 10, 15, 20 etc. were used. i x 
Over this commutator the contact arm of the 

aneroid barometer could move freely so that 
its position was determined by the pressure, 
that is by the height.) By means of this 
reference signal, as it was termed, it was ~ 0 voureeine 
possible to evaluate the altitude of the Figure 2. 
sounding equipment, as well as to obtain 

a check on any possible oscillator drift. The values of R,’ and C,’ were chosen 
to produce a squegging frequency of 190 c/s. and were used as the high reference 
value. 

The baro-switch also introduces on the 15, 30, 45 etc. contacts, a resistor R 
in parallel with R,. This latter combination produced, when no point discharge 
was taking place, a low reference of 140 c/s. Thus this low reference was subject 
to a change, the magnitude and direction of which would be dependent upon 
the magnitude and sign respectively of the potential gradient. However, the 
chief function of this reference was not so much a check on the oscillator as a 
check on the vertical course of the balloon. It served also as a height signal. 
It is possible should the ascent take place through a cumulo-nimbus cloud, as. 
was intended, to encounter such a strong downward current as to produce a 
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descent of the balloon. This would produce a change in the predetermined 
sequence of high and low references which would be observed on the record. 

The value of R, and C, were chosen to give an audio frequency of 50 Cis 
with no bias battery in circuit. When the small 1}-volt penlight cell was 
introduced in series with R,, a squegging frequency of 64 c/s. was obtained. 
Thus the cut-off was moved 14 c/s. further to the left. ‘This was an advantage 
should a large negative gradient occur. 


§4. THE COLLECTORS AND RADIATING AERIAL 


Figure 3 shows a typical collecting device. It consists of pins soldered to 
the ends of the wires which form the framework. Two of these collectors with 
pins in opposed senses were used at either end of a gauge 30 copper connecting 
wire, at the centre of which was located the modified radiosonde. In order to 


Figure 3. A typical collector, showing the manner in which the sharp points were connected. 


support the sounding equipment which includes the collectors and radiosonde, 
the thin wire was wound round some strong hemp string, while the upper collector 
was fastened with a similar string, 10 feet long, to the balloon. 

‘The transmitter aerial, a half-wave dipole, was fastened parallel to the collector 
wires, but separated from them by stand-off insulators of cork (for lightness) 
impregnated with polystyrene solution in toluene. 

The length of the collector wire was chosen according to the formula (2 — 1) 
half-wavelengths (where n is an integer) to give, as it was experimentally shown, 
a maximum signal strength. ‘The value of 7 was chosen as 5 which thus gave 
a length of 64 feet between the collectors for a frequency of 69 Mc/s. The 
64-foot connecting wire consisted of gauge 30 enamelled copper wire. For this 
length the minimum potential gradient capable of beginning point discharge 
was determined experimentally as 3v/cm. at N.T.P. Figure 1 shows schematically 
the collector wire connected to either end of the R,C, combination, while the 
actual radiating antenna is seen on the right of the circuit. | 
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Son LH ES CARIB RATT ON 


Calibration was carried out by means of an adjustable electrostatic field 
between two large parallel metal sheets. In this field were placed the collectors, 
actually used, with the pins parallel to the field direction. The collectors were 
connected by a short wire with the resistor R, in series, as the 64-foot length was 
impracticable. 

The value of the potential between 
two parallel planes passing through the 
points was determined from the relation 


potential = Ldv/dx volts, 


where L is the distance in cm. between 
the pin points and dv/d~ the field strength 
in volts/em. For points of given shape 
the value of Ldv/dx was shown to be 
constant for any value of L. Hence by 
calibrating this for a small given value 
of L, the calibration of frequency in 
cycles/sec. as a function of Ldv/dx was “9 2 4% 
the same for any value of L*. A cali- 
bration of a typical sonde using 7 points 
on both collectors is shown in figure 4. 
It will be observed that the minimum 
value of Ldu/dx necessary to initiate point 
discharge is approximately 6800 volts, 
corresponding, for the length of collector 
wire used, to 3-5 v/cm. ; it is thus equal Figure 4. 

to the value of M required in equation (1), 

at the existing pressure, temperature and humidity. 

The minimum value of Ldv/dx necessary to produce a discharge was shown 
to be independent of the number of points, but once discharge had set in, the 
point-discharge current became directly proportional to the number of points. 
Thus when the current for a given number of points is measured, that for any 
other number of points can be readily determined. 

Point-discharge current is also a function of pressure related by the formula 
(Tamm 1901) 


LBY IN THOUSANDS OF VOLTS 
WO 


ECOusLIOO mal 14 I 
DEFLECTION IN C/S. 


Pina = Va t02\ lows iP) 2 8 8 eee (2) 
wherez and 1, are the currents at pressures P and P, respectively and V = Ldv/dx. 
P, is the pressure at ground level at time of releasing sounding equipment. 
This cumbersome formula may be reduced, without affecting the accuracy 
required, to the approximate expression, 
Deine ae ay pee) See (3) 
which was used in these experiments. 
The height of the sonde was determined from the formula 
=OlOo a Pee ER Oe et (4) 
combining equations (3) and (4) we obtain 


Cet ee ed 
Pea eag 


* By increasing the value of Lit is possible to make the apparatus as sensitive as is desired. 
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Hence expression (1) can be modified to include a correction for height 
i=a(F?—2)eo"%, ee (5) 
By the use of this formula the actual value of F, i.e. du/dx, can be determined 
at any height. 
§6. RESULTS 

The equipment outlined above was used with only one modification, viz. 
thirteen points were used on the collector frames instead of seven. ‘This altered 
the calibration by the factor 13/7. With the temperature, pressure and humidity 
prevailing at the time of the following ascent, the striking gradient value was 
observed to be 3 v/cm.; this is the value used in figures 5 and 6. 
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Figure 5. 


fep)| 


Figure 6 shows a typical record obtained by releasing this apparatus beneath 
a cumulo-nimbus cloud with fracto-stratus formation below. ‘The field at ground 
level was 4 v/cm. positive, that is, downwards. It increased in a fluctuating 
manner to a maximum of about 5 v/cm. at 600 feet, and then slowly decreased 
to zero at 2000 feet. ‘The graph has been smoothed particularly in the lowest 
negative bulge between 2000 and 3500 feet which region probably corresponds to 
the fracto-stratus region. On entering the cloud at 3500 feet the potential gradient 
increased rapidly to a high positive value (off record) in the vicinity of 7000 feet. 
It then began to decrease again and at 10000 feet suddenly reversed to a high 
negative value (off record). At this stage the gradient was so large that the 
discharge current produced a voltage across the grid resistor R, which was well 
beyond the cut-off value. However, a signal was still detected in the form of the 
69-Mc/s. carrier wave. At 14000 feet, the cloud top, the gradient slowly became _ | 
less negative, swinging over to positive when leaving the cloud, and finally 
reaching another subsidiary maximum of about 6 v/cm. at 14600 feet. Above 
this the gradient decreased again according to the inverse square law until at 
about 20000 feet the point discharge current became zero. 

A reproduction of a part of the actual record obtained, in this case, when 
emerging from the cloud top is shown in figure 5, which corresponds to the upper 
partoffigure6. ‘This record shows, by the scattered points, the fluctuating manne 
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in which the gradient increases or decreases during an ascent. This fluctuating 

increase could be explained by slight variations in the state of ionization with 

height, or possibly a humidity variation which might produce the same result. 
Interpreting the curve of figure 6, we note that the cumulo-nimbus appeared 


to have a large positive distribution of electricity in the base and a large negative 
distribution in the top. This 


corresponds to a strong posi- 20 
tive gradient below and above Pie a leat 
the cloud, the latter of which Bd ya | 
was observed. ‘The former t 
was distorted by the presence i ea 
of the fracto-stratus formation - | 
on which it was assumed the _— 
main field of the cloud termin- eerGun TOP | 5 (ae So la (I 
ated, thus producing in effect a c Pal asa | | ial 
negative anda positive change >| | | pet 
in the upper and lower parts a t Tea | 
of the formation. 2 i 
<10 
§7. CONCLUSION 3 | 
The method outlined above - | | | [ 
isapplicabletomany problems ~~ : t 
ee to rea eeaiene 5 gees aaa eae ima oe eae 
it is excellent, as proved, for Er [ hk KE 
the measurement of electrical a 
distribution of clouds; it is sey a PS - 
possible, by lengthening L, to T | (nes 
measure the fluctuation in ical 
dv/dx with height on a fine += ieee 
day, when it is approximately Penal [ 
equal to lv/cm.; by slight <4 +3 4 5 6 
VOLTS PER CM 


modification to the sonde it 
is possible to include in the 
record both humidity and temperature and thus obtain reliable quantitative data 
on the true variation of the potential gradient; again, by a further slight modifica- 
tion to the circuit which is at present being attempted, it will be possible to obtain 
measurements of very large gradients, especially negative, which have hitherto 
been beyond the range of the apparatus. 


Figure 6. 
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A Note on the Lecher Wire Method of Measuring Impedance 


By MARJORIE WILLIAMSON 
Bedford College, London 


MS. received 18 August 1947 


ABSTRACT. Some time ago a new method of measuring impedance by means of a 
Lecher wire system was described by Dr. Williams in these Proceedings (1944). This 
method would seem to be very useful for the investigation of the properties of dielectrics 
since it gives both the real and the imaginary part of the dielectric constant, and it is also 
adaptable for use over a considerable range of frequencies. 

In the course of some measurements on these lines at a wavelength of about a metre a 
number of interesting points have arisen. 


across one end of the Lecher wires which are loosely coupled to an 

oscillator at the other end. ‘The standing-wave pattern is investigated 
by means of a pair of movable bridges consisting of vacuum thermojunctions 
of impedance Z, and Z, respectively, the values of which need not, however, 
be known. 

In the position shown in 
figure 1 the currents passing 
through the bridges are J, > 
and J, respectively. It has 
been shown (Williams 1944, 
equation (10)) that the 
modulus p of the ratio J,/J, is given by 


K,+ Kg sin 2(b+ Bs) 


Bie in this method, Z, the impedance to be measured, is connected 


p?=Ky+ ‘sinha fein? (bs) en nn (1) 

where f =27/A, the current reflection coefficient 
(24) =2)| (Zot 2) =e 7, (2) 
K,=(1/| Z,|*)[((R3+X2) cos 2Bs, + X_Z, sin 28s, 4 Z sin2Bs,],  -.. . (3) 
K,=(1/| Z,|*)[(R3 +X) sin’Bs, cosh 2a+ RZ, sin®Bs, sinh 2a],  ...... (4) 
K,=(1/|Z,)){(R3+X3) sin Bs, cos Bs, +XZy sin? By} (5) 


and Z, is the characteristic impedance of the Lecher wires. 
In order to simplify equation (1) the constant K, is made zero by adjusting 
the value of s,. The necessary condition is, from (5), 
(23+ Z,X, tan Bs,) sin Bs, cos Bs, =0, 
ine: tan'Bs;= = 752,45, i (6) 
the other solutions being of no practical importance. 


‘This critical value of s,, which we shall call s,, was determined experimentally ? 
in the original method. The impedance Z was removed and the lines short 
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circuited by a copper plate. In this case a and b should be zero and the equation 
(1) reduces to 


Re a Ae COseC“4as 2K COL bs, ees (7) 


K, and K; being unchanged and Ky simplified by the vanishing of the term 
with sinh 2a. 

If now K,=0, a graph of p? against s is symmetrical about its turning points. 
A finite value of K, introduces | 
an asymmetry and this, in the ; 
original paper, is the test for 
the critical setting. | 


aw 


For ka— UP argraphot p72? 
against cosec? fs is a straight 
line from which K,, which is 
necessary for the determina- . 
tion of a, is obtained. | | Ve 

It was noticed, however, \ | ‘ 


in our experiments, that to a \ | WE 


(p?, s) curve which appeared ake ii Oe 

; Se Bigs’ ¢ 
to be symmetrical there cor- See Lee 

2 are 0 20 40 | 2 3 4 
responded a (p?, cosec? fs) en ee 


graph which was a curve with 
two branches as shown in 
figure 2, i.e. Ky was still finite, and this 
suggested that there might be a more x 
sensitive test of the accuracy of the setting. 
It was thought of interest to examine the 
equation (7) more closely for values of s, not 
equal to s, in the hope of finding a quick 
and accurate way of determining s,. 

The form of the curve given by equation 
(7) is shown in figure 3 which is a graph 
of the simplified function 

p?=cosec? 0+ K cot 0 


from 6=0 to 0=7 for the arbitrary values 
K=1, K=%. For values of 0 greater than 
x (i.e. for increasing s) the curve is repeated. 

As @ increases the curve is traversed in 
the direction shown by the arrows if K is 
positive, and in the reverse direction if K Figure 3. 
is negative. Hence from the experimental 
curves it is quite obvious when the critical setting has been overshot. 
As s, approaches the value s,, i.e. as K tends to zero, the curves close in and 
approximate more and more nearly to the required straight line. 

A number of experimental curves were obtained for various values of s, 
and while the general behaviour corresponded to the above analysis in that the 
branches of the curves closed in and opened out again on passing through the 
critical setting, most of the curves showed unexpected peculiarities of curvature, 


26 


Figure 2. 
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and in some cases there was a crossing of the two branches. In no case could 
a reasonably good straight line be obtained. ; 
There is nothing in equation (7) to account for this, but a consideration of 
the general equation (1) shows that if a or b or both are finite, though small, 
the (p2, cosec? Bs) graph is considerably altered. 
Take first the case a=0, 6 finite; then 


p?=K,+ Ky, cosec? (b+ Bs) +2K; cot (b+ Bs). 


A graph of p? plotted against cosec?(b+ fs) would give the curve of figure a 
but experimentally the graph plotted is of p? against cosec?fs, and this shows | 
a crossing of the two branches at some point, since, at any particular value of p?, 
the value of Bs is too small. 

Figure 4 is a graph of the simplified function | 


p?=cosec?(6+10°)+K cot (8+ 10°) 


plotted against cosec?@ and should be compared with figure 3. 

The point of crossing over depends upon 
the value of b and comes nearer to the bend 
of the curve as 6 increases. ‘The width of 
the loop, however, that is the value of K, 
also affects the position of this point, so that 
it does not seem possible, as was hoped at 5} 
first, to get any information about a or 5b , 
from the form of this curve. 

Considering now the case a finite, b=0, 
which gives : 
p*(1 + sinh?a cosec? Bs) = K, + 

(K, sinh? a+ K,) cosec? Bs + 2K, cot Bs, 


the general effect of the finite value of acan 4 A 
be seen. ys 
Sa” 


left-hand side which is always greater than 
unity, the ordinates will be reduced by a 
fraction which increases with cosec?fs so 
that there will be a downward curvature 
superimposed upon both branches of the 
curve. At the same time there will be a change in the slope of the whole graph | 
due to the change in the constant which multiplies cosec? Bs on the right-hand side. |} 

Figure 5 shows the general curve, plotted for arbitrary values of the constants. 
with bothaandb finite. Figure 6 is one of the experimental curves for comparison. |} 

It was concluded then that the copper plate was not a sufficiently good short | 
circuit, and an estimate of the values of the constants a and b was made from the | 
nearest approach to a straight line. These values were not large. 

A copper plate of larger size was tried and gave a slight improvement, but. 
evidently the attachment of the plate to the lines was contributing to the reflection 
coefficient and this quantity could not be reduced to zero in our experiments. __ 

In measuring the reflection coefficient of an unknown impedance Z the 
quantity 5 is obtained from the displacement of the minimum of the (p?, s) curvel] 
when the short circuit is replaced by Z. To find aa straight line plot 1/(p?— K,) 


Owing to the multiplying factor on the oe 


= P 


4 5 


2 
cosec? @ 


Figure 4. 


A note on the Lecher wire method of measuring impedance 3.91 


against cosec? Bs is required and this can only be obtained at the critical separation. 
Hence, if the lines cannot be satisfactorily short circuited, the critical separation 
is not accurately obtainable, and there is a certain amount of error in the deter- 
mination of the reflection coefficient. The error is likely to be greater in a than 
in b, since the latter quantity is found more directly. Unfortunately it is just 
the quantity a which is the important one in the case of those dielectrics which 
show absorption. 
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Figure 5. Figure 6. 


The experimental work to test these points has been carried out by Miss E. B. 
Harriss and forms part of the thesis submitted by her for the M.Sc. degree of 
the University of London. ‘The author wishes to thank her for her help, and 
also to thank Professor Flint, with whom the problems have been discussed. 
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Slip at Grain Boundaries and Grain Growth in Metals 


By NN. EF eM OT TL, 
H. H. Wills Physical Laboratory, Bristol 


MS. received 6 Fanuary 1948 


si. INTRODUCTION 
N arecent paper from the Chicago Institute for the Study of Metals, Ting-Sui- 
Ke (1947) gives an expression for the rate of slip at the grain boundaries 
in pure polycrystalline aluminium; this is derived from measurements of the 
internal friction and some other anelastic phenomena. ‘The expression obtained 
for the velocity of slip at each boundary is 
Devote cCmsec, Theta (1) 
where a is the stress applied at the boundary in dynes/cm?, B = 34-5 kcal/gm.atom, 
and A ~18(cm/sec.)/(dyne/cm.). 
20-2 
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Ke also points out the following fact : let it be assumed for purposes of calcu- 
lation that between the grains there is a layer of fluid of thickness ¢ and viscosity 7- 
Then the velocity v of flow should be given by the formula 

v=ot/y. 1 a ee (2) 

Suppose then that we extrapolate the formula (1) deduced from experiment to 
the melting point of aluminium, and then compare the results with (2), taking for 
7 the viscosity of liquid aluminium at the melting point. We may deduce a value 
fort. This comes out to bea few angstroms only. This suggests, as is probable on 
other grounds, that the boundary is not a thick amorphous layer, but merely the 
surface of contact of two crystals fitting together as well as their different orienta- 
tions allow. 

In this note we shall put forward a theory of the slip process and, in particular, 
of the factors which determine the constants A and B. 


§2. SLIDING OF ONE SURFACE OVER ANOTHER 

Let us consider two close-packed surfaces of metal in contact and with the 
same orientation, and determine the activation energy required to slide one over 
the other. This quantity is required as a preliminary to the subsequent calcu- 
lations. 

If the distance between the planes is h, the stress o required to produce a 
displacement x is given by 

= Gril 

where G is the shear modulus for the direction of shear considered. ‘This expres- 
sion is valid only if x is small. For large x we may set as an approximation 


o = (Ga/27h) sin (27x/a), 


where a is the distance to the nearest stable position of equilibrium; the use of | 
the sine function gives an upper limit toc. Thus the activation energy per unit 
area is, to this approximation, 


i «dx =Gatj2nth. 


0 
If 2d is the lattice constant, then we may set 


d=a, b= Dano. 


and the number (1/w) of atoms per unit area is 1/,/3d?. Therefore the activation j} 
energy 7U required to slide m atoms of the top plane over the lower plane is 


nU = 3nGa?/4n?. 


For aluminium a=2-02 x 10-$cm., and since the metal is isotropic we may set for | 
G the shear modulus of the bulk material : 


G=2-5 x 10" dyne/cm. 
This gives U=0-066 ev. =c. 1:5 kcal. 


The smallness of this value is to be noted; it is much smaller than the heat of | 
vaporization (76 kcal.), and of the order of the heat of fusion (2:55 keal.). 


§3. THEORY OF INTERCRYSTAL MINES SLIP 


If two crystal planes are in contact, but cannot fit owing to different indices 
or orientation, one may suppose that the surface of contact is divided into islands 
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where the fit is reasonably good, separated by lines near which fit is bad. We may 
Suppose that on one of these islands the top plane can move relatively to the one 
underneath, and that this is the elementary act responsible for intercrystalline 
slip. With this hypothesis the rate of slip might be estimated as follows : 

Let each island contain n atoms and thus cover an area nw. In the absence of 
stress the activation energy for slip is nU. In the presence of a stress o it is 
nU + 30anw, it being supposed that 4a is the distance to the potential barrier 
separating one position of equilibrium from the next. Thus the number of times 
per second that each island moves a distance a in the direction of the stress is 

vexp { —(nU —onaw)/RT}, 
and in the opposite direction 
vexp { —(nU + onaw)/RT}, 
where v is the frequency of atomic vibration. It follows that the rate of slip is 
ve "UIRT sinh (onaw/RT). 
If o is not too great this becomes 
(va2nac/RT)e~ URE, 
If nU is to be identified with the observed 34-5 kcal., we must have 
W345) 1573. 
Oyitheve Wl sec ag 2 10° cm., w= /« 101° em? 7 — 600°, we get for A 
A =a" nok 65e10-*: 

In this treatment it has been assumed that o at the point of slip is equal to the 
applied stress. Some stress magnification is possible, in which case A would be 
somewhat greater. 

The objections to this treatment are : 

(a) The calculated value of A is too small by a factor 104-10°. 

(b) ‘The observed fact that at the melting point the slip is the same as that which 
would be given by a monomolecular layer of liquid appears in the theory as an 
accident. ‘The mechanism of flow in liquid aluminium cannot be anything like 
that sketched here, because the temperature dependence is 10—20 times smaller ; 
the viscosity of liquid metals depends on temperature according to the formula 

Boe, 
where W is of the order of the latent heat of fusion (Frenkel 1946). 

We must therefore modify our hypothesis. Let us suppose that the elementary 
act which allows slip to occur is the disordering of atoms round each island where 
fitisgood. The free energy F necessary to do this will approach zero at the melting 
point and vL at the absolute zero of temperature ; here L is the latent heat of fusion 
per atom. At any other temperature, let us assume F to be given by 

F=nL(1—-T/Ty), 
where 7',, is the temperature of melting. Let us assume also, since the disordering 
will result in a slip through a distance a, that a stress o will decrease or increase 
F by +3onwa. Then the rate of slip is now 
v= 2vaexp { —nL(1 —T/Ty)/RT} sinh (onwa/2k7T), 


which for small o reduces to 


_ vanwe nL “a (Fr) 
Crime Air) NRT | 
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For aluminium the latent heat L is 2:55 kcal/mol, so that m (the number of 
atoms to be disordered) should be about 14. 

The factor exp(nL/RTy) is exp (34.000/2 x 930) 108. One is thus able to 
explain the large factor outside the exponential. One also sees why the surface 
layer behaves like a layer of liquid at the melting point, since the activation energy 
for melting small volumes there disappears. 


§4. VELOCITY OF RECRYSTALCLIZATION 


Anderson and Mehl (1945) have shown that the velocity of growth of new 
crystals during recrystallization of cold-rolled aluminium is given by 


v=, e VF, 


where O&55kcal. and vg 10!4-10'6 cm/sec. Burgers and the present author 
(1947) have suggested that this very high value of uv, indicates that some process 
which depends on temperature through the factor e~°/"* triggers off the change of 
crystal form of a whole mosaic block, which then changes its crystal structure to 
that of the growing crystal without the intervention of any other thermally activated 
process. 

An alternative and perhaps more plausible explanation can be given along the 
lines of the last section. Let us assume that the crystal grows at each step by the 
disordering or melting of a group of n atoms of the cold-worked material which then 
takes up the structure of the growing crystal. We have at present no theoretical 
means of determining 7; as in the previous problem it must be the value for which 
the free energy of disordering, including surface energy, is a minimum. 

We should thus expect as in the last section a rate of crystal growth 


vaexp (O/RTy,)exp(—Q/RT)._--* ee (3) 


With va~10*?cm/sec., and the first exponential factor equal to 10”, this gives. 
1016e~ WRT Which is of the observed order of magnitude. 

This Aisi ak assumes that the local melting always results in an addition to the 
growing crystal and never the reverse process. ‘To test this assumption, let g | 
be the mean energy of cold work per atom. The condition for the validity of the | 
assumption is then ng>RT. Otherwise (3) should be multiplied by (1 — =e 
For heavily cold-worked aluminium according to Taylor and Quinney (1934) q 1s | 
30cals/mol. With n= 20 this gives ng/RT at 600°k. of the order 4. For lightly | 
rolled aluminium the value should be considerably less, but the order of magnitude | 
of (3) will not be affected. 
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EETTERS, TO CHE EDITOR 


Paramagnetic Resonance at Low Temperatures 
in Chromic Alum 


The use of centimetre wavelengths to determine the small splitting (due to the crystalline 
electric field) of the basic spin quadruplet of Cr+++ in potassium chromic alum has been 
described by Bagguley and Griffiths (1947), who find a value of 0:12 cm—! for the separation 
between the +1/2 and +3/2 levels. The same authors (unpublished) have obtained the 
value of 0-13; cm~ for the analogous splitting in ammonium chromic alum; independently 
Weiss, Whitmer, Torrey and Jen-Sen Hsiang (1947) give a value 0:15-+0-01 cm~! for this 
salt, obtained by a similar method. These values differ considerably from those obtained. 
at low temperatures by other methods: viz. 0:16-0:19 cm for potassium chromic alum 
(Bleaney 1939, Du Pré 1940, Casimir, de Haas and de Klerk 1939, Starr 1941, Broer 1947, 
Gorter et al. 1940, 1942) and 0-34 cm~ for ammonium chromic alum (Starr 1941). The 
application of centimetre wavelengths to the determination of the change of splitting with 
temperature was therefore undertaken by the authors. 

A fixed tuned cavity is used, resonant in the H, mode at about 3 cm. wavelength; it is 
excited by probes at the ends of small concentric line feeders. "The cavity is surrounded by 
a dewar vessel, containing liquid oxygen, liquid hydrogen etc.; its temperature is deter- 
mined by a vapour pressure thermometer, or by a thermocouple calibrated at various fixed 
points. 

The magnification factor Q of the empty resonator is about 13000 at 90° K. and 15 000 
at 20° K. The magnetic absorption in the crystal rises with 1/7, and it is therefore possible 
to obtain large effects from very small crystals. The crystal is placed in the resonator in a 
position of maximum magnetic field and minimum electric field (to avoid dielectric loss). 
The power transmitted through the cavity at resonance is observed by means of a silicon- 
tungsten rectifier while a variable magnetic field is applied, perpendicular to the R.F. magnetic 
field. ; 

As the temperature is lowered, the absorption curve for ammonium chromic alum remains 
similar to that observed at room temperature, except that the side peaks due to the Stark 
splitting move steadily inwards towards the main peak, showing that the splitting is de- 
creasing. At about 80° k., however, a transition sets in, which requires several minutes for 
its completion. It is accompanied by a considerable change in the resonant frequency of 
the cavity, and the absorption curve obtained on applying a magnetic field after the transition 
is complete is quite different from that obtained at higher temperatures. On warming up 
slowly, the reverse transition does not occur until about 100° k., after which the magnetic 
absorption curve reverts to that previously obtained at this temperature. These transitions 
are repeatable, but the exact transition temperature varies by a few degrees ; in general the 
hysteresis between the transition temperatures for cooling and warming decreases after the 
operation has been repeated a number of times. 

There is little doubt that this effect is the same as that found by Guillien (1939), who 
discovered a discontinuity in the dielectric constant, accompanied by jumps in the dielectric 
loss and specific heat, showing temperature hysteresis. The change in resonant frequency 
of the cavity observed by us is consistent in sign with this dielectric constant change, though 
comparatively small since the crystal projects only slightly into the R.F. electric field. ‘The 
effects are attributed by Guillien to “‘ co-operative relaxation of librational restraint on water 
or water-containing complexes with rising temperature’. If the water molecules invclved 
are those forming the octahedron around the Cr+*+* ion, the great influence of this transition 
on the Stark splitting of this ion can be understood, since Van Vleck (1939) has shown that 
the splitting can be explained by the slight distortion of the octahedron, together with other 
contributions due to the Jahn-Teller effect and to more distant ions. The transition has 
also been observed by Kraus and Nutting (1941), who found that the crystal often became 
practically opaque, as though shattered into an enormous number of tiny crystals; this is 
consistent with our observations. This shattering process does not, however, alter the 
magnetic absorption spectrum observed at higher temperatures. 
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Below the transition point, the magnetic absorption curve observed with the external 
magnetic field along the normals to the (111), (110), (100) crystal planes is quite Pane 
with that to be expected from the theory of Broer (1942), assuming a small trigonal fie 
superposed on the cubic field. Since the results at room temperature are 1n good agreement 
with this theory, one must assume that at low temperatures the crystalline electric field is 
differently aligned or has different symmetry. The quantitative interpretation of the 
magnetic absorption curve is therefore not possible at this stage. In this salt, however, there 
is considerable absorption at 3 cm. wavelength in zero magnetic field below the transition 
point, showing that the splitting is of this order. By examining the absorption curve in low 
fields, the splitting has been determined as 0°314 cm! just below the transition point, 
increasing slightly to 0-317 cm7? at 
20° x, ; the probable error in these 
values is less than 1%. Since, by 
Kramers’ theorem, an electric field 
cannot do more than decompose the 
spin quadruplet into two doublets, eet 
these values for the position of max1- 
mum absorption in zero field must 020) 
correspond to the Stark splitting, 
irrespective of the nature of the 
crystalline field. 

The variation of the Stark splitting 

down to 20° k. is shown in the figure O10 
for ammonium chromic alum. Above 
the transition point the splitting 
varies linearly with temperature, 
and is, presumably, associated with 
rotation of the water molecules. 
Below the transition point, the slight 
change of the splitting may be due 
to thermal contraction, and it seems 
unlikely that there will be any sub- 
stantial change below 20°K. The 
results of preliminary measurements 
on potassium chromic alum show 
that there is a similar decrease in the 
splitting on cooling downwards from room temperature, but no sudden transition occurs. 
Instead, abnormally weak side peaks appear in the absorption curve which, though not 
consistent with a trigonal field directed along the space (body) diagonal of unit cell, appear 
to correspond to a splitting of the same order as that obtained by other methods at low 
temperatures. Measurements are in progress at longer wavelengths to determine this 
splitting directly by a method similar to that used for the ammonium alum. 
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Variation of Stark splitting in ammonium chrome 
alum with temperature. 


(-] Bagguley and Griffiths. 
x Before shattering. 
© After shattering. 
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Milli-microsecond Pulse Generation by Electron Bunching 


Work is in progress in this laboratory on a project for coincidence counting of nuclear 
particles with a resolution of the order of 10~® seconds. For a study of transients of this 
‘duration, the first requirement is a technique of milli-microsecond pulse generation; 
these periods of time are below the limit of normal electrical pulse production methods. 
Equipment has been built for the production of these pulses by bunching electrons with 
alternating voltages of frequency 210 Mc/s., on the assumption that the mathematical 
theory of electron bunching applies (Webster 1939, Harrison 1947). 

The figure shows in diagram the electronic arrangement, which is in effect a klystron with 
the resonant catcher replaced 
by a non-resonant collector in 
the form of a coaxial line. 


210 Mds 
Oscillator 


Power from a 210-Mc/s. os- 

-cillator is fed into a coaxial line + 4 ea tee se 
resonator short-circuited at a Oe 
one end and capacitively term- oO i! co ave 
inated at the other by two pn [ 

grids 1 cm. apart, across which ee 

R.F. potentials up to 400 volts — °_| Bonehing 

-can be developed. Anelectron s=40Gn 


beam, after modulation by 


these grids, is magnetically 
focused on to the end of a se Rasen 4 


70-ohm coaxial line a distance aA at 

-s=40 cm. away. The outer oe 

conductor is terminated in a Bunching pulse generator and associated circuits in 
grid through which the beam block diagram form. 

passes before reaching the Inset: Theoretical pulse 


inner conductor 4mm. away. 

The bunched beam passing across the capacity between the beam and the conductor induces 
a voltage pulse which is then propagated down the line. The line is terminated in a probe, 
which feeds a length of 10-cm. waveguide through a “‘ door knob” transformer, and a 
milliammeter by which the D.c. beam current, of the order of 2 ma., is read. ‘The waveguide 
is terminated by a 10-cm. mixer with local oscillator, and a 60-Mc/s. I.F. amplifier with 
diode current meter (Collie, Hasted and Ritson 1948). ‘This serves as a narrow band power 
measuring device with which it is possible to pick out a harmonic component of the pulse. 
The peak alternating voltage in the bunching pulse generator, that is the peak energy in- 
crements given to the electrons, can be measured directly by back-biasing the electrode on 
which the beam impinges until no current is observed either with or without the alternating 
‘voltage at the buncher grids. 

The theory of electron bunching predicts that the mth harmonic of the alternating 
-component of the current is proportional to J, (nx); x is the bunching parameter, given by 
the relation x = V,s7f/V vp where vp is the velocity of the electrons accelerated to a potential 
V», and V, is the alternating voltage of frequency f. With the local oscillator tuned to a 
wavelength of 10 cm., maximum power in the 14th harmonic (of the order of 10~* watts) 
was obtained with values of V;=250-270 v. at Vp>=2000 v., and V,=140 v. at Vy =1350 v.; 
both of these results correspond to x=1:3, nx=16°8. The first maximum of the 14th 
Bessel function occurs at nx=16, so that the agreement with theory is satisfactory. ‘The 
theoretical pulse, the shape of which it is hoped to verify with viewing equipment, is shown in 
the inset to figure 1 to be of the order of 2 10~!° seconds long. Responses have been 
plotted for varying values of nx, and lie on a curve of roughly the same form as that of the 
Bessel function, although observation of maxima beyond the first has not yet been attempted. 

From the power observed in the 14th harmonic, and the agreement with theory, we 
-should expect that this device will be capable of supplying pulses of variable width and of 
.duration short enough for testing wide-band coincidence counting circuits. 
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REVIEWS OF BOOKS 


Radar System Engineering, edited by Louis N. Ripenour. Pp. xviiit+748: 
with 646 figures. (New York: McGraw Hill Book Co. Inc., 1947.) 
375. Ody nets 


During World War II a tremendous research and development effort went into the 
development of radar and related techniques, and resulted not only in hundreds of radar sets. 
being available for military use but also in the collection of much valuable data and new 
techniques in the electronics and high-frequency fields. This effort was a splendid example 
of co-operative research under war-time conditions, for each new advance was rarely the 
work of an isolated worker, but generally the combined efforts of a very large group of 
people: British, American, Army, Navy, R.A.F., civilian, Government, and industry, 
working in different laboratories, and often at widely separated locations. A very large and 
important part of this combined effort was undertaken by the Radiation Laboratory of the 
Massachusetts Institute of Technology, and they have now played a further important part 
by undertaking the heavy responsibility of preparing for publication the results of this 
collective work. As originally planned the series comprised twenty-three volumes dealing 
with theory and techniques, but later it was decided to extend the series by the addition of a 
* further five books dealing with radar and applied systems. 

The present book, which is the first of the series to be published, is also the first of the 
five dealing with radar and applied systems. It is intended to serve as a general treatise 
and reference book on the design of radar systems : it deals, in fact, with the principles and 
practice of radar. Much of the work has been written by Professor Ridenour himself, but- 
he has been assisted in this monumental task by many other writers, thirty-two of whom 
have contributed to this first volume. 

‘The aim of the book, as stated in the preface, is to cover the collective efforts of the many 
workers, both in the United States and in Britain, Canada, Australia and other British 
Dominions, but in spite of this, the book tends to be rather overweighted with American 
applications of the new techniques, and this may detract from its value to the British reader.. 
Furthermore, although various references to the original literature are given throughout the 
book, the references mentioned are usually to obscure reports (e.g. M.I.T. Radiation 
Laboratory Reports) not generally available in this country. The book would have been of 
more value to the British reader if reference could have been made to the excellent literature 
on the subject published-last year by our Institution of Electrical Engineers. However, 
these are small points, and it can be said at once that this is a very good book, which should 
be of real value to all those wanting to make themselves conversant with the principles and 
practice of the subject. 

Chapter I is a general introduction and explains clearly and concisely ‘“‘ what radar does ’” 
and “ how it works’. It then goes on to give details of the components of a radar system, 
the types of radar developed during World War II and the performance obtained. The 
author departs a little from the usual custom by referring at the outset in his explanation to a. 
microwave pulse radar system, rather than to the simpler and historically earlier metric-- 
wave systems. ‘This insistence on the importance of microwaves from the start is quite 
reasonable. Considerably more work has gone into pulse radar than into any other kind,. 
and by far the greater part of this work has been concerned with microwave pulse radar. 
Moreover, it is now recognized that for almost all radar purposes microwaves offer advan-- 
tages. ‘This first chapter is very well written and the explanations are good, but one small. 
mis-statement has been made. On page 11, in referring to super-refraction conditions, the 
author states “that such phenomena are relatively rare and essentially unpredictable ”’.. 
This statement would be true in referring to this country ; it is probably also true in refer- 
ring to conditions in North America, but it is obviously not true in relation to conditions in 
more tropical countries. The Royal Air Force have observed super-refraction conditions 
repeated for several weeks in the year, for several years in succession, in certain parts of 
India, and so for certain parts of the world, at least, the increase in range of a radar set 
beyond the optical horizon is a predictable phenomenon. 
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Chapter II contains a treatment of the radar equation according to which the range (R): 
of a radar set is a function of the transmitted power (P), the gain of the antenna (G), the 
wavelength (A), the radar cross-section of the target (c) and the minimum detectable received 
signal (S)._ The author writes down his formula as 


s= IRE to GH 
PERRO. Uo - 


and proceeds to discuss in a very thorough fashion the consequences of the equation, and 
the choice of the various parameters for particular purposes. He discusses beams of special 
shapes, the limitations of the receiver, receiver noise and the statistical problem, and the 
modification of the range equation called for when “‘ free-space ’’ propagation conditions 
are not applicable. He discusses one case in detail, namely the calculation of the maximum 
radar range of an aerial target when the propagating path is over a flat, perfectly reflecting 
surface, and shows that in this case the ‘‘ free-space’? formula becomes 
Gna 
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and he discusses, in detail, how accurately the flat, perfectly reflecting conditions laid down 
in his analysis are fulfilled in practice, and therefore how reasonable it is to use the new 
equation. ‘The treatment, so far as aerial targets are concerned, is thorough and excellent, 
but the author fails to point out the limitation of his analysis to ship targets, and except for a 
casual reference to figure 2.12, the reader is left with the impression that equation (B) may 
be used for ship targets. This is very unsatisfactory. "The reviewer would have liked to 
have seen the important case of the ship target considered in detail with proper comment on 
the special conditions which arise in this case (namely that the target is not placed in a 
constant incident field, but in one which varies across the surface of the target), and the 
modifications to the concept of target cross-section made necessary. A satisfactory treat- 
ment of this problem was published recently by Wilkes and Ramsay (Proc. Camb. Phil. Soc..,. 
1946, p. 220). 

Chapter III gives a detailed treatment of the properties of radar targets, and after a 
consideration of targets of simple shapes, e.g. sphere, corner reflector etc., some discussion 
is included on methods of target shaping and the use of absorbent materials to diminish 
cross-section. The chapter then goes on to discuss the radar returns from more compli-. 
cated targets such as aircraft, ‘“‘ window ”’, rain, land, and the sea, and this is illustrated with 
many excellent photographs of the radar returns actually obtained. ‘The use of radar 
*“ land-painting ”’ is considered very fully in connection with the important applications of 
this technique to navigation problems. 

Chapter IV discusses the limitations of pulse radar, and Chapter V gives a detailed 
treatment of the many c.w. radar systems, which includes many critical inter-comparisons: 
of the performance of c.w. and pulse systems. The c.w. systems considered include 
simple Doppler systems, range-measuring Doppler systems, frequency-modulation range 
measuring systems, multiple-target frequency-modulation systems, pulse-modulated! 
Doppler systems etc. One small mis-statement occurs sufficiently frequently to be 
annoying, namely that a frequency-modulation system will work down to zero range. ‘There 
is a minimum range for a F.M. system but it is not zero, but equal to cdf/2F Af, where df is 
the beat frequency, F is the modulation frequency and Af is the frequency swing. This 
reduces to c/2Af when 8f=F. The reviewer would have liked to have seen this point dealt 
with properly. The corresponding figure for a pulse system is 3c7, where 7 is the pulse 
length. For a normal F.M. system the limitation is usually of the order of 1 metre, whereas 
for a pulse system the limitation in range might be as great as 100 metres ; hence all the 
author’s statements about the advantages of c.w. systems for the accurate measurement of 
very short ranges are obviously true. 

Chapter VI gives details of the presentation of radar information, and most readers will 
find this a valuable introduction to the various applications of radar methods. Chapter VIT 
considers the employment of radar data, and gives examples of radar organization, viz. radar 
in the R.A.F. Fighter Command and the U.S. Tactical Air Commands. It is very pleasing 
to the reviewer to find this important material included. The technical design of a radar set 
is relatively straightforward to the specialist at least, but the employment of radar data in the: 
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best possible way is far less well understood, and Professor Ridenour is quite correct when 
he states that the major improvements to be looked for in the use of radar during the next 
few years will lie, for the most part, not in technical radar design, but in the field of fitting the 
entire radar system, including its operational organization, to the detailed needs of the use 
and the user. This obviously applies to both military and peace-time applications. 

Chapter VIII continues with a consideration of radar beacons, and gives a good intro- 
duction to the subject. The whole of volume 3 of the series is also devoted to this subject. 
Chapters IX to XIV take up the leading design considerations for the various important 
components that make up a radar set, and deal respectively with aerials, magnetrons, R.F. 
components, radar receivers, indicators and power supplies. The treatment is thorough and 
written throughout by specialists in their subjects. The reviewer would have liked the needs 
.of the British reader to have been considered, by finding a few more references to British 
practice, and he would have expected the 70-paze chapter on the magnetron developments to 
have mentioned that the all-important cavity resonator magnetron was invented by Professor 
Randall, our recent Duddell medallist, but bearing in mind the emphasis on co-operative 
-efforts, perhaps this is not so unfair. Chapter XV is one of the most valuable in the book and 
gives some examples of radar system design practice. The authors (R. G. Heib and R. L. 
Sinsheimer) consider the complete design problem from a consideration of the formulation 
.of requirements, to the choice of the various radar parameters, and to the detailed problems 
of the design of the various component parts of the equipment. ‘This chapter has been 
-carefully written and will repay careful study. Finally the book closes with a consideration 
of two new and important ancillary techniques—moving target indication and the trans- 
mission of radar displays to a remote indication by radio means. 

The book is well produced and illustrated and can be recommended to any serious 
student of radio. It is hoped that some of the criticisms made above will be dealt with in 
future editions. DENIS TAYLOR. 


Radioactivity and Nuclear Physics, by J. M. Cork. Pp. x+313. (New York: 
Van Nostrand Co, Inc., 1947.) 


This book will be easily read—and with a fair measure of assent—by anyone familiar 
with the subject, yet it is a tantalizing and a disappointing book. It lacks logical form : 
there is a great deal that passes for information, very little which seeks to be educative. 
The serious student will not find that the book answers his questions, though it will certainly 
give rise to many, in a bewildering sequence. All the fashionable terms are freely used— 
-curium is mentioned on p. 3, fission on p. 6, the relativistic change of mass with velocity on 
p. 14, and the meson and the neutrino on the following page. But that surely is the wrong 
place for them, unless a very casual acquaintance with words is all that the reader requires. 

“It has been the aim of this text to offer an authentic historical development of each 
topic ”’: thus the preface. This much must be granted, that the aim is evident throughout; 
nor is its success negligible (most obvious perhaps in Chapter 3 which describes the evolution 
-of large-scale apparatus for particle acceleration), but the historian is at fault, sometimes, and 
at others is uncritical. It would seem strange that the earliest reference to Rutherford’s 
work in radioactivity is of date 1902, and the almost literal statement that the experiment of 
Ellis and Wooster (1927) was planned in order to test the neutrino hypothesis is a peculiar 
reversal of history. Yet the reader should be grateful that the historical approach is followed 
consistently enough, whenever occasion offers, and the copious references which appear as 
footnotes to the text are a good feature of the book. 

At the end appendices give tables of isotopes (stable and unstable), of fundamental 
constants, of numerical data concerning electrons, protons and a-particles over a considerable 
range of energy, and of the characteristic x-ray excitation energies of the elements titanium 
to uranium. ‘These tables also are good (incidentally the table of isotopes includes some 
data regarding recently discovered isotopes of uranium and the transuranics which, so far 
as the present reviewer is aware, have not hitherto been published either in this country or in 
America). But why exclude rhenium and bismuth from the table of X-ray excitation 
energies? This is rather a frivolous question, but it should set the user of the tables on his 
guard—for there are, perhaps inevitably, some errors as well as omissions. 

From what has been said in praise of the industry which has gone to the compilation of 
lists of references and useful tables of data it will be evident that it is the greater pity that — 
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the general level of the text is not better. But criticism cannot be avoided : the text too- 
frequently appeals to the memory of the student, not to his understanding. And even his 
memory is not encouraged by a repetitiveness which so often introduces inconsistency : 
the maximum energy of the f-particles of radium E is given as 1050 kev. on p. 124 and as 
1:17 Mev. on p. 133, the discovery of the positron is attributed to Anderson in 1933 on patos 
—and in almost identical language the date is given as 1932 on p. 124 (incidentally, on the 
former page the relativistic theory of the electron is dated 1934—six years too late), the 
energy release in the reaction ' ?B (n, a) 7Li is quoted as 2:5 Mev. on p. 177 and as 3-0 Mev. 
on p. 181, that for the reaction 7H (d, n) He as 3-17 Mev. on p. 178 and 3-16 Mev. on p. 198 
(the “ best” value is in fact 3:23 Mev.), that for °Be (d, 7) 1B as 4:1 Mev. on p. 178 and 
3°74 Mev. on p. 198, and on p. 210 it is stated that the specific ionization due to an electron 
is independent of its energy when’the energy is high, whereas on p. 229 we have the correct 
statement that the ionization increases with increasing energy in this range. These are pure 
inconsistencies, which diminish confidence when they are detected; of the outright errors, 
which the unwary student might not find it so easy to detect, it will suffice to quote one in full 
(p. 125)—** The upper limit of the spectrum for the ®*Cu positrons is 0:66 Mev. compared 
with 0-58 Mev. for the negatives. The difference may be due to the action of the positive 
nucleus on the departing particles, being an attraction for the negatron and a repulsion 
of the positron”. It must be admitted that there is a superficial air of verisimilitude in 
this wholly misleading statement. N. FEATHER. 


Tables of Spherical Bessel Functions. Vol. II. Prepared by the Mathematical 
Tables Project of the National Bureau of Standards. Pp. xx+328. (New 
York : Columbia University Press, 1947.) $7.50. 


The tabulation of the functions 1/(7/2x) . J (x) is continued *, and this volume covers 
the ranges 
+2v=29 (2) 43; x=0 (0-01) 10 (0-1) 25; 
+2v=45 (2) 61; KON O=1)e2 5 
Apart from a few entries near the zeros, the values are given to eight or more significant 
figures for x<<10, and to seven for x>10. Second differences, modified second differences, 
or second and fourth differences, are given to facilitate interpolation, except in ranges from 
x—=0 where the functions vary very rapidly. To enable interpolation in these latter ranges. 
(and also in the v direction) to be performed, values of 
(oe) 2 eo aT ( vp i) J (x) 
are given for 
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Dy=1 O61; x=10 (0-1) 25; 
y= 29 (2) 33" x=0 (0-1) 9-5 (0-05) 10 (0-1) 25; 
2 y= 35 (2).61; x=0 (0-1) 25; 


the first group to nine decimal places, the second to (mostly) seven significant figures. 
Finally the zeros of J,(x) and J,’(x) obtained by inverse interpolation from the main 
tables of this and the preceding volume, and the value of each function at the zero of the 
other, are given to six (occasionally more) decimal places. 
The functions here tabulated are important in, for instance, solutions of the wave 
equation in spherical polar coordinates, and are to be welcomed as extending the range of 
the tables already given in volume I. WwW. G.B. 


Table of the Bessel Functions J((z) and J,(z) for Complex Arguments. Second 
edition. Prepared by the Mathematical Tables Project of the National 
Bureau of Standards. Pp. xliv+403. (New York: Columbia University 
Press, 1947.) $5. 


This is substantially a reprint of the first edition, but some errors of lettering in the 
diagrams have been corrected and explicit relations between the functions tabulated along 
the 45° ray and the ber and bei functions are now included, with other minor revisions, in the 
introduction. W. G. B. 

* cf. Proc. Phys. Soc., 1947, 59, 259. 
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Time and Thermodynamics, by A. R. UBBELOHDE. Pp. vi+110. (Oxford: 
University Press, 1947.) 6s. net. 


The main theme of this well-produced little book is that, in the metrical study of the 
external world, vision (i.e. the use of any of the so-called ethereal vibrations) is of prime 
importance, and since vision depends on the absence of radiation equilibrium (in a constant 
temperature enclosure the world would appear as an undifferentiated uniformity) the 
conditions of such absence (e.g. the continual increase of entropy) are ultimately the con- 
ditions of obtaining metrical knowledge of the universe. ‘‘ And since this is so it may be 
anticipated that the trend of events, or direction in time, which is indicated by thermo- 
dynamics has its foundations primarily in the thermodynamics of radiation ”’ (p. 77). 

The book is disappointing, both in its main thesis and in the development of it. "There 
are other conditions, besides absence of radiation equilibrium, of obtaining metrical know- 
ledge ; for example, the temperature of the body must not exceed, say, 110° F., the atmo- 
sphere must not be deprived of oxygen, and so on. It is not clear why, among 
these conditions, the absence of radiation equilibrium is ‘‘ prime’. When we come to time 
(which, by the way, is very much the junior partner, notwithstanding the order of names 
in the title), the book is even less satisfactory. It appears to be held that entropy enables us 
to measure spontaneous change. ‘“‘ Historically, the way in which the concept of entropy 
first aroused an interest outside its immediate technical applications was by providing a 
quantitative measure for the trend of spontaneous happenings in the particular instance of 
the direction of spontaneous heat flow from one body to another ” (p. 78). This, of course, 
is nonsense. ‘The second law of thermodynamics merely requires that in an irreversible 
change from one equilibrium state to another the entropy must increase. It says nothing at 
all about the rate of increase, and therefore provides no quantitative measure of the trend. 

Still more unexpectedly, the book abounds in gross errors. Perhaps the most sur- 
prising is in the discussion (pp. 83-6) of the expansion of a gas from a volume v, to a volume 
Uy. Itis stated that “ the ratio of the probability W, that the molecules will all remain in the 
volume v, to the probability W, of their occupying uniformly the larger volume v, which is 
made available to them, is given by W,/W,=v,/v,’’, and even the recital of an example of 
patrons choosing theatre seats, in which the absurdity is equally obvious, has not opened 
Professor Ubbelohde’s eyes to it. Again, “ the unfolding of the solar system on Laplace’s 
dynamical theory would be completely reversed by reversing the momenta of all the masses 
taking part’ (p. 79). Does Professor Ubbelohde really think that if the bodies revolved the 
other way the planets would return to a Sun which would expand back to the original 
nebula? Again, infra-red solar radiation “‘ is largely responsible for the sun’s heat, and 
for the comfortable use of glasshouses ” (p. 61). Professor Ubbelohde surely knows how 
small a proportion of the sun’s heat comes from infra-red radiation, and how negligible its 
part is in the greenhouse effect. Once more, the uncertainty principle is stated as 
AE’. At >h, “‘ where A is Planck’s constant in radiation theory, and AE is the minimum 
uncertainty in F, At the minimum uncertainty in t”’ (p. 75). These increments are, of 
‘course, any uncertainties at all. The least value for each is zero and the greatest infinity, 
and only the product has a minimum value. There are other faults of the same type, of 
which we need add only the inexcusable transformation by which probability conceptions, 


first properly introduced (p. 83) as terms in which the laws of thermodynamics can be 1 


ce 


“reformulated ”’, automatically become “‘ the statistical foundations of the Second Law of 
“Thermodynamics ”’ (pp. 87, 97, 101 ; reviewer’s italics). 

It is difficult to persuade oneself that this book has really come from so competent and 
well-informed an author as Professor Ubbelohde, but one must believe the Oxford Univer- 
sity Press. He is plainly guilty of contempt of subject. Fortunately a reviewer is in the 
position of jury and not of judge, and we gladly refrain from passing sentence. 

HERBERT DINGLE. 


Fundamentals of Photography with Laboratory Experiments, by Pau E. Boucuer. 
2nd edition. Pp. xii+394+liii. (New York: D. Van Nostrand Co. Inc., 
1947) é 


The first edition of this work was published in 1940 ; the present edition has been 
-broadened, so the preface tells us, by the inclusion of more material on colour photography 
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‘and by chapters on ultra-violet and infra-red, x-ray, and motion-picture photography. 
The book is intended to supply lecture and laboratory notes for a one-semester course in 
photography for college students. The major part consists of some nineteen chapters 
ranging from the historical development of photography to cinematography. There follow 
a series of twenty-five laboratory exercises which follow the same order as the chapters. 

The book appears far too ambitious. It contains too much and as a result in the treatment 
‘of specific matters it is superficial, and many careless statements occur. The sketch of the 
historical development of photography is very unsatisfacotry. Fox Talbot is left out 
entirely and Herschel is not mentioned. Niepce is mentioned, but his really important 
early work is overlooked. In the next chapter is a heading—‘‘ The Reduction of the Latent 
Image’. One knows something of the development of the latent image, but we do not 
usually talk about its reduction. Very many more examples could be cited of careless 
statements, and it has to be remembered that this work is intended as a college text-book. 
It does not deal satisfactorily with the fundamentals of photography, while at the same time 
dealing very superficially with many other matters. It is surely unnecessary to provide a 
-description and a wiring diagram of an X-ray generator when no mention is made of the use 
of intensifying screens in x-ray work. There is no mention of stereoscopic photography, 
but quite a lot about certain aspects of ultra-violet and infra-red that might well have been 
left out. 

The general plan of the book is good and most of the information sound and well 
arranged ; it is a great pity that it should be marred by such imperfections as are detailed 
above. H. W. GREENWOOD. 


Atomic Energy, by R. R. Nimmo. Pp. 201. (London: The Pilot Press Ltd., 
1947.) 9s. 6d. 


“The average citizen cannot be expected to understand clearly how an atomic bomb .. . 
works ”’ says Professor Smyth in the preface to the official United States booklet on Atomic 
Energy. On the other hand, the Pilot Press and their editor, Dr. Lovell, believe that the 
““ non-specialist public” is capable of keeping abreast of such topics as atomic energy, and 
they publish this book by Dr. Nimmo to prove their point. We shall not attempt an opinion 
on this controversial question here, but simply examine the content of Dr. Nimmo’s book. 
It contains an excellent account of nuclear fission, chain reactions and uranium piles, 
leading up to the atomic bomb and atomic energy plants. A fairly conservative attitude is 
indicated towards the future of atomic power plants. ‘Ten thousand tons, or more, of 
uranium would be needed yearly to provide the power requirements of Great Britain, and 
there is little up-to-date information available, apparently, regarding the size of uranium 
deposits. Generally speaking, all the information on atomic energy which has been 
released, and all the deductions which a nuclear physicist can make from the official state- 
ments, are included in the book. It can be recommended to any scientist as a competent 
survey, carefully written by an expert in this field. The Pilot Press turn out a pleasing 
volume and are welcomed into this sphere of publication. ‘Their production would have 
been improved, however, if they had conformed to the standard practice of printing mathe- 
matical symbols in italics. W. COCHRANE. 


The Evolution of Modern Physics, by C.'T. Cuase. Pp.ix+203. (New York: 
D. van Nostrand Co. Inc., 1947.) 14s. 


There is always much pleasure to be derived from accounts of the many striking incidents 
and events in the history of physics. Galileo’s work on falling bodies ; Fresnel’s triumph 
over the academic doubters ; Réntgen’s discovery of x-rays ; Becquerel’s discovery of 
radioactivity—these, and many others, never fail to delight. If the writer manages, at the 
same time, to show the use of the scientific method on these occasions, it will benefit the 
seneral reader, who is patently aware of the achievements of science, but less familiar with 
the discipline of thought which science requires. : 

Professor Chase’s book contains a well-written and connected account of the history of 
physics, with a reasonable emphasis on the methods and thought involved, and with practi- 
cally no mathematics. Some points dealing with the atom and nucleus are rather concen- 
trated and one looks in vain for the apt analogy that would help over an awkward fence. 
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There are no diagrams or sketches in the book but there are a dozen plates, mainly. of. 
complicated pieces of apparatus not referred to in the text. The book is, however, a useful 
one, since it maintains a high standard. It will help the reputation of physics and will lead 
people to an understanding of the physicist’s work. W. COCHRANE. 


The Electron Microscope, by V. E. Cossterr. Pp. viiit+128. (London: Sigma 
Books Ltd., 1947.) 7s. 6d. 


The electron microscope is hardly more than 15 years old, yet half a dozen books have 
been written on it. These books range from the exhaustive monograph to the short 
introduction for scientists who have no intention to specialize in the subject. The author 
of the present book now brings principles and uses of the electron microscope “ before that 
wide and increasing public which has a lively interest in the advance of science ”’. 

The content of the book may be divided into three parts of roughly 40 pages each = 
(1) General introduction, (2) Principles and designs of electron microscopes, (3) Present use 
and future possibilities of electron microscopes. 

The treatment starts in a very general way with some principles of vision, with refraction, 
prisms and lenses ; then it proceeds to the eye, to telescopes and microscopes. In rather 
quick succession follow explanations on the electron, its refraction and deflection, and on 
general features of electron lenses. Two excellent chapters on electron microscopes contain 
details of magnetic and electrostatic types, of shadow- and of scanning-microscopes. ‘Then 
follows a brief survey of electron microscopical methods and preparation techniques for the 
investigation of powders, fibres, bacteria, viruses etc. Sectioning and straining methods 
are discussed and replica methods for the investigation of surfaces are described. A con- 
cluding chapter contains information on limitations and possibilities of the electron micro- 
scope. ‘These points are of great interest, but perhaps the author has here exceeded the 
frame of a booklet which is intended for a generally interested public. 

The little volume is well illustrated by 40 figures in the text and 12 plates containing a 
selection of beautiful electron microscopical photographs. The name of the author— 
who has also written a standard monograph on Electron Optics and who is now in charge of 
electron microscopic developments in the Cavendish Laboratory in Cambridge—vouches 
for a masterly treatment of the subject. The book should be warmly recommended to 
everybody who is interested in electron microscopy. O. KLEMPERER. 


A University Text-Book of Physics, Vol. I: Properties of Matter, by PoYNTING 
and ‘THOMSON and Topp. Pp. 278. (London: Charles Griffin and Co., 
Lids, 19472)" #205: 


Of the famous text-book by Poynting and Thomson, the most famous volume was 
the first, on properties of matter. It is appropriate, therefore, that the first volume to 
appear of the revision now being undertaken should be that on properties of matter. The 
original work is so well planned, and so good in its treatment, that one feels reluctant 
to see too much revision; if an author desires to treat the whole subject differently, it is 
better that he should write a completely new book. Fortunately, Prof. Todd, the reviser 
of this volume, and an old student of both senior authors, has treated the original text 
with due respect, and has revised only to the extent that was called for by the developments | 
in the interval since first publication. His chief additions are in the chapters on gravitation, 
surface tension, and viscosity, and he has certainly made the book a suitable introduction 
for present-day students. Like others of Griffin’s post-war publications, the. book is | 
printed on thinner paper than they used to adopt, and the result is to make a pleasanter | 
book to handle. mts S. 


Frontiers of Astronomy, by Davin T. Evans. Pp. 175. (London: Sigma 
Introductions to Science, 1946.) 6s. net. 


This is a very good little sketch of the present state of astronomy, giving concise information 
concerning its chief problems, methods and results. The book is written ina racy, conversa- 
tional style and is illustrated by a number of diagrams in keeping therewith. A remarkable | 
amount of knowledge is compressed into a small space without introducing obscurity, and |] 
although there are, unfortunately, a rather large number of slips, they are mostly of slight 
importance, and the reader can be satisfied that the impression he gets is an essentially true 
one. As present prices go the bock is strikingly good value for the money. H.D. 
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in 1947 at 24s., in paper covers. Present price the same. 

Catalogue of the 31st Exhibition of Scientific Instruments and Apparatus. Pp. 298. Published in 
1947 at 6s.,in paper covers. Present price 3s. 6d. 

Catalogue of the 30th Exhibition of Scientific Instruments and Apparatus. Pp. 288. Published 
and reprinted in 1946 at 3s., in paper covers. Present price 2s. 

Report on Defective Colour Vision in Industry, by a Committee of the Colour Group. Pp. 52. 
Published in 1946 at 3s. 6d., in paper covers. Present price the same. 

Science and Human Welfare. A Conference held by the Association of Scientific Workers, The 
Physical Society and other bodies. Pp. 71. Published in 1946 at 2s. 6d., in paper covers. 
Present price 1s. 6d. 

Report on the Teaching of Geometrical Optics. Pp. 86. Published in 1934 at 6s. 3d.,in paper 
covers. Present price the same. 

Report on Band Spectra of Diatomic Molecules, by W. Jevons, D.Sc., Ph.D., Pp. 308. Published 
in 1932 at 18s., in paper covers, and 21s., bound in cloth. Present prices 9s. in paper covers 
and 18s. bound in cloth. 


Discussion on Vision. Pp. 327. Published in 1932 at 12s. 6d., in paper covers. Present price - 


6s. 6d. 


tian on Audition. Pp. 151. Published in 1931 at 7s. 3d., in paper covers. Present price 
S. 


Discussion on Photo-eleciric Cells and their Applications. Pp. 236, Published in 1930 at 12s. 6d., 
in paper covers. Present price 6s. 6d. 


A Discussion on the Teaching of Geometrical Optics. Pp. 84. Published in 1929 at 4s. 6d., in paper 
covers. Present price 2s. 6d. 
The Decimal Bibliographical Classification (Optics,Light and Cognate Subjects), by Professor A. F. C. 


PoLiarD, D.Sc. Pp. 109. Published by the Optical Society in 1926 at 7s. 6d., bound in 
cloth. Present price 4s. 


Motor Headlights. Pp. 39. Published by the Optical Society in 1922 at 2s. 6d., in paper covers. 
Present price 1s. 6d. 


Report on Series in Line Spectra, by (the late) Professor A. Fowirr, C.B.E, Sc.D., F.R.S. 
Pp. 182, with 5 plates. Published in 1922 at 12s. 6d., in paper covers. Present price 6s. 6d. 


A Discussion on the Making of Reflecting Surfaces. Pp. 44. Published in 1920 at 5s., in paper 


covers. Present price 2s. 6d. 


Reports on Progress in Physics. Vols. IV (1937, reprinted 1946) and X (1944-45), Bound in cloth. 
Price 30s. each. 


The Proceedings of the Physical Society. Vols. 1 (1874-75) -— 59 (1947), excepting a few parts which 
are out of print. Prices on application. 


The Transactions of the Optical Society. Vols. 1 (1899-1900) - 33 (1931-32), excepting a few parts 
which are out of print. Prices upon application. 


Orders, accompanied by remittances, should be sent to 
THE PHYSICAL SOCIETY | 
i Lowther Gardens, Prince Consort Road, London S.W.7 


